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1.  INTRODUCTION 


During  the  period  1  May  1980  through  30  April  1981,  the  Electronics  Research 
Laboratory  (ERL),  University  of  California,  Berkeley,  is  receiving  support 
for  its  basic  program  from  the  Joint  Services  Electronics  Program  ( JSEP) , 
sponsored  by  the  Departments  of  the  Air  Force  (Air  Force  Office  of  Scientific 
Research),  Army(Array  Research  Office),  and  Navy  (Office  of  Naval  Research), 
under  contract  F49620-79-C-0178. 


This  annual  report  is  an  updating  report;  hence  it  covers  the  period  1  Sept¬ 
ember  1979  through  31  August  1980.  Significant  Accomplishments,  where  applic¬ 
able,  appear  in  the  project  reports  (section  A).  In  particular,  the  three 
most  significant  accomplishments,  as  determined  by  the  Director,  are  sunmar- 
ized  as: 


i.  Project  EM-81-1  -  Multiple  Scattering  of  Conducting  Bodies 

We  have  positively  demonstrated  that  the  Unimoment  Method  we  have 
developed  is  adaptable  to  multiple  scattering  of  arbitrary  shaped 
scatterers.  The  computation  is  also  proven  to  be  accurate  and  e- 
conomical.  The  results  are  contained  in  a  PhD  thesis  submitted  by 
Dr.  John  F.  Hunka,  entitled,  "Electromagnet ic  scattering  by  two 
bodies  of  revolution.*1 

i 

ii.  Project  SSD-81-1  -  Advanced  Lithographic  Techniques  for  Microcircuits 

A  breakthrough  in  the  problem  of  negative  resist  swelling  and  its  in¬ 
herent  electron  beam  resolution  has  been  made.  We  have  demonstrated 
that  a  non-swelling  positive  photoresist  (AZ1350J)  may  be  used  as  a 
negative  electron  resist  by  means  of  image  reversal.  Resolution  be¬ 
low  1  ym  in  films  1.2  \m  thick  was  observed  at  a  dose  of  25  yC/cnr  at 
20keV.  The  high  etch  resistance  in  plasma  etching  is  also  retained. 


iii.  Project  ISS-81-3  -  Computer  Optimization  of  Electronic  Circuits 

The  significant  accomplishments  of  our  research  can  be  summarized  as 
follows:  a)  We  have  developed  a  prototype  interactive,  optimization 

based  computer  aided  design  package,  capable  of  accepting  semi- infin¬ 
ite  constraints  on  frequency  and  time  responses,  but  not  tolerancing 
or  post  manufacturing  tuning.  This  package  serves  as  a  test  bed  for 
the  development  of  ideas  and  concepts  for  a  much  more  sophisticated 
package  which  is  now  being  constructed,  b)  We  have  shown  that  de¬ 
vice  models  can  be  obtained  which  make  it  practical  to  solve  the  op¬ 
timal  design  centering,  tolerancing  and  tuning  problem,  c)  We  have 
developed  a  number  of  optimization  and  inequality  solving  algorithms 
for  computer  aided  design  which  either  solve  problems  that  were  hith¬ 
erto  unsolvable,  or  which  converge  much  more  rapidly  than  currently 
available  algorithms. 

Further:  the  paper  by  D.  J.  Allstot,  R.  W.  Brodersen  and  P.  R.  Gray  entitled, 

*'M0S  switched  capacitor  ladder  filters,*1  IEEE  Jour.  Solid-State  Circuits, 
vol.  SC-13,  no.  6,  Dec.  1978,  pp.  806-814,  received  the  1980  IEEE  W.  R.  G. 
Baker  Prize  Paper  award.  The  work  reported  was  supported  in  part  by  JSEP. 


INTRODUCTION  -  continued 


The  report  contains: 

A.  1980-1981  JSEP  Projects:  as  of  31  August  1980 

B.  Consolidated  List  of  Published  JSEP  Papers  through  31  August  1980 

During  this  time  period  1  September  1979  -  31  August  1980,  there  are  currently 
an  additional  14  papers  being  processed:  in  press,  in  review,  submitted  and 
being  prepared.  Further,  there  have  been  13  Conference  papers  presented  (or 
to  be  presented). 


2.  RESEARCH  PROJECTS 


RESEARCH  PROJECTS  -  continued 


Project  No.  Title  and  Faculty  Investigator (a) 
It.  SOLID  STATE  ELECTRONICS,  concluded 
B.  Devices,  concluded 


SSD-81-5 


5>  Research  on  Electronic 
Systems  Composed  of  Poly¬ 
mer  Films  and  Planar  De¬ 
vices  for  Transducer 
Applications 


D. 

R. 


W.  Hess 
S.  Muller 


III.  QUANTUM  ELECTRONICS 
QE-81-1  ( 

QE-81-2 

TV.  INFORMATION  SYSTEMS 


^  Thin-Film  and  Guided- 
Wave  Active  Optical 
Devices 

Millimeter  &  Infrared 
Heterodyne  Mixing  & 
Detection 


(terminated) 


erization  of  Image 
Objects  , 


S. 

T. 
J. 

T. 

S. 

T. 


Wang 

K.  Gustafson 
R.  Whinnery 

K.  Gustafson 
E.  Schwarz 
Van  Duzer 


ISS-81-1 

*5  Large-Scale  and  Non- 
Linear  Circuits  Study 

L.  0.  Chua 

E.  S.  Kuh 

1SS-81-2 

r" 

^  Control  of  Large 

Systems 

C.  A.  Desoer 

P,  P.  Varaiya 

ISS-81-3 

^  Computer  Optimization 
of  Electronic  Circuits 

R.  G.  Meyer 

E.  Polak 

A,  Sangiovanni- 
Vincentelli 

ISS-81-4 

(terminating) 

VLSI  Circuits  for 

Future  Computing 

Systems 

A.  M.  Despain 

D,  A.  Patterson 

C.  H.  Sequin 

ISS-81-5 

(terminating 
in  part) 

Analog  Adaptive 

Filters  for 

Speech  Processing 

R.  W*  Brodersen 

D.  G.  Messerschmitt 

ISS-81-6 

An  Analytical  Charact- 

D.  J.  Sakrison 

B.  Consolidated  List  of  Published  JSEP  Papers  through  31  August  1980 
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University  of  California,  Berkeley  Electronics  Research  Laboratory 

Joint  Services  Electronic s  Program  September  30,  1980 

I.  Basic  Research  in  Electromagnetics  Coordinator:  Professor  K.  K.  Mei 

Basic  research  in  electromagnetics  at  EI^L  presently  consists  of  two  parts.  One 
is  the  development  and  application  of  computational  techniques  to  electromag¬ 
netics,  and  the  other  is  the  basic  research  in  theory  and  experimentation  of 
electromagnetics  at  infrared  frequencies.  These  two  research  activities  are 
described  in  the  following. 

Computations  in  Electromagnetics 


Research  in  electromagnetic  theory  was  traditionally  done  within  the  regime  of 
applied  mathematics.  Since  the  advent  of  high  speed  digital  computers,  we 
have  seen  more  and  more  emphasis  on  the  computational  research  being  done  in 
applied  electromagnetics.  In  the  early  1960*8  the  main  thrust  of  electromag¬ 
netic  related  computation  was  in  integral  equation  approaches.1  The  integral 
equation  was  a  natural  choice  in  radiation  problems  because  the  radiation  con¬ 
ditions  at  infinity  could  easily  be  satisfied  by  using  proper  Green's  functions. 
Although  the  integral  equation  method  is  not  too  well  adapted  to  problems  in¬ 
volving  fat  bouies  and  dielectric  bodies,  early  investigators  were  still  quite 
happy  with  it,  because  even  with  its  limitations  it  was  possible  to  solve 
many  problems  which  were  hitherto  unsolvable.  As  more  and  more  problems  were 
solved,  researchers  began  to  look  for  alternatives  to  the  integral  equation 
approach,  which  might  be  better  than  the  Integral  equation  method  in  solving 
fat,  dielectric  or  even  inhomogeneous  bodies.  This  effort  is  represented  by 
the  unimoment  method  which  was  developed  by  K.  K.  Mei  and  his  associates  at 
ERL^,  and  the  several  publications  that  follow** indeed  show  that  the  dif¬ 
ferential  equation  approach  can  be  very  efficient  in  regions  where  the  inte¬ 
gral  equation  method  is  clumsy  and  Impractical.  This  method  is  currently  be¬ 
ing  used  to  investigate  the  scattering  by  buried  dielectric  obstacles,  a  project 
now  supported  outside  the  JSEP. 


The  unimoment  method  is  not  only  good  in  solving  scattering  by  dielectric  bodies, 
it  is  also  an  excellent  technique  for  solving  scattering  by  perfectly  conducting 
bodies.  A  simple  illustration  shows  how  it  works.  Consider  the  scattering  by 
a  conducting  spheroid  as  shown  in  Fig.  1.  Outside  a  circumscribing  sphere,  the 
solution  is  expressed  in  classical  spherical  harmonics  with  undetermined  coef- 
feclents,  and  inside  the  sphere  the  solutions  are  generated  by  solving  the  fi¬ 
nite  element  equations.  The  expansion  coefficients  on  both  sides  of  the  sphere 
arc  determined  by  enforcing  the  continuity  conditions  on  the  sphere.  It  is  In¬ 
teresting  to  see  that  the  finite  element  mesh  area  becomes  larger  if  the  spher¬ 
oid  becomes  thinner,  and  it  becomes  smaller  if  the  spheroid  becomes  fat.  That 
means  the  unimoment  method  is  more  efficient  (more  precise)  for  a  fat  scatter¬ 
ing  than  a  skinny  scatterer.  The  integral  equation  method  works  just  the  oppo¬ 
site  way,  i.e.,  the  surface  area  of  the  scatter  Increases  as  the  scatter  be¬ 
comes  fatter  and  hence  requires  more  computation. 


We  have  already  succeeded  in  using  the  unlmoment  method  to  solve  scattering  by 
dielectric  bodies^,  and  it  is  apparent  that  it  should  work  even  better  when  ap¬ 
plied  to  a  single  fat  metal  scatterer.  So,  in  the  JSEP  we  decide  to  apply  it 
fo  scattering  by  multiple  bodies,  and  initial  success  is  nltrc*ady  observed^. 

A  summary  of  this  paper  follows. 
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Numerical  codes  utilizing  the  unimoment  technique  (  K.  K.  Mei, 

IEEE  Trans.  Ant.  &  Prop.,  AP-2^,  760-766,  Nov.  1974)  have  been 
developed  which  compute  the  scattering  by  metal  and  inhomogeneous 
dielectric  bodies  of  revolution.  The  research  described  In  this 
paper  uses  these  codes  to  solve  the  two-body  scattering  problem. 

First,  interior  field  solutions  (i.e.  interior  to  the  unimoment 
surface)  for  each  scatterer  are  solved  with  the  unimoment  routines 
and  stored  in  the  form  of  system  solution  matrices.  These  matrices 
contain  the  "individual  isolated  scatterer  information"  for  each 
scatterer  and  this  information  need  be  obtained  only  once  for  each 
scatterer.  Second,  an  iterative  scheme  is  used  to  solve  for  the 
scattering  coefficients  for  each  successive  scattering.  Rotation 
and  translation  addition  theorems  are  employed  to  rotate  and  trans¬ 
late  incident  and  scattered  fields  between  coordinate  systems  for 
each  scatterer.  The  iterative  scheme  has  been  tested  on  the  special 
case  of  two  perfectly  conducting  spheres  and  has  proven  to  be  well 
behaved  and  accurate.  (Presented  at  USNC/URSI  Meeting  at  Boulder, 
Colorado,  November,  1978.) 

Other  applications  of  our  computational  skill  are  considerable,  but  we  only 
wish  to  investigate  fundamental  problems.  One  of  the  fascinating  problems 
is  that  of  the  infrared  and  optical  antennas.  At  optical  frequencies  the 
material  behaves  very  differently  from  its  known  properties  at  microwave 
frequencies.  The  compu rational  skill  we  have  developed  for  dielectric 
bodies  can  be  used  to  obtain  physical  insight  into  the  operation  of  infrared 
antennas.  It  will  be  used  to  verify  experimental  observations  and  to  pro¬ 
vide  design  data  to  improve  such  antennas. 


solve  scattering  by  a  spheriod 


The  scattering  by  two  bodies  of  different  shapes  has  been  shown  to  be  success¬ 
ful.  The  results  are  now  available? * However,  the  two  body  configurations 
are  still  limited  to  the  rase  that  the  two  bodies  do  not  penetrate  the  ter¬ 
minating  spheres.  We  are  working  to  remove  this  limitation. 
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Electromagnetics  at  Infrared  and  Optical  Frequencies 

Radio,  microwave,  and  light  are  all  electromagnetic  waves.  The  only  differ¬ 
ence  between  them  is  the  frequency.  However,  both  in  theory  and  application 
they  are  dealt  with  quite  differently.  For  example,  the  guided  transmission 
of  radio,  is  mostly  done  by  coaxial  lines ,  of  microwave  by  metal  waveguides, 
and  of  light  by  beams  or  dielectric  waveguides.  There  are  borderline  fre¬ 
quencies  where  different  transmission  methods  can  be  used.  That  Indeed  is 
the  case  at  aub-mil lime ter  wave  to  low  infrared  waves.  Of  particular  inter¬ 
est  is  the  utilization  of  metal  as  a  guiding  system  at  low  infrared  region. 

It  is  understood  that  metal  is  highly  lossy  at  such  frequencies,  however,  it 
still  has  the  unique  ability  of  concentrating  electromagnetic  power  to  a  very 
small  region,  much  smaller  than  a  wavelength,  which  is  highly  desirable  for 
detection.  In  this  part  of  elect  romagnet  ic  research  we  shall  study  the  scat¬ 
tering  by  impurities  and  transmission  in  optical  fibres,  the  application  of 
metal  as  antennas,  and  field  coupling  from  dielectric  to  metal  waveguide  at 
low  infrared  frequencies.  It  is  worth  noting  that  all  the  technology  de¬ 
veloped  in  this  study  will  eventually  be  applicable  to  millimeter  waves  also. 
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Scientific  Objective 

Multiple  scattering  of  conducting  bodies  is  a  problem  of  both  academic 
and  practical  interest.  The  knowledge  of  multiple  scattering  will  help  many 
design  techniques  In  microwaves;  it  may  alao  provide  insight  into  many  phy¬ 
sical  phenomena  in  field  theories.  However,  scattering  by  electromagnetic 
wave  is  a  difficult  problem  even  without  the  complication  of  multiple  scat¬ 
tering.  Recent  development  of  the  Unimoment  Method  has  advanced  the  art  of 
scattering- related  computation  to  fairly  large  sized  three  dimensional  scat- 
terers^.  One  interesting  extension  of  the  new  method  Is  to  solve  the  problem 
of  multiple  scattering. 

The  objective  of  this  program  is  to  solve  the  scattering  of  two  metal 
bodies  of  revolution.  The  geometries  of  the  two  scatterers  could  be  differ¬ 
ent  and  the  distance  between  them  arbitrary.  The  illuminating  plane  wave 
could  be  arbitrarily  oriented  and  polarized.  The  bistatic  scatterings  of 
the  two  body  scatterings  will  be  our  computational  objective. 

Because  of  the  difficulties  involved  in  theoretical  verification  of  the 
computation,  an  experimental  study  of  the  multiple  scattering  will  be  an  in¬ 
tegral  part  of  the  program. 

State  of  the  Art 

Scattering  by  two  spheres  was  first  investigated  by  Trinka^  in  1935,  and 
then  by  Cermogenova^  in  1963,  for  small  spheres.  Scattering  by  two  arbitrarily 
shaped  bodies  was  treated  by  Zitron  and  Karp^  and  was  limited  to  objects  sepa¬ 
rated  far  apart.  The  latest  attempt  to  solve  the  two  sphere  scattering  prob¬ 
lem  rigorously  was  by  I.iang  and  Lo°.  While  the  multiple  scattering  mechanism 
was  presented  by  Zitron  and  Karp-*  for  arbitrarily  shaped  obstacles,  they  can 
only  obtain  results  for  scalar  fields.  In  our  proposed  program,  the  limita¬ 
tion  on  the  objects  is  that  they  be  axially  symmetric,  in  addition  to  the  un¬ 
avoidable  limitation  of  size  (diameter  less  than  2 X ) . 

Progress  and  Publications  Since  Last  Major  Proposal 

1 .  Experimental  System 

Considerable  experimental  results  have  been  obtained  on  high-frequency 
scattering  by  multiple  spheres.^  These  were  primarily  back  scatter  measure¬ 
ments  because  blstatic  measurements  involved  considerable  difficulties;  con¬ 
sequently  a  refinement  of  the  equipment  was  undertaken. 

A  semi-automated  radar  cross-section  measurement  system  has  been  developed 
where  it  Is  possible  to  obtain  scattering  data  to  be  used  as  verification  of 
the  numerical  results.  The  system  consists  of  an  indoor  image  ground  plane 
with  a  CW  microwave  system  operating  in  the  x-band.  It  utilizes  the  well- 
known  null  balance  technique  for  extracting  the  scattered  field.  A  process 
computer  is  linked  to  the  system  via  a  network  analyzer  which  Is  used  as  a 
coherent  detector.  The  measurement  system  employs  a  reciprocal  measurement 
technique  (RMT)  which  enables  one  to  obtain  bistatic  scattering  data.  The 
process  computer  stores  and  removes  inherent  errors  In  the  system  and  pro¬ 
cesses  the  data  in  a  variety  of  formats. 
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The  measurement  syatem  has  proven  successful  in  verifying  previous  num¬ 
erical  analysis  endeavors.1  The  accuracy  of  the  system  has  been  shewn  to  be 
reliable  for  cross-sections  on  the  order  of  o/X  -  1. 

2 .  Computation 

The  program  of  scattering  by  two  metal  bodies  of  revolution  has  been 
completed.  * " 

Progress  and  Publications  Since  Last  Major  Proposal 
A.  Under  J SEP 

J.  F.  Hunka,  "Electromagnetic  scattering  by  two  bodies  of  revolution/1  Ph.D. 
dissertation.  University  of  California,  Berkeley,  1979. 

K.  K.  Mel,  J.  F.  Hunka  and  S.  K.  Chang,  "Recent  developments  of  the  uni¬ 
moment  method,"  International  Symposium  on  Recent  Developments  In  Classi¬ 
cal  Scattering,  Columbus,  Ohio,  June  1979.  To  be  published  as  a  chapter 
In  Recent  Developments  In  Classical  Scattering,  edited  by  V.  K.  Varadan, 
Academic  Press. 

J .  F.  Hunka  and  K.  K.  Mei,  "Scattering  by  two  bodies  of  revolution/’  IEEE 
Antennas  and  Propagation  Symposium  and  URSI/USNC  1980  meeting,  Quebec  City, 
Canada,  June  1980. 


B.  Under  Other  Sponsorship 

H.  Kao  and  K.  K.  Mei,  "Scattering  by  advanced  composite  bodies  of  revolution," 
USNC/lIRSI  Conference,  Seattle,  Washington,  June  1979. 

S.  Coen  and  K.  K.  Mei,  "inverse  scattering  of  layered  media,"  USNC/URSI  Con¬ 
ference,  Seattle,  Washington,  June  1979. 

S.  K.  Chang  and  K.  K.  Mei,  "Generalized  Sommer f eld  Integrals  and  Field  Expan¬ 
sions  in  Two-Medium  Half-Spaces,"  IEEE  Trans,  on  Antennas  and  Propagation, 
vo 1 •  A  p.  28,  No.  4  pp.  504-512,  July  1980. 

Significant  Accomplishments 

We  have  positively  demonstrated  that  the  Unimoment  Method  we  have  developed 
la  adaptable  to  multiple  scattering  of  arbitrary  shaped  scatterers.  The  com¬ 
putation  is  also  proven  to  be  accurate  and  economical.  The  results  arc  con¬ 
tained  in  a  Ph.D.  thesis  submitted  by  Dr.  John  F.  Hunka,  entitled  "Electro¬ 
magnetic  scattering  by  two  bodies  of  revolution." 

Abstract 


Powerful  numerical  programs  have  been  developed  utilizing  Mel’s  Unimoment  « 

method  which  can  solve  for  the  electromagnetic  scattering  by  metal  and/or  in¬ 
homogeneous  dielectric  bodies  of  revolution.  The  purpose  of  this  research  Is 
to  use  these  existing  programs  to  solve  the  general  three  dimensional  two- 
body-of- revolution  scattering  problem,  thereby  extending  the  unimoment  tech¬ 
nique  to  include  a  class  of  problems  not  axially  syianetrlc.  This  extension 
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of  the  technique  and  the  experience  gained  in  solving  the  two  body  problem 
may  soon  enable  the  analysis  of  the  N-body  of  revolution  problem  and  eventually 
develop  the  capability  to  compute  the  scattering  by  complex  composite  structures 
made  up  of  many  simple  bodies  of  revolution. 

This  dissertation  presents  a  systematic  procedure  for  solving  the  two  body 
of  revolution  scattering  problem  by  interfacing  with  a  unimoment  program. 

Thia  research  presents  a  detailed  investigation  of  spherical  vector  wave 
function  expansions  and  addition  theorems.  A  sea tterer- independent,  modal 
incident  field  formulation  is  developed  as  a  means  of  interfacing  the  mul¬ 
tiple  scattering  and  unimoment  programs.  Also,  an  experimental  scattering 
measurement  system  is  used  for  the  verification  of  the  multiple  scattering 
computations. 

An  iterative  or  "successive  scattering"  approach  is  used  whereby  the 
scattered  field  from  one  body  is  used  as  an  incident  field  for  the  other 
scatterer.  The  process  mandates  the  need  for  suitable  field  expansions 
and  change  of  bases  for  representing  fields  in  the  coordinate  systems  de¬ 
scribing  the  two  bodies.  To  this  end,  spherical  vector  wave  expansions 
and  rotational  and  translational  addition  theorems  are  developed.  For 
each  iteration,  the  scattering  coefficients  are  computed  by  a  system  matrix 
operation  on  the  known  incident  field  coefficients.  The  system  matrix  is 
generated  by  the  unimoment  program  only  once  for  each  scatterer  and  contains 
all  the  geometrical  and  electrical  properties  of  the  scatterer.  A  scattered 
energy  criteria  is  established  to  terminate  the  iteration  process. 

The  verification  of  the  multiple  scattering  program  both  numerically 
and  experimentally  is  discussed  in  depth.  Selected  case  studies  of  two 
body  configurations  are  presented  to  demonstrate  the  utility  of  the  program 
and  to  lay  a  foundation  for  future  study  and  analysis. 
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Scientific  Objective 

Antennae  and  guides  in  the  infrared  and  optical  region  of  the  electro¬ 
magnetic  spectrum  are  expected  to  have  applications  for  radiation  detection 
astronomy,  communications,  and  other  fields.  For  example,  it  has  baen  shown 
that  the  use  of  antennas  can  provide  order-of -magnitude  Improvements  in  infra¬ 
red  detector  performance  [9).  Antenna-coupled  devices  have  also  been  used  to 
obtain  accurate  frequency  measurements  at  frequencies  as  high  as  the  visible, 
and  to  obtain  a  measurement  of  the  speed  of  light  100  times  more  accurate 
than  thoae  preceedlng  (11]. 

The  analysis  and  design  of  infrared  antennaa  Involves  certain  problems 
which  are  different  from  traditional  work  on  lower-frequency  antennas.  Some 
of  these  differences  are  as  follows:  (a)  The  electric  properties  of  metals 
are  altered  at  high  frequencies;  (b)  Typical  antenna  dimensions  are  comparable 
to  wavelength,  hence  for  optical  antennas  novel  antenna  designs  are  required 
to  accomodate  available  microfabrication  techniques;  (c)  correspondingly, 
because  of  the  extreme  smallness  of  optical  wavelengths,  antenna  dimensions 
may  at  times  be  large  compared  with  wavelength.  Excluding  reflector  type 
antennas,  this  is  a  case  not  often  encountered  in  traditional  antenna  researcti. 

Guiding  the  signals  at  optical  frequencies  poses  fundamental  problems  of 
electromagnetics,  whether  optical  waveguides,  strip-lines  or  other  transmission 
systems  are  to  be  used.  In  the  optical  region  the  use  of  dielectric  fiber 
waveguides  in  the  realization  of  transmission  systems  has  grown  markedly. 

There  are  a  great  number  of  theoretical  and  practical  electromagnetic  problems 
that  have  not  as  yet  been  solved  due  to  the  complexity  of  the  hybrid  propaga¬ 
tion  modes  In  the  fiber.  The  particular  basic  problems  of  coupling  between 
guides  and  the  scattering  from  obstacles,  such  as  local  impurities,  are  of 
fundamental  importance  with  respect  to  the  application  of  fibers  to  communica¬ 
tion  systems  and  the  development  of  new  devices. 

State  of  the  Art 


Basic  electromagnetic  problems  have  appeared  ever  since  infrared-  and 
optical-f requency  devices  have  been  in  existence.  The  two  areas  of  concern 
In  this  proposal  Involve  antennas  and  guides  at  these  frequencies. 

Infrared  and  Optical  Antennas  -  It  has  been  recognized  for  some  time  that 
metal- insulator-metal  (MIM)  "optical*’  diodes  (also  known  as  "MOM"  or  "MBM" 
diodes)  can  be  coupled  to  the  electromagnetic  field  by  means  of  an  antenna 
which  is  Identical  with  the  cat  whisker  of  the  diode  [1].  Matarrcae  and  Evenson 
have  shown  that  the  dependence  of  the  output  of  a  whisker  diode  on  its  orienta¬ 
tion  while  illuminated  by  an  Infrared  laser  beam  can  be  explained  by  considering 
the  structure  as  a  simple  long-wire  antenna  (2).  They  have  published  antenna 
patterns  at  33 7p •  Professor  Schwarz  and  his  coworkers  have  obtained  well- 
defined  receiving-antenna  patterns  at  wavelengths  as  short  as  3 p,  also  by  inruns 
of  metal-whisker  antennas  (3].  Professor  Wang  has  applied  the  theory  of  linear 
antennas  to  derive  an  expression  for  the  optical  voltage  variation  with 
frequency  (4).  Planar  evaporated  type  of  antennas  have  been  attempted  but  as 
yet  have  not  met  with  complete  success. 
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Infrared  and  Optical  Guides  -  Until  now  optical  waveguides  have  been 
primarily  of  the  fiber  type.  The  objective  of  thle  proposed  study  Is  to 
Investigate  the  processes  of  mode  conversion  and  radiation  that  come  Into 
play  when  an  object  of  different  dielectric  constant  la  introduced  Into 
either  the  fiber  or  the  fiber  cladding.  Investigations  thus  far  have  been 
restricted  to  the  scattering  from  point  inhoaogenelt les  [5, 6, 7, 8].  Such  a 
treatment  allows  one  to  approximate  the  effects  by  using  a  modified  Rayleigh 
approach  which  examines  the  current  induced  on  the  obstacle.  Only  the  HE^ 
mode  has  been  treated  with  any  detail;  the  results  are  thus  limited  in  use¬ 
fulness.  The  proposed  research  would  deal  with  the  problem  in  a  more  general 
sense.  The  hybrid  modes  will  be  examined  without  reference  to  a  specific  mode. 
Also,  the  dimensions  of  the  scatterer  are  not  restricted  to  the  Rayleigh  region. 

Other  types  of  optical  transmission  lines  such  as  those  using  metallic 
conductors  or  using  portions  of  the  substrates  have  not  been  suff iciently 
investigated  and  seem  to  offer  substantial  new  possibilities.  Here  a  good 
deal  of  basic  research  is  required.  Metallic  conductors  are  comparatively 
lossy,  but  there  are  important  applications,  such  as  coupling  together  elements 
of  antenna  arrays,  where  only  short  distances  are  required.  Meta  1 -conductor 
guides  have  the  virtue  of  closely  confining  energy,  for  easy  coupling  to  diodes 
or  antennas.  Dielectric  waveguides  such  as  ridge  waveguides  are,  on  the  other 
hand,  promising  for  use  in  applications  where  longer  distances  of  propagation 
are  required. 


Progress  and  Publications 

Considerable  progress  was  made  in  design  and  analysis  of  antennas  suit¬ 
able  for  us  In  mil lime ter- wave  Integrated  circuits.  Tapered-dielectric- 
waveguide  antennas  have  been  fabricated  monolithically  in  silicon  and  tested 
at  wavelengths  of  1.3  and  0.119  railimiters.  The  results  have  appeared  as  an 
invited  paper  entitled,  "Antennas  and  Waveguides  for  Far-Infrared  Integrated 
Circuits,"  by  D.  B.  Rutledge,  S.  E.  Schwarz,  T.  L.  Hwang,  D.  J.  Angelakos, 

K.  K.  Mei  and  S.  Yokota,  IEEE  Jour.  Qant.  Elect.,  vol.  QE-16,  pp.  508-516, 
May  1980.  An  abstract  of  this  work  follows: 

Antennas  and  waveguides  for  the  wavelength  range  0.1-3 
mm  are  considered.  Emphasis  Is  placed  on  those  designs 
which  lend  themselves  to  Integration  with  each  other 
and  with  other  components  such  as  diodes.  The  general 
properties  of  FIR  antennas  are  reviewed.  A  novel  sili¬ 
con  waveguide  antenna  is  discussed,  and  its  design ,  sim¬ 
ulation,  fabrication,  and  performance  at  119pm  are  de¬ 
scribed.  This  antenna  has  a  highly  symmetrical,  single- 
lobed  beam  with  3  dB  beamvidths  of  35  and  38°  in  the  E- 
and  H-planes,  respectively.  The  gain  (measured  in  micro- 
wave  simulation)  is  12.8  dB.  This  antenna  is  well- 
suited  for  integration  with  Schottky  diodes. 

The  related  subject  of  FIR  waveguides  is  discussed.  Ex¬ 
periments  with  metal  transmission  lines  at  119pm  are  de¬ 
scribed  and  dielectric  guides  related  to  the  waveguide 
antenna  are  also  considered.  Using  components  such  as 
these  it  may  soon  be  possible  to  construct  receiver  front 
ends  for  this  wavelength  range  in  integrated  circuit  form. 
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Progress  and  Publications,  cont#d 

This  work  is  continuing  under  other  sponsorship.  The  principles  devised 
earlier  are  now  being  applied  to  integrated  antenna-diode  mixer  structures. 

A  different  type  of  antenna  was  devised  for  use  in  multi-mode  detector 
arrays.  Arrays  of  400  microbolometer  detectors  (array  size,  1  cm^)  were  de¬ 
signed  and  tested  at  230  GHz.  Theoretical  analysis  baaed  on  the  EMF  Method 
was  in  good  agreement  with  experiment.  These  results  are  reported  In 
"Planar  Multi-Mode  Detector  Arrays  for  Infrared  and  Mil lime ter- Wave  Applic¬ 
ations,"  by  D.  B.  Rutledge  and  S.  E.  Schwarz,  to  appear  in  IEEE  Jour.  Quant. 
Elect.  Abstract  of  the  paper  is  as  follows: 

A  new  type  of  detector  array  is  described.  By  means  of  a 
suitably  designed  metallic  network,  many  detector  elements 
(each  individually  small  compared  to  wavelength)  are  as¬ 
sembled  into  an  impedance-matched  termination  for  radiation 
incident  normally  on  the  plane  of  the  device.  Residual  re¬ 
actance  is  tuned  out  by  means  of  a  movable  backshort.  An 
array  of  400  bismuth-film  microbolometers  with  a  total  area 
of  one  cm2  has  been  tested  at  215  GHz.  A  coupling  effici¬ 
ency  of  approximately  60%  was  observed.  The  detector  has  a 
D*  of  4  x  10®  cm  Hz*'2/W  at  room  temperature  with  response 
time  on  the  order  of  2  x  10~^  second.  Similar  arrays  of 
Schottky  and  SIS  diodes  can  probably  be  constructed. 

This  work  is  not  at  the  moment  being  continued  but  may  be  revived  as 
required  for  progress  toward  major  program  objectives. 

Both  of  the  abovementioned  publications  acknowledged  partial  JSEP  sup¬ 
port  • 

During  the  reporting  period,  David  B.  Rutledge  completed  his  doctoral 
thesis,  entitled  " Submillimeter  Integrated  Circuit  Antennas  and  Detectors.** 
Dr.  Rutledge  has  now  taken  a  position  as  Assistant  Professor  of  Electrical 
Engineering  at  the  California  Institute  of  Technology.  The  abstract  of 
Rutledge's  dissertation  is: 

Integrated-circuit  techniques  promise  to  make  possible 
reliable  submillimeter  devices  and  arrays.  This  thes¬ 
is  discusses  three  different  integrated-circuit  antennas. 

The  first  is  an  evaporated  V- antenna  in  a  sandwich  of 
crystal-quartz  substrates.  This  was  the  first  evaporat¬ 
ed  submillimeter  antenna  to  have  a  predictable  pattern, 
and  it  has  been  successfully  tested  in  a  plasma  interfer¬ 
ometer  at  119  pm.  The  second  is  a  dielectric-waveguide 
antenna,  made  by  anisotropic  etching  of  silicon.  This 
antenna  has  been  Integrated  with  a  Schottky  diode.  The 
final  antenna  is  a  multi-mode,  tunable  array  of  400  mic- 
robolometers.  This  approach  gives  efficient  coupling 
between  a  source  and  the  detectors.  The  array  achieves 
a  D*  of  2  x  10®  cm  Hz^'2/W“^  at  1.4  mm  with  a  modulation 
frequency  of  1  MHz. 
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Significant  Accomplishment 8 

The  novel  structures  devised  in  this  reporting  period  have  clear  ap¬ 
plications  ,  both  demonstrated  and  potential.  The  single-mode  antennas  are 
promising  candidates  for  inclusion  in  developing  millimeter- wave  IC  techno¬ 
logy,  and  suggest  avenues  for  progress  in  waveguides  for  this  application. 
The  multi -mode  antenna  structure  has  already  been  shown,  when  used  with 
simple  microbolometers,  on  FIR  detector  as  sensitive  as  the  common  pyro¬ 
electric  detectors  and  several  orders  of  magnitude  faster.  This  technique 
can  be  adapted  to  produce  large-area  detector  arrays  based  on  other  kinds 
of  small  detectors,  such  as  diodes  or  SIS  junctions. 

Contacts  with  Other  Agencies 

We  have  been  in  frequent  contact  with  scientists  of  the  U.S.  Army 
Harry  Diamond  Laboratory,  U.S.  Army  Research  Office,  and  the  U.S.  Army 
Night  Vision  Laboratory  at  Fort  Belvoir.  Related  work  in  our  group  is 
presently  sponsored  by  ARO,  and  a  proposal  to  NVL  (involving  practical 
applications  of  novel  techniques  developed  under  JSEP)  is  pending. 

Other  related  work  in  our  group  is  supported  by  the  National  Science 
Foundation.  The  latter  organization  conducts  annual  meetings  on  the 
subject  of  Optical  Communications.  These  have  often  supplied  useful 
motivations  for  research  to  be  conducted  under  JSEP.  We  have  also  had 
useful  contacts  with  scientists  of  the  Aerospace  Corporation. 
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General 


The  main  objectives  of  the  proposed  research  program  are:  (1)  to  characterize 
the  carrier  lifetime  profiles  and  the  deep  levels  In  ion  implanted  semi¬ 
conductors  for  an  understanding  of  the  lifetime  degradation  mechanism,  gath¬ 
ering  information  for  advance  device  design,  and  search  for  improved  tech¬ 
nologies,  and  (2)  to  explore  the  use  and  study  the  properties  of  ternary  and 
quaternary  III-V  compounds  for  advanced  and  improved  Electronic  devices  and 
to  develop  the  necessary  growth  and  processing  technologies  for  these  com¬ 
pounds.  The  goal  of  the  second  objective  is  to  find  an  optimal  approach  to 
and  to  develop  the  MBE  and  LPE  technology  for  integrating  electronic  and  op¬ 
tical  devices  for  fast  communication  and  signal  processing. 

Ion  implantation  is  a  technology  of  great  promise  for  its  ability  to  dope  a 
semiconductor  with  precise  control  of  the  doping  profile.  Among  the  proper¬ 
ties  of  the  ion- implanted  semiconductors,  the  least  known  and  understood  is 
the  carrier  recombination  lifetime.  The  lifetimes  within  the  implanted  lay¬ 
ers  have  not  been  quantitatively  measured  (leaking  current  of  an  implanted 
diode  mostly  reflects  the  lifetime  of  the  substrate).  We  propose  to  measure 
the  lifetime  profiles  with  sub-micron  resolution  using  a  novel  technique, 
which  also  determines  the  built-in  electric  field  profile.  Por  a  basic  under¬ 
standing  of  the  lifetime  degradation,  we  propose  to  correlate  the  lifetime 
profiles  with  the  densities  of  deep  levels  measured  by  transient  C-V  measure¬ 
ments.  Using  these  characterization  techniques  as  tools  we  plan  to  search  for 
the  optimal  annealing  cycle,  to  examine  the  effects  of  laser  annealing,  and  to 
develop  a  novel  technique  of  eliminating  the  defect-causing  stress  in  implanted 
semiconductors  by  co-injp lan ting  two  dopants  of  different  atomic  sizes. 

One  important  area  of  materials  research  is  III-V  compounds  on  account  of  their 
unique  properties  for  a  wide  range  of  device  applications,  such  as  lasers,  photo- 
detectors,  Gunn-effect  devices,  end  field-effect  transistors.  One  fruitful  area 
of  research  is  to  examine  the  unique  possibility  of  combining  various  TTI-V  com¬ 
pound  electronic  and  optical  devices  for  communication  and  signal  processing  and 
to  develop  the  needed  technical  capability  for  fabricating  integrated  opto-ilec- 
tronic  circuits.  In  our  laboratory,  we  have  accumulated  extensive  knowledge 
and  experience  on  the  (GaA£)As  system  from  our  work  on  GaAs- (GaAt)As  lasers. 

So  far,  we  have  used  the  LPE  process  to  fabricate  the  basic  structures  for  opto¬ 
electronic  integration.  However,  the  LPE  technology  has  its  limitations  in  pre¬ 
cise  control  of  film  thickness  and  doping  profiles.  We  propose  that  we  study 
and  develop  the  MBE  technology  as  a  means  of  fabricating  optical  and  elec¬ 
tronic  devices.  Besides  high  electron  mobility  and  direct-gap  band  structure, 
another  advantage  of  III-V  compounds  is  that  the  compounds  are  miscible  so 
that  the  material  parameters  can  be  optimized  in  a  given  application.  The 
(CaIn)(AsP)  quarternary  system  is  especially  interesting  for  being  possible 
to  get  perfect  lattice  match  to  the  InP  substrate.  We  further  propose  that  our 
LPE  work  be  extended  to  the  (Gain) (AsP)  quarternary  system.  I.lquid-phase-epi- 
taxy  growth  conditions  will  be  studied  to  control  the  amount  of  lattice  mismatch 
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and  experiments  will  be  performed  to  examine  the  junction  properties  sen¬ 
sitive  to  interface  defects,  such  as  the  surface  recombination  velocity 
and  interface  state  density  at  the  heteroj unction,  and  to  correlate  the 
measured  properties  with  the  degree  of  lattice  mismatch.  We  believe  that 
such  studies  will  contribute  to  our  fundamental  understanding  of  semiconductor 
heterojunctions.  Lattice  defects  such  as  dislocation  lines  terminating  at 
the  he teroj unction  affect  bulk  as  well  as  Interface  properties.  Lifetime 
and  quantum  efficiency  in  the  (Gain)  (AsP)  LPE  grown  layers  will  be  measured, 
and  thermal  expansion  coefficients  of  the  quaternary  compound  will  be  deter¬ 
mined.  The  proposed  research  program  is  designed  to  study  the  effects  and 
causes  of  defects  introduced  by  lattice  mismatch  so  that  a  basic  picture  of 
semiconductor  he teroj unctions  may  emerge.  An  integral  part  of  the  research 
program  is  the  development  of  suitable  technologies  for  processing  the  quater¬ 
nary  compound,  needed  for  device  fabrication  and  surface  passivation.  Experiments 
such  as  chemical  etching  and  anodic  oxidation  in  which  we  have  done  exten¬ 
sive  work  for  the  GaAs-(GaAi)As  system  will  be  extended  to  the  InP- (Gain) (AsP) 
system. 
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Scientific  Objective 

In  order  to  characterize  the  ion -implanted  semiconductors ,  we  shall  de¬ 
velop  a  novel  technique  of  measuring  the  carrier  lifetime  profile  with  sub- 
micron  resolution  and  the  interface  recombination  velocity,  in  the  presence 
of  arbitrary  and  unknown  electric  field.  The  lifetime  profile  will  be  com¬ 
pared  with  the  density  distributions,  energies,  and  cross  sections  of  the 
deep  levels  that  may  exist.  Using  these  characterization  techniques,  we 
shall  explore  the  possibilities  of  improving  the  lifetimes  by  optimizing 
the  annealing  cycle,  laser  annealing,  or  a  novel  method  of  eliminating  the 
stress  in  heavily  implanted  semiconductors. 

State  of  the  Art 


The  impact  of  ion  implantation  on  the  semiconductor  technology  has  been 
strong  and  broad  and  is  expected  to  deepen  even  more  rigorously.  Basic  un¬ 
derstanding  of  some  of  the  ion  implantation  phenomena  has  not  kept  pace  with 
the  expanding  applications.  In  particular,  the  recombination  lifetimes  within 
the  ion  implanted  layers  have  not  been  studied  with  any  accuracy.  Neither 
have  the  recombination  levels  responsible  for  the  lifetime  degradation  been 
characterized. 


It  is  known  that  large  numbers  of  lattice  defects  are  created  by  ion  im¬ 
plantation  as  the  energy  of  the  ion  is  dissipated  by  displacing  the  host  atoms 
from  their  regular  lattice  sites.  The  density  and  distributions  have  been 
studied* »  .  The  area  density  of  lattice  disorder  measured  by  ion  scattering-* 
increases  with  the  implantation  doses  linearly  until  the  top  surface  layer  be¬ 
comes  amorphous.  Each  implanted  ion  can  create  several  thousand  lattice  de¬ 
fects^  when  the  implantation  is  performed  at  room  temperature . 


Annealing  can  reduce  the  number  of  lattice  defects,  but  typically  some 
of  the  disorder  remains  even  after  high  temperature  (900°C)  annealing*5.  Rem¬ 
nant  lattice  defects  can  have  significant  effects  on  the  electronic  proper¬ 
ties  of  the  material.  Lattice  damage  is  directly  and  indirectly  responsible 
for  the  degradation  of  carrier  lifetime  in  ion  implanted  materials.  Lifetime 
in  silicon  damaged  by  low  doses  of  implanted  ions  may  be  restored  by  anneal¬ 
ing  at  900°C.  At  higher  doses,  some  degradation  remains  even  after  1000°C 
anneal-**  .  Since  the  defects  may  be  electrically  active,  the  carrier  density 
is  influenced  by  the  defect  density.  In  fact,  even  change  of  conductivity 
type  (p  or  n  type)  during  annealing  is  not  uncommon*.  In  some  semiconduc¬ 
tors  the  compensation  is  such  that  the  carrier  concentration  does  not  vary 
with  the  implantation®. 


Most  of  the  lifetime  studies  performed  to  date  measured  the  reverse  leak¬ 
age  current  of  ion  Implanted  pn  diodes®"*^  or  the  generation  current  In  the 
MOS  diodesll.  However,  the  generation  lifetime  cannot  be  equated  with  the 
recombination  lifetime.  Besides,  the  junction  leakage  current  does  not  yield 
reliable  information  about  the  lifetime  in  the  implanted  layer.  The  leakage 
current  is  determined  by  the  lifetime  in  the  space  charge  depletion  layer, 
which  lies  almost  entirely  in  the  unimplanted  substrate.  The  damage  in  the 
implanted  layer  could  actually  improve  the  lifetimes  in  the  depletion  region 
by  gettering  the  lifetime  killers^^.  Generation  lifetime  longer  than  200 

\i9  has  been  measured^. 
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To  our  knowledge,  only  one  group  has  attempted  to  study  the  lifetimes 
within  the  ion  implanted  layers1-**1^.  The  diodes  used  were  phosphorus  dif- 
fued  planar  diodes  on  5fl-cm  p  type  substrates.  The  Junctions  were  about  1 
pm  deep.  To  generate  damage  beneath  the  junction,  1.5  MeV  silicon  ions  and 
2  MeV  carbon  ions  were  Implanted.  The  ion  range  was  estimated  to  be  1.4  pm. 
Even  Chat  study  had  the  same  serious  drawback.  Assuming  a  0.1  ps  lifetime, 
the  diffusion  length  is  about  20  pm.  This  means  that  most  of  the  injected 
carriers  are  stored  in  the  undamaged  material  outside  the  sub-micron  damaged 
region.  The  measured  lifetime  is  probably  not  characteristic  of  the  ion- 
implanted  layer.  Watanabe  et  al.^*  recently  demonstrated  that  the  lifetime 
in  an  ion- implanted  Si  sample  might  vary  over  several  orders  of  magnitude 
over  a  sub-micron  distance.  What  is  needed  is  a  technique  to  quantitatively 
measure  the  lifetime  with  a  resolution  much  better  than  a  diffusion  length. 


Naturally,  the  lifetimes  in  the  implanted  layers  have  not  been  corre¬ 
lated  with  the  deep  energy  levels  that  may  be  present.  Kimerling*”  has  de¬ 
termined  the  deep  energy  levels  introduced  by  proton  implantation  by  a 
junction  capacitance  technique.  These  results  may  not  apply  to  samples  im¬ 
planted  with  other  ions  such  as  boron  and  phosphorus,  for  many  defects  arc 
associated  with  the  implanted  impurit ies^? .  Schulz^*  has  used  ion  implan¬ 
tation  as  a  tool  to  study  the  deep  levels  introduced  by  various  implanted 
impurities  in  silicon.  His  method  of  determining  the  density  and  energy 
of  deep  levels  is  relatively  simple.  However,  it  can  not  determine  the 
spatial  distribution  of  the  deep  levels. 


Progress  and  Publications  Since  Last  Major  Proposal 

All  publications  reported  in  this  section  were  results  of  JSEP  support¬ 
ed  research. 


We  have  proposed  and  analyzed  a  unique  technique  that  can  measure  the 
nonuniform  diffusion  length  L(x)  in  the  presence  of  an  unknown  electric 
field  Ex(x)  with  a  sub-micron  resolution.  This  technique  is  well  suited 
to  the  study  of  ion- implant ed  layers,  which  have  spatially  varying  life¬ 
times  and,  usually,  high  built-in  electric  field.  Previously,  diffusion 
length  or  lifetime  could  not  be  measured  with  better  resolution  than  a 
diffusion  length.  The  principle  of  the  technique  has  been  published 
[C.  Hu,  "Determination  of  Nonuniform  Diffusion  Length  and  Electric  Field 
in  Semiconductors,"  IEEE  Trans,  Elec,  Dev,  ED-25,  822-825,  1978],  Fig.  la 
shows  one  of  the  possible  test  configurations.  The  region  of  interest  is 
sandwiched  by  two  reverse  biased  p-n  or  Schott ky  junctions  at  x  °  0  and 
x  *  w.  The  y  *  0  plane  is  a  cleaved  surface,  perhaps  passivated  with  a 
thin  oxide  layer.  The  electron  beam  scans  from  xe  -  0  to  xs  «  w.  The 
diffusion  length,  and  the  electric  field  can  bo  obtained  from  the  meas¬ 
ured  currents  IQ(xs)  and  Iw(xs). 


-  V‘»">'"oV  -  *,V> 


K  (x  )  -  —  { In [ I  2(I  /I  )' )) 

x  s  q  1  o  w  o 
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Primes  and  double  primes  Indicate  the  first  and  second  derivatives  with  re¬ 
spect  to  x8.  One  may  appreciate  the  novelty  of  the  results  more  when  con¬ 
sidering  that  there  are  no  analytic  solutions  to  the  reverse  problem,  l.e. 
there  are  no  analytic  aolutlons  for  IQ(x8)  and  ^(Xg)  when  L(x)  and  E(x) 
are  general  known  functions. 


•|  electron 
%  beam 


Fig.  1 


A  Ph.D.  student,  Clifford  Drowley  has  carried  out  the  experiment.  He 
has  fabricated  samples  in  the  form  of  Fig.  la.  Results  indicated  that  sur¬ 
face  recombination  at  the  y  -  0  plane  complicates  the  interpretation  of  the 
data.  Another  Ph.D.  student,  Min-hwa  Chi,  is  using  the  configuration  shown 
In  Fig.  Ih,  where  x8  is  varied  by  changing  the  electron  beam  voltage.  In 
that  case,  the  electron  generated  is  not  confined  to  a  single  x8  as  shown 
by  Fverhart  of  this  laboratory  and  the  induced  current  must  be  convoluted 
with  the  carrier  generation  function.  Results  so  far  are  encouraging. 

Meanwhile,  we  have  extended  the  theory  of  the  method  to  include  ac  elec¬ 
tron  beam  excitation.  The  beam  induced  currents  would  contain  both  magnitude 
and  phase  information.  This  expectedly  further  broadens  the  method's  capa¬ 
bility.  t(xs  and  D(t8)  may  be  independently  measured  and  only  one  junction 
is  needed,  not  two.  We  have  succeeded  in  obtaining  a  $5,820  grant  from  the 
University  and  ordered  a  Princeton  Applied  Research  Model  5021  50  MHz  lock-in 
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amplifier  to  implement  this  improved  technique. 


Hr.  Drovley  has  also  been  studying  methods  of  annealing  ion* imp lan ted 
Si.  He  has  successfully  demonstrated  annealing  by  focused  light  from  a 
high  intensity  Xenon  arc  lamp.  This  may  be  an  alternative  to  furnace  an¬ 
nealing  and  laser  annealing.  Small  samples  as  As-implanted  (  100>  Si  (10*^ 
ions/cm^,  200  keV,  resulting  in  an  amorphous  surface)  have  been  recrystal ized . 
Steady  state  temperatures  greater  than  1090°  C  have  been  attained;  rc- 
crystallization  is  quite  rapid  at  these  temperatures  and  occurs  within  10 
seconds  at  700°  C.  TEM  investigations  have  verified  the  recrystall izat Ion 
of  the  amorphous  layers,  and  the  regrown  layers  exlblt  6-10  times  smaller 
sheet  resistivity  than  unannealed  samples.  A  publication  on  this  study  is 
being  prepared. 


19-21 

Mr.  Chi  has  set  up  the  transient  C-V  spectroscopy  experiment  for  profil¬ 

ing  the  trap  densities  and  the  energy,  as  well  as  spatial  distributions.  In 
the  course  of  this  work,  we  have  recognized  that  the  low  speed  of  capacitance 
measurement  and  the  resulting  need  for  sample  cooling  is  a  major  drawback  of 
this  technique.  As  a  result,  we  have  investigated  a  recently  proposed  method, 
transient  Vg  spectroscopy  that  achieves  the  same  end  by  monitoring  the  trans¬ 
ient  change  in  the  threshold  voltage  of  a  MOSFET  after  a  step  change  in  the 
substrate  bias^.  ^  paper  entitled,  "Errors  in  Threshold  Voltage  Measure¬ 
ments  of  MOSFET  for  Dopant  Profile  Determinat ion"  has  been  accepted  for  public¬ 
ation  by  Solid  State  Electronics.  This  new  technique  has  been  demonstrated  in 
our  laboratory. 


Ion-impJ ant ed  semiconductors  are  often  heavily  doped.  We  have  studied  the 
degradation  of  lifetime  in  bulk  heavily  doped  silicon.  This  degradation  is 
generally  attributed  to  Auger  recombination.  This  attribution  has  never  been 
unambiguously  proved.  We  theoretically  examined  an  alternative  mechanism  - 
recombination  through  neutral  donors  or  acceptors.  Somewhat  surprisingly, 
tills  mechanism  could  explain  the  T  «•  l/n^  relationship,  the  weak  temperature 
dependence  of  T  and  the  deduced  neutral  trap  capture  cross-section  is  consist¬ 
ent  with  the  value  measured  at  low  temperatures.  This  and  the  Auger  mechan¬ 
isms  can  be  distinguished  by  studying  heavily  compensated  semiconductors. 

Tills  work  Is  reported  in  C.  Hu  and  W.  G.  Oldham,  Appl.  Phys.  Lett.,  pp.  636- 
639,  Oct.  15,  1979,  an  abstract  of  which  follows: 


Carrier  recombination  in  heavily  doped  semiconductors 
via  shallow  donor  or  acceptor  states  is  analyzed.  The 
results  are  in  general  agreement  with  experimental 
lifetime  observations  including  the  l/n£  or  l/p5  depend¬ 
ence,  and  the  insensitivity  to  temperature  and  to  the 
dopant  used.  Capture  cross  sections  of  about  10*20  cm^ 
needed  to  fit  the  lifetime  data  are  reasonable  for  neut¬ 
ral  traps  and  are  consistent  with  low- temperature  capture 
cross  sections  reported  for  shallow  dopants. 


Significant  Accomplishments 

Wo  have  implemented  and  improved  a  new  deep  level  transient  spectroscopy 
technique  that  uses  MOSFET  structure  and  does  not  require  capacitance  meas¬ 
urement.  \t  Is  a  simpler  technique  and  is  capable  of  measuring  microsecond 
transients  easily.  We  have  demonstrated  the  annealing  of  ion- implanted  Si, 
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Significant  Accomplishments,  cont'd. 

using  focused  light  from  a  CW  arc-lamp.  Compared  to  a  laser,  an  arc-lamp 
is  not  only  much  cheaper  but  also  has  300  times  higher  energy  efficiency, 
the  wide  band  spectrum  also  minimizes  the  interference  effects  of  Si02 
over layers.  We  have  pointed  out  recombination  through  neutral  donors  as  a 
possible  explanation  for  the  degradation  of  lifetime  in  heavily  doped  silicon. 

Interaction  with  Other  Work  Units 

The  characterization  techniques  to  be  developed  under  this  project  will  be 
quite  useful  to  the  nIII-V  Semiconductor  Materials11  project  (SSM-81-2) . 

Results  regarding  the  lifetimes  in  ion- implanted  material  and  improved 
technologies  will  be  helpful  to  all  1C  technology,  devices  and  circuits 
research. 
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Scientific  Objective 

We  propose  to  study  the  properties  of  the  quaternary  compound 
(GaIn)(AsP)  and  the  ternary  compound  Gao. H7In0.53As  concerning  LPE 
growth,  lasing  characteristics,  ohmic-contact  formation,  and  Schottky- 
barrier  height.  The  quaternary  compound  is  important  for  its  spectral 
coverage  best  suited  for  fiber  optical  communications,  and  will  be  stu¬ 
died  as  laser  material  ♦  The  ternary  compound  is  promising  for  its 
theoretically  predicted  high  mobility  and  peak  drift  velocity.  Our 
objectives  are  to  find  the  dominant  factors  affecting  LPE  growth,  lat¬ 
tice  mismatch,  impurity  concentration ,  contact  resistance  and  Schottky- 
barrier  height  in  these  materials,  and  to  correlate  the  material  proper¬ 
ties,  such  as  lattice  mismatch,  doping  profile  and  mater i al -composition 
profile,  with  the  device  performance,  such  ns  lasing  chnrae ter  1st ics , 
turn-on  voltage,  and  1-V  character i sties ,  and  with  the  electrical  meas¬ 
urements,  such  as  carrier  concentration  and  mobility. 

State  of  the  Art 

III-V  compound  semiconductors  possess  many  distinct  and  unique 
advantages  as  an  electronic  and  optical  material.  Most  of  the  compounds 
of  smaller  gap  energy,  such  as  InAs  and  InSb,  have  very  small  electron 
mass  and  thus  exceedingly  high  electron  mobility.  Some  of  the  compounds 
of  larger  gap  energy,  such  as  InP  and  GaAs,  can  be  made  semi-insulating 
by  compensation  either  through  deep  donors  or  through  deep  acceptors. 
The  compounds  are  almost  completely  miscible  except  in  some  systems 
where  a  small  miscibility  gap  exists.  The  mixed  compounds  can  be 
lattice-matched  to  a  given  substrate,  for  example  (Ga -j_xInx)  ( As -j_yPy)  to 
InP,  by  properly  choosing  the  relative  composition  signified  by  x  and  y. 
The  miscibility  and  the  possibility  of  lattice-match  of  these  compounds 
enable  us  to  build  electronic  and  optical  devices  with  materials  dif¬ 
ferent  from  the  substrate  and  thus  to  take  full  advantage  of  the 
material  properties,  for  example  (OaJn)As  on  InP,  with  the  former  for 
its  high  mobility  and  the  latter  for  its  semi-insulating  property. 

Most  of  the  compounds  are  direct-gap  semiconductors  and  hence 
attractive  candidates  as  materials  for  lasers  and  photodetectors.  In 
compounds  where  the  energy  gap  is  direct,  the  energy  separation  between 
the  central  valley  (the  direct  gap)  and  the  satellite  valleys  (above  the 
central  valley)  of  the  conduction  band  varies  from  compound  to  compound. 
This  variation  offers  various  possibilities  for  mixed  compounds.  For 
mixed  compounds  with  relatively  small  separation,  for  example  0.3*  ev  in 
GaAs,  we  expect  the  material  to  exhibit  the  so-called  transferred- 
electron  or  Gunn  effect.  For  mixed  compounds  with  relatively  large 
separation,  we  expect  that  the  material  should  be  able  to  attain  a  much 
higher  drift  velocity  than  the  value  limited  by  the  transferred-electron 
effect.  The  possibility  of  a  high-drift  velocity  makes  these  materials 
attractive  candidates  for  FET.  In  any  case,  on  a  time  scale  shorter 
than  the  inter val ley-scatter ing  time,  the  electron  drift-velocity  is 
expected  to  continue  to  grow  with  increasing  applied  field  to  values 
beyond  those  limited  by  the  transferred-electron  effect.  This 


phenomenon,  called  velocity  overshoot,  makes  compound  semiconductors 
extremely  attractive  for  ultra-high  frequency  applications. 

Among  III-V  compounds,  the  GaAs/(GaAl ) As  system  has  been  most 
extensively  studied  for  both  laser1  and  FET?  applications.  For  the 
laser,  a  great  deal  of  attention  has  been  focused  on  stabilizing  the 
mode  by  incorporating  a  two-dimensional  waveguide  in  the  structure,  and 
remarkable  progress  has  been  made  on  the  Ga As/(Ga A] ) As  laser.  Many 
lasers  exhibiting  linear  light-versus-current  relation,  stable 
transverse  and  lateral  mode  pattern,  and  single  longitudinal  mode  opera¬ 
tion  have  been  developed  at  several  laboratories.  These  include  the 
bur ied-heterostructure  (BH)  laser3»^#  the  channeled-substrate-planar 
(CSP)  laser^t^  and  the  laser  utilizing  thickness  variation  in  the  grown 
layer  to  provide  waveguiding  in  the  lateral  direction?**10.  Among  the 
lasers  referenced  here  are  the  IRW  laser?,  the  CJSF  laser?  and  the  nar¬ 
row  DCC-CSP  laser^  being  developed  and  studied  in  our  laboratory. 
Recently,  the  research  interest  has  shifted  from  the  GaAs/ ( GaAl ) As  laser 
toward  the  quaternary  (Gain) ( AsP)/InP  laser  which  covers  the  spectral 
range  where  optical  fibers  have  exhibited  minimum  loss  (around  1.55 
micrometers)  or  zero  dispersion  (around  1.3  micrometers).  Because  the 
development  of  the  quaternary  laser  is  relatively  recent,  many  problems 
including  mode  stability  and  temperature  sensitivity  have  not  been  fully 
investigated  and  are  yet  to  be  overcome. 

For  the  field-effect  transistor,  the  emphasis  is  on  low  power  and 
high  speed,  and  a  power/speed  product  around  20  f.J  and  a  propagation 
delay  around  100  ps  have  been  reported  by  several  laboratories  (see 
papers  appearing  in  the  special  issue  of  MTT,  Reference  2). To  further 
advance  the  FET  technology,  several  important  problems  need  to  be  over¬ 
come  and  have  received  prominent  attention.  One  of  the  steps  used  in 
the  fabrication  of  the  FET  is  ion  implantation.  It  is  found,  however, 
that  considerable  degradation  of  the  doping  profile  takes  place  during 
annealing  of  Cr  and  S  implants 1 1“*1 3,  Although  anodic  oxidation  has  been 
tried  on  compound  semiconductors1^,  the  interface-state  density  is  con¬ 
siderably  higher  than  that  at  Si/Si0?  interface.  For  lack  of  a  suitable 
insulator,  most  FET's  are  made  of  Schottky-barr ier  gates  in  the  form  of 
MFSFFT.  However,  the  barrier  height  formed  on  Si  and  GaAs  is  generally 
not  as  high  as  we  hope  to  have.  Theoretical  studies1^*^  have  shown 
that  the  degree  to  which  the  barrier  height  at  a  metal-semiconductor 
interface  is  influenced  by  the  metal  work  function  is  strongly  dependent 
on  the  ionicity  of  the  semiconductor.  Because  Si  and  GaAs  are  covalent 
the  barrier  height  in  these  two  materials  is  not  expected  to  change  much 
with  the  metal,  as  confirmed  by  experimental  data1*’.  Therefore,  an 
interface  material,  such  as  slliclde  for  silicon,  must,  be  found  for  com¬ 
pound  semiconductors  between  the  metal  and  the  semiconductor. 

Recently,  the  ternary  compound  InQ^Ga^  i^As  has  generated  consid¬ 
erable  interest,  not  only  as  a  photodetector  for  the  quaternary 
(GaIn)(AsP)  laser,  but  also  as  a  FET1?.  Theoretical  calculations1^  of 
the  low-field  mobility  and  the  peak  drift  velocity  in  the  ternary  com¬ 
pound  show  that  mu  =  9700  cm^/v-s  and  vm  =  2. ft  x  10?  cm/s  are  higher 
than  the  corresponding  values  in  GaAs  and  InP,  even  with  alloy 
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scattering  being  taken  into  account.  The  quaternary  (C.a  In)(AsP)  and  the 
ternary  Ino,  53^0, 247 As  compounds  offer  exciting  research  opportunities 
for  their  outstanding  properties  and  for  their  relatively  unexplored 
status.  Our  research  efforts  will  be  concentrated  at  first  on  those 
aspects  for  which  we  have  acquired  expertise  and  then  extended  to  areas 
which,  we  think,  are  important  and  in  which  we  can  make  significant  con¬ 
tributions  . 

Progress  and  Publications 

We  started  our  work  on  LPE  growth  of  InP  films  on  Inp  substrates 
and  fabrication  of  InP  p-n  diodes  about  one  and  one-half  years  ago  as  a 
M.S.  thesis  project.  During  the  initial  phase,  three  problems  were 
encountered:  (1)  thermal  damage  of  the  InP  substrate,  (2)  large  turn-on 
voltage  (greater  than  2  v)  and  (3)  high  contact  resistance.  The  work  is 
described  in  the  M.S.  thesis  of  K.  D.  Pedrotti,  "Liquid  Phase  Epitaxial 
Growth  of  Indium  Phosphide  for  Device  Applications,"  Department  of 
Electrical  Engineering  and  Computer  Science,  University  of  California, 
Berkeley,  December  1979-  Since  then  these  problems  have  been  overcome, 
and  a  Ph.D.  student  started  the  work  on  LPE  growth  of  the  quaternary 
(Gain) (AsP)  films  on  InP  substrates  about  eight  months  ago.  Using  a 
graphite  cover  plate,  we  have  reduced  thermal  damage  to  an  insignificant 
degree  and  thus  been  able  to  grow  the  first  InP  film  as  well  as  subse¬ 
quent  quaternary  films  of  excellent  quality  when  examined  under  the 
Nomarski  interference  microscope.  Using  evaporated  and  sintered  Au/Zn 
film  as  the  p-side  contact,  turn-on  voltage  around  1.4  v  and  contact 
resistance  about  8  ohms  for  a  15  micrometer  x  300  micrometer  stripe- 
geometry  laser  are  obtained.  The  contact  resistance  can  be  further 
reduced  by  Zn  diffusion.  This  process  will  be  tried  soon.  We  are  happy 
to  report  that  even  without  Zn  diffusion,  we  have  been  able  to  achieve 
lasing  action  in  a  majority  of  the  stripe-geometry  lasers  tested.  The 
lasers  have  a  surprisingly  small  threshold  current  as  low  as  180  mA  and 
show  stable  far-field  radiation  pattern  in  the  fundamental  transverse 
mode.  The  lasers  also  have  a  surprisingly  small  number  of  longitudinal 
modes  which  are  centered  around  1.2*48  micrometers.  We  are  very  pleased 
with  the  progress  we  have  made  during  a  relatively  short  period  of  time. 

Besides  the  quaternary  compound,  we  also  made  significant  progress 
on  GaAs-related  work.  Another  M.S.  student,  Y-S.  Kim,  has  developed  a 
novel,  simple  method  of  measuring  diffusion  length  and  surface  recombi¬ 
nation  velocity  using  a  pn  junction  structure  such  as  a  solar  coll. 
Chopped  light  that  may  be  broadband  is  incident  from  the  bulk  sit**  of 
the  test  cell.  The  photocurrent,  I,  is  measured  as  a  function  of  the 
effective  bulk  thickness  of  the  cell,  W,  which  is  varied  by  varying  the 
junction  reverse  bias.  The  diffusion  length,  L,  and  the  surface  recom¬ 
bination  velocity  of  the  bulk  material  can  be  deduced  from  the  data 
without  curve  fitting.  For  example,  =  Id?(1/I)d2W.  This  work  has 
been  reported  (Y.  S.  Kim,  C.  Hu  and  C.  Drowley,  Record  1*4th  Photovoltaic 
Specialist  Conference,  pp.  596-600,  Jan.  1980). 

In  support  of  the  development  of  III-V  compound  power  devices,  a 
comprehensive  analysis  of  power  MOSFET's  was  made  and  presented  (C.  Hu, 
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"Parametric  Study  of  Power  MOSFET's,"  1979  IEEE  Power  Electronics  Spe¬ 
cialists  Conference,  San  Diego,  June  1979*  pp.  385-395.)  It  was  shown 
that  the  bulk  ohmic  drop  is  the  principal  contributor  to  the  on-state 
voltage  —  a  conclusion  favorable  to  GaAs  high  electron  mobility. 
I urthermore,  using  the  technique  of  variational  calculus,  we  have 
derived  the  optimum  doping  profile  that  minimizes  the  ohmic  resistance 
while  maintaining  a  required  breakdown  voltage.  The  minimized  resis¬ 
tance  is  proportional  to  the  2.5th  power  of  the  breakdown  voltage.  This 
work  appeared  in  IEEE  Trans .  on  Elec .  Dev.,  vol .  ED— 6,  2*43*  March  1 97  9  - 

Interaction  With  Other  Work  Units 

Close  interaction  is  expected  with  the  work  units  on  MBF  and  on 
ion-implanted  semiconductors,  and  frequent  consultation  is  also  antici¬ 
pated  with  other  work  units,  mainly  in  the  IC  and  device  area,  for 
electrical  measurements. 
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Scientific  Objective 

Recent  advances  in  thin-film  growth  technology  have  opened  exciting 
opportunities,  not  only  in  the  development  of  new  device  structures,  hut 
also  in  our  quest  for  a  better  understanding  of  material  properties. 
The  objectives  of  this  research  program  are  (1)  to  master  the  MBE  growth 
technology  for  III-V  compound  semiconductors  and  to  develop  the  associ¬ 
ated  processing  techniques,  such  as  etching  and  masking,  suitable  for 
the  MBE  process,  and  (2)  to  analyze,  to  characteri ze ,  and  to  correlate 
the  properties  of  the  materials  grown  under  different  growth  conditions. 
Our  goals  are  to  seek  an  optimal  approach  and  to  develop  the  needed 
technology  for  fabrication  of  devices  of  micron  and  submicron  dimensions 
where  interface  and  surface  properties  are  expected  to  play  a  crucial 
role. 

State  of  the  Art 

Over  the  last  several  years,  significant  progress  has  been  made  in 
the  thin-film  technology  based  on  epitaxy  for  the  fabrication  of  optical 
and  microwave  semiconductor  devices.  The  three  basic  methods  are 
liquid-phase  epitaxy1  (LPE),  vapor-phase  epitaxy2  (VPE),  including  metal 
organic  chemical  vapor  deposition2  (MO-CVD,  a  form  of  VPE),  and  molecu¬ 
lar  beam  epitaxy3  (MBE).  Among  the  three,  LPE  is  the  most  extensively 
used,  especially  in  the  fabrication  of  double-heterostructure  (DH) 
injection  lasers.  While  the  LPE  technology  has  3ome  distinct  advantages 
over  the  VPE  and  MBE  technology,  it  has  its  limitations.  One  obvious 
disadvantage  is  the  lack  of  control  on  the  film  thickness  and  the  impur¬ 
ity  profile  due  to  the  fast  growth  rate.  Another  serious  drawback  is 
the  impossibility  of  monitoring  film  growth  and  having  the  growth  pro¬ 
cess  temporarily  stopped  for  analytical  examination  of  the  grown  film. 
The  MBE  technology  provides  the  answer  to  these  problems  and  seems  well 
suited  as  a  research  tool  for  the  fabrication  and  study  of  semiconductor 
devices  of  micron  and  submicron  dimensions. 

One  problem  for  which  we  must  find  a  solution  if  we  are  to  improve 
the  quality  of  epitaxially  grown  films  of  compound  semiconductors  is  the 
incorporation  of  unwanted  impurities  during  the  growth  process.  For 
example,  a  dominant  shallow  donor**  was  found  in  LPE-grown  (GaAl)As  and 
this  impurity  was  thought  to  be  sulphur  introduced  into  the  melt  with 
aluminum.  Contaminants^  such  as  transition  metals  were  reported  to  be 
present  in  MBE-grown  GaAs.  The  total  contamination  incorporated  in  a 
MPE-grown  film  may  be  written  as 

C  =  (4/rg)  (d2N/dt  ds)  (1) 

where  C  is  the  impurity  concentration,  (;  is  the  sticking  coefficient,  r-. 
is  the  growth  rate,  and  rj?N/dt  ds  in  the  number  of  contaminant  gas 
molecues  striking  the  growt  h  surace  per  cm2  per  see.  From  kinetic 
theory  of  gases. 


d2N/dt  ds  =  3.5  x  1022P/(MT) 
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where  P  is  the  partial  pressure,  M  is  the  molecular  weight,  and  T  is  the 
background  gas  temperature.  Take  CO  as  an  example  which  has  M  =  28  AMU. 
At  a  typical  growth  rate  r/dg/u  =  1  ^m/hr,  the  product  6  must  be  smaller 
than  7  x  in  order  to  ensure  C  <  10^  cm~3.  This  simple  calcula¬ 
tion  tells  us  that  we  must  be  extremely  careful  with  contaminants  with 
high  sticking  coefficients  and  high  vapor  pressure.  It  is  important, 
therefore,  for  us  to  have  analytical  equipment,  such  as  CMA  and  STMS,  to 
perform  chemical  analysis  of  the  grown  film  and  to  use  electrical  meas¬ 
urements,  such  as  transient  capacitance-voltage  spectroscopy,  to  charac¬ 
terize  the  incorporated  impurities. 

Progress  and  Publications 

Under  JSEP 

This  proposal  was  approved  recently  as  a  supplementary  proposal 
submitted  last  year,  and  laboratory  work  commenced  only  after  the 
arrival  of  the  MBE  system  in  June  of  1980.  In  preparing  for  work  pro¬ 
posed  here,  we  have  done  extensive  literature  surveying  to  get  familiar¬ 
ized  with  the  MBE  growth  process  and  to  know  the  relative  merits,  as 
well  as  limitations,  of  the  MBE,  LPE,  and  VPE  processes.  Before  decid¬ 
ing  on  the  specifics  of  the  MBE  system  we  would  like  to  have,  we  talked 
with  resarchers  in  the  field  and  arranged  several  visits  to  laboratories 
to  sec  demonstrations  of  the  MBE  growth  process.  Finally  we  decided  on 
the  Riber  Model  MOP  MBF  system  and  the  system  was  delivered  and  assem¬ 
bled  in  late  June.  In  the  meantime,  we  have  requested  and  obtained 
financial  support  from  the  Department  of  Electrical  Engineering  and  Com¬ 
puter  Sciences  and  the  Electronics  Research  Laboratory  for  the  use  of 
industrial  funds  towards  the  renovation  and  air-conditioning  of  the  MBE 
room  and  the  purchase  or  a  glove  box  and  a  laminar-flow  hood.  These 
preparations  are  necessary  in  order  to  provide  a  clean  ennvonment  for 
sample  storage,  handling,  and  preparation  in  the  MBE  room.  The  MBE  sys¬ 
tem  was  baked  and  a  vacuum  in  the  middle  to  upper  10“^  range  was 
obtained.  A  separate  vacuum  system  was  built,  to  bake  the  crucible  and 
furnace  assembly  prior  to  their  installation  in  the  MBE  system.  While 
the  initial  work  is  being  carried  out  by  Mr.  Maurice  Bales,  a  highly 
qualified  laboratory  personnel  having  extensive  experience  m  the 
ultra-high  vacuum  system  and  m  the  analytical  instruments,  graduate 
students  are  asked  to  watch  and  help  in  the  work  process  so  they  can 
learn  the  procedure,  and  more  importantly,  get  a  feel  of  the  extreme 
care  needed  m  the  operation.  We  hopp  and  expect  to  start  the  growth 
process  in  late  September. 

Interaction  with  Other  Units 

Close  interaction  is  expected  with  the  work  in  mi crofnbr i ontionn 
under  ;;SD-R2-1,  SS0-82-.\  with  thf  work  on  material  ^harar-top  i  zat  i  on 
under  3SM-82-1  and  3SM-82-2,  and  with  the  work  on  injection  lasers 
underQK-32-1 . 
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Since  1960  solid  state  device  research  and  integrated  circuit  research  have 
been  carried  out  in  a  cotranon  facility  in  the  Electronics  Research  Laboratory. 
We  have  long  maintained  and  observed  the  obvious  benefit;  everyone  enjoys  a 
more  complete,  flexible  experimental  capability.  Further,  as  a  consequence 
of  sharing  experimental  facilities,  there  has  been  a  beneficial  contact  be¬ 
tween  researchers  in  apparently  quite  different  fields.  One  result  is  a  sig¬ 
nificant  transfer  of  technique,  accelerating  progress  in  experimental  tech¬ 
nology  for  all  concerned.  More  subtle  results  are  the  development  of  larger 
research  goals  in  common. 

An  excellent  example  of  the  interaction  of  different  research  groups  is  the 
research  in  electron-beam  lithography  first  begun  by  Professor  Van  Duzer  md 
students  about  5  years  ago.  This  work  was  motivated  by  the  need  for  very 
fine  geometries  in  Josephson  junction  structures.  This  background  stimulated 
work  by  Professor  Neureuther  and  students  in  wider  problems  of  lithography, 
and  more  recent  work  undertaken  by  Professors  Neureuther  and  Oldham  to  ex¬ 
plore  advanced  lithographies  for  integrated  circuits. 

The  IC  fabrication  facility  continues  as  the  focal  point  for  the  collabora¬ 
tion  between  the  device,  design  and  technology  research  activities  in  solid 
state  electronic  devices.  During  the  last  year  the  laboratory  has  been  up¬ 
graded  in  capability  and  cleanliness,  with  the  latest  add  1 1  Ion  being  a  CANON 
FPA  141  projection  printer  funded  through  NSK  and  University  sources. 

Hie  proposals  contained  herein  are  all  advanced  device  basic  studies  which 
use  the  conmon  laboratory  facilities  in  the  experimental  phase.  They  have 
in  common  the  state  of  the  art  techniques  required  for  successful  fabrica¬ 
tion.  Such  processes  as  ion  implantation,  high  resolution  lithography,  plasma 
etching,  and  ion  milling  will  be  used,  as  well  as  more  standard  processes 
such  as  diffusion  and  thin  film  deposition.  The  research  proposals  also  have 
in  common  certain  intermediate  and  long-term  goals.  In  general  the  inter¬ 
mediate  goals  are  the  basic  study  and  development  of  simple  devices  and  de¬ 
vice  models  which  will  point  the  way  toward  the  long  term  goal  of  more  power¬ 
ful,  cheaper,  denser  circuits. 

The  roost  basic  studies  in  this  proposal  are  the  investigations  in  advanced 
lithography.  Progress  in  all  the  device  studies  depends  in  large  measure 
on  advances  in  the  definition  of  the  structures.  The  e-beam  studies  are 
aimed  at  an  analytic  understanding  of  the  factors  controlling  the  ultimate 
performance  of  e-beam  lithography  systems.  The  results  of  such  studies 
would  make  possible,  for  example,  the  application  of  the  extensive  Monte 
Carlo  and  experimental  energy  deposition  characterization  t  o  a  wide  class 
ot  new  materials.  Including,  e.g.,  losephson  superconductor  HtruriurrH.  Tlx* 
proposed  x-ray  and  photolithographic  studies  similarly  promise  to  aid  in  the 
advance  of  the  integrated  circuit  and  transducer  state  of  the  art. 
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The  proposal  on  Josephson  device  studies  Is  aimed  at  exploring  the  fundamental 
limits  of  the  performance  of  these  devices  in  large  scale  electronic  applica¬ 
tions.  Van  Duzer  and  co-workers  have  demonstrated  that  they  can  perform  sub- 
micron  lithography  in  a  modified  scanning  electron  microscope.  They  have 
also  constructed  a  number  of  novel  working  Josephson  structures.  They  now 
propose  to  combine  the  technologies  and  make  really  dense  circuit  cells  to 
explore  the  problems  in  this  new  domain  of  device  size.  A  number  of  funda¬ 
mental  theoretical  problems  will  be  encountered  and  resolved  (e.g,f  storage 
of  a  single  flux  quanta  appears  impossible  in  the  tiniest  loops)  in  the 
course  of  these  studies;  it  promises^  to  be  a  very  exciting  time  in  the 
Josephson  effect  research. 

Another  proposal  concerns  the  exploration  of  novel  processes  for  high  per¬ 
formance  LSI.  This  is  a  continuation  of  a  program  which  has  already  demon¬ 
strated  a  unique  CMOS  process  for  fabrication  of  analog  and  digital  circuits 
on  the  same  chip.  In  the  present  proposal  emphasis  is  placed  on  the  search 
for  simpler,  denser  processes  which  make  use  of  such  tricks  as  self-alignment 
to  make  possible  a  generation  of  higher  performance,  lower  power  integrated 
circuits.  One  largely  unexplored  circuit  approach  uses  silicon  MESFET's  as 
the  active  devices.  In  this  research,  various  methods  of  fabricating  MF.SFET 
integrated  circuits  will  be  explored.  The  potential  reward  is  very  large, 
as  MESFET's  are  capable  of  functioning  at  subnanosecond  speeds  but  with  a 
power-delay  product  orders  of  magnitude  better  than  present  approaches. 

The  proposal  for  noise  studies  in  integrated  circuits  is  aimed  at  finding  new 
ways  of  achieving  low-noise  amplification  in  integrated  circuits.  The  use  of 
ion  implantation  to  achieve  special  low-noise  bipolar  devices,  and  to  fabri¬ 
cate  monolithic  JFET  structures  is  one  avenue  which  will  be  pursued.  Another 
area  with  potentially  very  large  payoff  is  low-noise  MOS  devices  for  all  MOS 
analog  circuits.  It  appears  that  a  significant  portion  of  the  rapidly  ex¬ 
panding  telecommunications  IC  market  can  be  satisfied  with  MOS  LSI  provided 
low-noise  analog  amplification  can  be  realized.  One  approach  to  be  pursued 
is  the  use  of  buried  channel  MOS.  By  performing  a  second  implant  In  the 
channel,  the  electrons  can  be  confined  to  the  bulk,  which  will  have  the  ef¬ 
fect  of  eliminating  the  noise  from  surface  states. 

The  final  proposal  in  the  device  area  deals  with  the  need  for  all -elect roni c 
transducers.  As  electronics  has  become  cheaper,  an  increasing  bottleneck  in 
its  application  to  control  problems  is  the  lack  of  inexpensive  transducers. 
Such  transducers  often  can  be  based  on  the  piezoelectric  effect.  The  pur¬ 
pose  of  the  project  is  to  develop  a  novel  method  for  the  production  of  piezo¬ 
electric  polymer  thin  films  directly  from  monomer  sources,  thereby  elimin¬ 
ating  many  of  the  constraints  inherent  in  incorporating  such  films  into 
solid  state  devices. 


Significant  Progres.s 

Progress  in  the  program  supported  by  JSEP  includes; 

(1)  Measurement  of  tin*  absorption  spectra  of  several  electron-beam  resists 
an  d  pho  t  o  r e s i s l s ; 
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(2)  Demonstration  of  a  capacitance  measurement  to  determine  the  remaining 
resist  thickness  vs.  time  during  development; 

(3)  A  new  technique  for  noise  reduction  in  switched  capacitor  filters  has 
been  developed. 

(4)  Observation  that  depletion-mode  devices  exhibit  a  region  of  buried- 
channel  conduction  in  which  the  input-referred  flicker  noise  is 
greatly  reduced. 

(5)  An  experimental'*'  low-noise  chip  has  been  fabricated  containing  high- 
frequency  bipolar  devices,  3-terminal  and  4-terminal  high-frequency 
JFETS  as  well  as  compatible  I^lT  logic. 
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Scientific  Objective 

The  complexity  of  large-scale  integrated  circuits  continues  to  grow 
rapidly.  As  the  dimensions  of  individual  devices  are  reduced  to  the  sub¬ 
micron  range,  mask-making  and  mask  alignment  become  more  difficult  and 
costly.  We  wish  to  Investigate  simple  new  fabrication  processes  which 
might  provide  higher  performance  devices  and  circuits  than  those  obtain¬ 
able  with  existing  processes,  while  employing  fewer  masking  steps.  Other 
desirable  features  in  potential  new  processes  are  a  maximum  degree  of  inher¬ 
ent  self-alignment  of  device  features  and  isotropic  etching  and  impurity  do¬ 
ping  processes  to  give  improved  control  of  device  geometries. 

State  of  the  Art 


Standard  fabrication  processes  for  large-scale  Integrated  (LSI)  circuits 
are  generally  specialized  for  digital,  analog,  or  charge- coup led  device  (CCD) 
functions. ^  The  processes  which  provide  highest  circuit  density  (e.g. , 
I^L  and  silicon-gate  NMOS)  typically  have  performance  limits  such  that  digi¬ 
tal  logic  speeds  faster  than  5  ns  are  very  hard  to  achieve,  and  analog  cir¬ 
cuit  functions  are  realized  only  with  serious  limitations  on  accuracy  and 
signal-to-noise  ratio.  The  low  power  consumption  of  complementary  MOS  cir¬ 
cuits  and  the  unique  signal-processing  capability  of  CCD's  are  obtained  only 
at  the  expense  of  complex,  specialized  fabrication  processes.  To  date,  new 
fabrication  processes  such  as  Ion  implantation,  local  oxidation,  self-aligned 
masking,  plasma  and  ion  etching,  and  electron-beam  pattern  definition  have 
been  used  in  the  evolutionary  improvement  of  existing  integrated  circuit  pro¬ 
cesses,  without  enhancing  their  mutual  compatibility^ we  propose  research 
on  new  LSI  processes  which  exploit  the  new  fabrication  techniques  to  achieve 
better  circuit  versatility  and  performance. 


Progress  and  Publications 

A  Si  MESFET  model  for  circuit  analysis  has  been  derived  [8].  It  in¬ 
cludes  the  following  second  order  effects: 

o  Distributed  Channel  Charge 
•  o  Drain  Feedback 

o  Substrate  Bias  Effects 
o  Source-Drain  Punch-Through 
o  Variable  Capacitance  Effects 
©  Temperature  Effects 

The  model  is  mostly  empirical  and  has  been  implemented  in  the  Integrated 
circuit  simulation  program  SPICE  2. 

In  recent  years,  there  has  been  an  increased  interest  in  GaAs  MBS PETS 
for  high-speed  low-power  integrated  circuits.  However,  an  accurate  device 
model  for  circuit  analysis  is  still  lacking  as  drift  velocity  saturation 
dominates  the  characteristics  of  these  devices.  Analytical  models  proposed 
previously  use  a  coarse  approximation  of  the  mobility  as  a  function  of  the 
electric  field  and  yield  the  drain-source  current  in  implicit  form.  Such 
solutions  are  unacceptable  for  CAD  models  as  they  require  iteration  within 
the  device  model. 
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Progress  and  Publications,  contfd 


We  have  found  an  explicit  solution  for  the  drain-source  current  and 
saturation  voltage  by  approximating  the  mobility  y,  by 


U(E) 


E 


where  y  is  the  low-field  mobility,  VCAT  is  the  saturated  drift  velocity, 

and  E  is  the  electric  field.  This  approximation  for  the  mobility  is  widely 
used  in  two-dimensional  device  simulation  programs  like  CADET.  A  new  gen¬ 
eral  model  for  JFETs  and  MEFETs  that  contains  this  approximation  of  the 
mobility  and  that  also  includes  the  previously  mentioned  second  order  ef¬ 
fects,  sub threshold  conduction,  and  short-length  and  narrow-width  effects, 
has  been  derived.  The  model  is  mostly  physical  and  a  publication  of  it 
is  in  preparation.  As  the  model  includes  velocity  saturation,  it  is  also 
valid  for  III-V  compound  MESFETs  like  GaAs. 

The  MESFET  model  in  SPICE  2  has  been  used  to  design  the  following 
circuits : 


1  x  IK  BIT  SHIFT  REGISTER 


Speed 
Power 
Area 

Organization 
#  of  Elements 
On-chip  Clock  Drivers 
Fully  ECL  compatible 


0-500  M  bits/S 
1.2  W 

3x4  mm ^ 

SPS 

30,000 


6  x  6  BIT  AND  10  x  10  BIT  PARALLEL  MULTIPLIERS 


Speed 

Power 

Area 

Speed 

Power 

Area 

Organization 
//  of  Elements 


31 

0.5 

1.2 

39 

0.7 

1.5 

CSA, 


ram 


ns 
W 

x  1.  3 
ns 
W 

x  1.6  mm^ 
no  carry 


} 


6  x  6  BIT 


10x10  BIT 


look-ahead 


2000-3000 


TTL  Data  Compatible 


-  SEVERAL  SMALL  CIRCUITS 


3  Divide-by-two* s 
6  Ring  Oscillators 
2  Drivers  (ECL  and  TTL) 

The  shift  register  serves  as  a  test  vehicle  for  the  yield  of  this  new 
technology.  Together  with  the  multipliers,  the  circuits  also  demonstrate 
the  low  speed-power  product  of  MESFETs  in  practical  applications  and  the 
possibility  of  interfacing  to  other  technologies  like  ECL  and  TTL.  The 
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Progress  and  Publications,  cont'd 

smaller  circuits  serve  to  support  a  computer  study  of  Si  MESFET  circuit 
techniques  with  experimental  data.  The  work  on  these  circuits  is  being 
done  in  cooperation  with  Tektronix,  Inc.,  Oregon.  Presently  the  circuits 
are  being  fabricated  and  the  first  parts  are  expected  in  the  Fall  of  1980. 

It  should  be  noted  that  the  above  circuits  use  the  ',all-depletion,, 
circuit  technology  and  that  application  of  the  “enhancement-depletion11 
technology  will  significantly  reduce  the  power  and  area  of  the  above  cir¬ 
cuits.  The  latter  circuit  technology  requires,  however,  a  tight  control 
of  process  parameters. 

In  order  to  investigate  process  parameter  variations,  a  device  char¬ 
acterization  program  has  been  written.  It  runs  on  an  automatic  device 
characterization  system  that  uses  a  Tektronix  4052  as  controller  and  en¬ 
ables  the  user  to  (quickly)  obtain  statistical  parameters  for  a  small 
dimension  MESFET  model  for  circuit  analysis.  All  data  is  stored  on  disk 
and  is  used  as  a  feedback  for  the  process  and  for  further  development  of 
the  new  JFET/MESFET  model  for  CAD. 

Significant  Accomplishments 

A  Si  MESFET  model  has  been  derived  [1]  and  has  been  implemented  in 
the  integrated  circuits  simulation  program  SPICE  2.  This  model  contains 
distributed  charge  effects,  drain  feedback,  substrate  bias  effect,  source- 
drain  punch-through,  variable  capacitance  effects,  and  temperature  effects. 

A  rewer  model  has  been  derived  which  includes  velocity  saturation, 
subthreshold  conduction,  and  small  dimension  effects  in  addition  to  the 
above.  The  model  is  mostly  physical  and  applies  to  GaAs  JFETs  and  MESFETs 
as  well.  A  publication  of  the  model  is  in  preparation. 

The  SPICE  2  MESFET  model  has  been  used  to  design  a  1  x  1  k  bit  stat¬ 
ic  shift  register  that  is  an  order  of  magnitude  faster  than  commercially 
available  registers  of  that  size,  and  6x6  bit  and  10x10  bit  parallel 
multipliers  that  have  an  order  of  magnitude  lower  power  than  commercially 
available  multipliers  with  similar  speeds.  The  work  on  these  circuits  is 
being  done  jointly  with  Tektronix,  Inc.,  Oregon,  and  experimental  results 
are  expected  in  the  Fall  of  1980. 


Interactions  with  Other  Work  Units 


During  the  course  of  this  research,  interactions  with  the  groups  work¬ 
ing  with  Professors  W.  G.  Oldham,  A.  R.  Neureuther,  P.  R.  Gray,  R.  G.  Meyer 
and  R.  W.  Brodersen  were  maintained. 
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Scientific  Objective 

We  plan  to  explore  the  fundamental  limitations  in  lithography  as  micro- 
circuit  feature  sizes  decrease.  The  emphasis  in  this  work  is  on  the  exposure 
and  development  of  ultra-high  resolution  resists  for  optical  x-ray  and  e-beam 
lithography.  Exposure  models  and  development  models  will  be  formulated  based 
on  experimental  characterization  of  resist  properties.  A  major  emphasis  will 
be  placed  on  the  prediction  and  characterization  of  final  resist  cross-sec¬ 
tions.  The  limits  of  minimum  useful  exposure  are  of  key  Interest. 

State  of  the  Art 


Lithography 

The  first  major  advance  in  understanding  of  photolithography  came  in  the 
series  of  papers  by  Dill  and  co-workers,  of  which  references  1  and  1  are  repte- 
sentative.  A  simple  linear  model  for  the  single  wavelength  exposure  and  de¬ 
velopment  of  positive  photoresists  of  the  diazo  type  enabled  quantitative  pre¬ 
dictions  of  the  effects  of  exposure  and  development  parameters  on  resist  dis¬ 
solution  rate.  Very  recently  the  model  was  applied  to  a  multiwavelength  ex¬ 
posure^.  However,  the  basic  model  was  not  extended  to  account  for  nonlinear 
Interaction  between  the  wavelengths,  effects  which  we  have  observed  in  our 
characterization  at  Berkeley.  Because  all  practical  exposure  systems  are 
necessarily  polychromatic,  it  is  essential  to  modify  the  basic  resist  model 
to  correctly  account  for  multiwavelength  exposure  interactions. 

Negative  photoresists,  although  the  dominant  factor  In  1C  fabrication,  have 
not  been  modeled  in  any  comparable  fashion  to  the  positive  resists.  The  de¬ 
velopment  process  is  very  complex  and  It  is  only  recently  that  quantitative 
studies  of  development  kinematics  have  been  reported.  The  most  recent  and 
ambitious  work  is  that  of  Tu  and  Ouano^  which  models  polymer  dissolution  In 
general.  The  particular  experimental  case  examined  Is  etching  of  polystyrene 
in  MEK,  but  the  results  appear  applicable  to  classical  negative  photoresists. 

A  much  sinpler  negative  resist  model  which  ignores  the  kinematics  entirely 
was  proposed  by  Ko  and  Oldham  and  is  in  the  process  of  being  tested  at  this 
t  ime* . 


X-ray  and  e-beam  lithography  have  Intrinsically  such  a  very  fine  resolu¬ 
tion  capability  that  an  accurate  resist  model  would  appear  as  essential.  How¬ 
ever,  the  weak  source  intensity  problem  and,  in  the  case  of  e-beams,  the  prox  ¬ 
imity  effect  both  fo»*ce  a  very  detailed  look  at  exposure  and  development  mech¬ 
anisms.  In  fact  in  both  technologies  a  direct  tradeoff  exists  between  expo¬ 
sure  time  and  resolution0*  .  A  figure  of  merit  such  as  the  exposure  timc-line- 
width  product  depends  on  the  resist  material  and  the  exposure  and  development 
parameters. 

A  simple  exposure-development  model  for  the  classic  PMMA  resist  was  pro¬ 
posed  by  Creeneich®’  .  This  model  has  been  rather  successful  In  predicting 
x-ray  and  e-beam  produced  resist  prof ilee^*  ,  However,  the  development 
model,  based  purely  on  molecular  weight,  has  recently  been  questioned^.  It 
is  suggested,  e.g.,  that  resist  porosity  as  produced  by  gas  evolution  dining 
exposure  plays  a  role  in  development  kinematics  .  In  this  case  post  exposure 
processing  and  developer  formulation  might  be  key  parameters  affecting  fln.il 
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profiles.  Resist  characterization  and  derivation  of  special  cases  of  the 
kinematic  model^  are  required  to  attack  this  Important  problem. 

There  are  very  basic  problems  to  be  solved  in  the  use  of  resists  at  the 
limits  of  low  exposure.  Statistical  effects  ultimately  limit  the  usefulness 
of  any  resist^.  in  the  case  of  x-rays  an  effective  photo  electron  range  has 
been  determined  for  various  x-ray  wavelengths^^ * ^ #  a  complete  distribution 
function  should  be  determined  relating  shape  and  relative  numbers  of  photo 
electron  energy  deposition  profiles  to  x-ray  wavelength.  A  statistical  model 
of  exposure  could  then  he  developed  which  yields  the  local  three  dimensional 
energy  deposition  to  exposure. 

In  addition  to  the  stat ist ical  e l  f  ects  ,  long  range  electron  scattering 
plays  a  significant  limiting  role  in  high  resolution  e-beam  lithography. 

Recent  work  on  resist  exposure  by  electrons  can  be  divided  into  two  approaches, 
;inalytle  modelling  based  on  scattering  and  energy  loss  formulae”,  and  Monte  Carlo 
calculations  carried  out  in  computers.  The  latter  have  been  refined  over  the 
past  several  years,  and  yield  results  in  good  agreement  with  experiment 
They  have  the  disadvantage  that  each  new  physical  embodiment  requires  a  new 
calculation,  and  each  calculation  is  expensive  and  time  consuming.  The  re¬ 
sults  cannot  be  generalized  as  easily  as  analytic  formulations,  but  since  the 
Monte  Carlo  procedure  simulates  what  actually  occurs  in  the  solid,  it  is  po¬ 
tentially  as  accurate  as  the  cross-sec t ions  for  scattering  that  are  used,  pro¬ 
vided  enough  trajectories  are  <  .ilculated.  (We  are  developing  Monte  Carlo  pro¬ 
grams  on  a  different  grant  fur  comparison  purposes.) 

Analytic  modelling  of  the  electron  Interaction  with  the  solid  can  give  re¬ 
latively  simple  formulae  which  arc  surprisingly  accurate  given  the  approxima- 
t Ions  on  which  the  modelling  is  based.  These  easy-to-inte rpret  formulae  pro¬ 
vide  a  gooJ  physical  "feel"  tor  the  problem,  and  can  provide  accurate  numer¬ 
ical  results  In  some  cases.  Most  Incorporate  the  single  backscatter  model  of 
Kverhait^,  the  forward-se.it  tering  model  used  by  Nosker^®,  and  the  comb  in  a- 
t  ion  of  these  explored  by  uieeneich  and  Van  Dyzer^  and  by  Hawryluk,  Hawryluk 
and  Smtth^.  Recently,  McAifee*^  and  others^"  have  extended  Everhart’s  hack- 
scattering  model.  This  extension  suggests  an  important  avenue  for  future 
studies  -  the  formulation  of  a  simple  analytic  model  for  energy  deposition 
profiles. 


)'r 


ogress  Since  last  Major  Proposal 


Hu-  emphasis  in  the  second  year  has  been  primarily  on  electron  beam 
lirhogr.iphv  with  additional  work  on  developing  automatic  tools  for  charact- 
'■rf/lng  positive  type  photoresist.  Studies  of  electron  energy  dissipation 
in  layered  media  including  a  new  physical  model  based  on  the  random  nature 
inelastic  scattering  events  were  completed.  A  comparison  of  direct  and 
■  onr  minus  v  slowing  down  apprnximat ion  Monte  Carlo  methods  for  simulating 
>•  eitron  scattering  was  made.  Flectron  beam  alignment  signals  were  studied 
experimentally  as  well  as  with  suitably  modified  versions  of  the  Monte  Carlo 
energy  deposition  simulation  programs.  An  experimental  Investigation  of  the 
use  Of  posit  ivie  photoresist  A7.1350J  as  a  non-swelling,  high-resolution  neg- 
attve  electron  beam  resist  was  carried  out  in  collaboration  with  Siemens 
I  orporat ion .  Several  Improvements  in  the  accuracy  of  an  automatic  system 
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for  measuring  the  exposure  parameters  of  positive  photoresist  were  made 
and  a  method  of  determining  the  development  parameters  through  the  in 
situ  measurement  of  absorption  was  explored.  The  results  of  these  ef¬ 
forts  will  now  be  described  in  somewhat  greater  detail. 

A  new  model  for  electron  scattering  was  introduced  into  Monte  Carlo  sim¬ 
ulation  algorithms.  The  new  model  allows  random  inelastic  scattering  e- 
vents  and  does  not  assume  that  the  electron  is  continually  losing  energy. 
Thus  it  is  statistically  possible  for  some  electrons  to  pass  through  thin 
films  without  loss  of  energy  as  is  observed  experimentally.  Electron  dis- 
ipation  in  layered  media  was  also  considered  particularly  as  it  applies  to 
obtaining  very  high  resolution  on  thin  substrates.  Both  theory  and  exper¬ 
iment  showed  that  higher  resolution  in  electron  beam  lithography  can  be 
achieved  by  using  thinner  electron  sensitive  polymer  films,  higher  beam 
voltages,  and  thinner  substrates.  More  details  are  given  in  the  PhD 
thesis  of  I.  Adesida,  entitled  ’‘Electron  Energy  Dissipation  in  Layered 
Media,"  and  described  by  the  following: 

Abstract 

The  influence  of  electron  scattering  on  the  resolution  of 
electron  beam  lithography  has  been  studied  theoretically 
and  experimentally.  Theoretical  calculations  of  the  spat¬ 
ial  distribution  of  electron  energy  dissipation  in  a  very 
thin  electron  sensitive  polymer  film  coated  on  various 
thicknesses  of  silicon  substrates  have  been  made  using  two 
different  Monte  Carlo  models.  Experimentally,  we  have  ex¬ 
posed  elementary  lithographic  patterns  on  the  layered 
structures  mentioned  above.  Both  theory  and  experiment 
showed  that  higher  resolution  in  electron  beam  lithography 
can  be  achieved  by  using  thinner  electron  sensitive  poly¬ 
mer  films,  higher  voltages,  and  thinner  substrates. 

Fundamental  characteristics  of  electron  penetration  through 
thin  films  of  aluminum  and  polymethylmethacrylate  (PMMA) , 
such  as  energy  distribution  of  transmitted  electrons,  were 
theoretically  calculated  using  two  Monte  Carlo  approaches. 

One  approach  is  conventional  in  that  the  energy  loss  of  the 
electron  along  its  trajectory  is  computed  using  the  contin- 
uous-slowing-down  approximation  (CSDA) .  The  essential  feat¬ 
ures  of  the  second  approach  are  the  inclusion  of  the  random 
nature  of  inelastic  scattering  events,  and  also  the  exten¬ 
sion  of  Gryzinski’s  semi-empirical  expression  for  core  elect¬ 
ron  excitation  to  valence  electron  excitation  through  the  use 
of  an  appropriate  mean  binding  energy.  A  detailed  comparison 
with  available  experimental  results  showed  that  the  second 
Monte  Carlo  approach  gave  a  more  realistic  answer  for  the 
transmitted  electron  energy  distribution. 

In  addition,  both  Monte  Carlo  approaches  were  used  to  calcul¬ 
ate  the  electron  energy  dissipation  profiles  in  a  very  thin 
electron  sensitive  polymer  film  1000  X-thick  PMMA)  on  sil¬ 
icon  substrates  for  energetic  incident  electrons.  Energy  dis¬ 
sipation  profile  calculations  were  made  for  dot,  lino,  and 


Progress  Since  Last  Major  Proposal,  cont 1 d 


parallel  line  exposures.  The  effect  of  varying  the  sub¬ 
strate  thickness  and  exposure  dosage  were  studied  in  de¬ 
tail.  The  Gaussian  scattering  parameters  n„,  6^,  and 
were  also  evaluated  theoretically  with  the  CSDA  model. 

For  the  experiments,  different  thicknesses  of  silicon  sub¬ 
strates  were  fabricated  using  an  anisotropic  etching  method. 
Layered  samples  of  thin  PMMA  film  on  thin  films  of  silicon 
were  exposed  with  a  stationary  beam  (dot  exposure)  and  a 
scanned  beam  (line  exposure).  After  exposure,  the  samples 
were  developed,  and  the  resulting  relief  patterns  measured 
in  a  scanning  electron  microscope.  Experimental  results 
were  compared  in  detail  with  results  of  both  Monte  Carlo 
approaches.  Agreement  between  experiment  and  both  Monte 
Carlo  approaches  was  good,  which  indicates  that  both  ap¬ 
proaches  can  be  used  to  predict  energy  dissipation  profiles 
for  electron  beam  lithography.  Results  also  showed  that 
higher  resolution  in  electron  beam  lithography  and  a  re¬ 
duction  in  proximity  effect  can  be  attained  with  a  reduc¬ 
tion  of  electron  scattering  in  the  exposed  sample.  The 
reduction  in  electron  scattering  is  achieved  by  using 
thinner  films  of  electron  beam  resists  and  thinner  sub¬ 
strates  coupled  with  higher  incident  electron  beam  voltages. 

To  study  electron  beam  alignment  signals  theoretically,  the  CSDA  Monte 
Carlo  energy  deposition  algorithms  were  modified  to  produce  information 
about  electron  trajectories  in  the  presence  of  a  tapered  step  in  the  sub¬ 
strate  material.  This  nonplanar  substrate  extension  allows  for  the  possi¬ 
bility  that  electrons  which  leave  the  substrate  at  one  point  may  re-enter 
the  substrate  at  another  position.  The  algorithms  were  also  modified  to 
collect  position  and  momentum  statistics  within  the  substrate  as  well  as 
the  angular  and  energy  statistics  of  backscattered  electrons.  The  stat¬ 
istical  information  within  the  substrate  is  useful  in  developing  approxi¬ 
mate  analytical  methods  of  characterizing  alignment  signals.  Good  agree¬ 
ment  has  been  obtained  between  experimental  and  simulation  results  for 
54.6  degree  anisot Topically  etched  registration  marks  in  silicon.  Some 
progress  has  been  made  toward  generating  a  normalized  or  universal  curve 
for  characterizing  alignment  signals  as  a  function  of  the  Bethe  range  of 
the  electrons.  The  abstracts  of  two  publications  on  this  subject  follow. 

The  first  paper,  ’’Monte  Carlo  Simulation  for  Electron  Beam  Microfabrication 
by  Lin,  Adesida  and  Neureuther,  Appl.  Phy.  Lett.,  vol.  36(8),  15  April  1980 
is  abstracted  as  follows: 

Abstract 

Monte  Carlo  calculations  based  on  a  continuous-slowing- 
d own- approximation  (CSDA)  model  have  been  carried  out  to 
study  the  backscattered  electron  image  formation  of  a 
34.6°  tapered  etched  step  registration  mark  in  silicon. 

It  is  found  that  the  energy  signal  gives  slightly  better 
contrast  and  greater  average  slope  than  the  number  signal. 

The  contrast  Increases  with  step  depth  and  Is  very  close 
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to  its  maximum  value  of  1.4  when  the  step  depth  is  in¬ 
creased  to  40%  of  the  Bethe  range  of  thp  primary  beam.  Ad¬ 
ditional  step  depth  will  degrade  the  sigual  shape.  A  uni¬ 
versal  curve  of  the  average  signal  slope  as  a  function  of 
the  normalized  horizontal  mark  transition  is  also  presented. 

The  second  paper,  entitled,  "A  Study  of  Alignment  Signals  for  Electron  Lith¬ 
ography,”  by  Lin,  Adesida  and  Neureuther,  will  be  presented  at  the  9th  Interna¬ 
tional  Conference  on  Electron  and  Ion  Beams  Science  and  Technology,  at  St. 
Louis,  Missouri,  in  October  1980  and  is  abstracted  as  follows: 

Abst  ract 

High  quality  alignment  signals  are  required  in  electron 
lithography.  In  this  paper,  the  Monte  Carlo  method  is 
used  to  study  the  alignment  signals  formed  by  electrons 
that  are  backscattered  from  tapered  etched  silicon  marks. 

Effects  of  the  mark  geometry  and  the  incident  beam  voltage 
are  explored.  Universal  curves  of  both  the  contrast  and 
the  average  signal  slope  are  obtained.  The  results  indic¬ 
ate  that  although  a  step  of  depth  greater  than  40%  of  the 
Bethe  range  of  the  primary  beam  yields  a  better  contrast, 
too  large  a  step  will  result  in  the  decrease  of  the  aver¬ 
age  slope.  The  quality  of  the  alignment  signal  can  be  im¬ 
proved  if  the  detector  is  positioned  to  collect  the  back- 
scattered  electrons  with  small  takeoff  angles.  The  sum 
and  difference  signals  from  detectors  in  different  azimuth¬ 
al  angles  are  also  examined.  The  results  show  that  in  de¬ 
tecting  a  topographical  discontinuity,  the  difference  sig¬ 
nal  is  preferred  rather  than  the  sum  signal.  Experimental 
confirmation  of  the  theoretical  predictions  shows  good 
agreement  with  the  universal  curve  of  contrast  versus  the 
normalized  step  depth.  The  alignment  marks  were  fabricated 
by  using  anisotropic  etching  techniques  in  silicon,  and  the 
measurements  were  conducted  in  an  ETEC  Autoscan  SEM. 

Interest  in  the  characterization  of  positive  type  photoresists  has  led  to 
the  consideration  of  diazo  sensitized  novalak  based  photoresist s  as  nega¬ 
tive  electron  beam  resists  using  image  reversal.  Experimental  exposures 
were  carried  out  at  Siemens  Corporation  in  Munich,  Germany,  in  a  collabor¬ 
ative  study  by  W.  G.  Oldham  during  his  Sabbatical  leave.  It  was  found 
that  if  the  photo  active  compound  is  first  selectively  destroyed  in  a 
vacuum  exposure  by  an  electron  beam,  a  subsequent  flood  UV  exposure  in  a 
wet  atmosphere  will  convert  the  remaining  areas  into  a  highly  soluble  state. 
This  in  effect  produces  a  negative  imaging  type  resist  capable  of  submicron 
resolution  at  moderate  doses  (25  uC/ctt/  at  20keV)  without  the  usual  limita¬ 
tions  due  to  swelling  during  development. 

This  work  will  be  published  in  IEEE  Electron  Devices  Letters,  Oct.  1980, 
and  is  entitled,  ”A  High  Resolution  Negative  Electron  Resist  by  Image  Re¬ 
versal,”  by  Oldham  and  Hieke.  An  abstract  of  this  paper  follows: 
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Progress  Since  Last  Major  Proposal,  contfd 

Abstract 

The  use  of  AZ1350  family  photoresists  as  negative  elect¬ 
ron  resists  is  described.  Conventional  photolithographic 
technology  is  used  to  coat  and  process  the  resist,  with 
the  exception  of  an  e-beam  exposure  for  patterning.  A 
flood  UV  exposure  is  used  for  image  reversal.  Using  1.5 
Um  initial  thickness,  the  exposure  threshold  for  6  s  de¬ 
velopment  in  1  :  1  AZ  :  H^O  developer  is  7  uC/cm^.  The 
resist  contrast  under  these  conditions  is  1.3;  and  the 
sensitivity  is  about  25  pC/cro^  (70%  thickness  remaining). 

Useful  resolution  on  SiC^/Si  and  Al/Si02/Si  substrates 
is  demonstrated  to  be  at  least  0.5  pm,  and  the  resist  is 
shown  to  mask  the  plasma  etching  of  Al. 

Continuing  work  on  developing  automatic  systems  for  determining  resist  par¬ 
ameters  advanced  on  two  fronts.  The  exposure  measurement  system  was  revised 
to  include  signal  averaging  to  improve  the  S/N  ratio  and  was  modified  to 
allow  hard  bleaching  to  more  accurately  determine  the  nonbleachable  absorp¬ 
tion.  A  new  approach  for  determining  the  resist  development  rate  by  means 
of  measuring  the  absorption  of  the  remaining  resist  in  situ  was  explored. 
Test  runs  on  an  experimental  apparatus  gave  extremely  noisy  data  apprently 
due  to  the  turbulence  of  the  developer  solution  and  the  absorption  of  the 
dissolved  by-products.  These  experiments  are  described  in  "Characteriza- 
tion  of  Positive  Photoresist s , "  the  MS  thesis  of  S.  Mehotra,  which  is  ab¬ 
stracted  herewith: 


Abstract 

In  this  report,  a  model  for  characterizing  positive 
photoresist  is  discussed.  A  computer-aided  system 
for  resist  measurement  is  described.  Experimental 
techniques  to  determine  the  exposure  parameters  us¬ 
ing  the  model  proposed  by  Dill  and  co-workers  are 
explained.  A  new  approach  for  measuring  the  devel¬ 
opment  rate  as  a  function  of  the  photoactive  active 
inhibitor  concentration  is  proposed. 


Publications 


JSEP  Support 

"Electron  Energy  Dissipation  in  Layered  Media1'  -  PhD  thesis  of 
llosanrai  Adesida.  University  of  California,  1979 

"High  Resolution  Electron-Beam  Lithography  on  Thin  Films"  -  1. 
Adesida,  T.  E.  Everhart,  and  R.  Shimizu.  Jour.  Vac.  Sci.  Tech., 
vol .  16(6),  pp.  1 7 A 3— 1748.  Nov. /Dec.  1979. 

"Monte  Carlo  Simulation  of  Registration  Signals  for  Electron 
Beam  Microfabrication"  -  Y.  C.  Lin,  I.  Adesida  and  A.  R,  Neu- 
reuther.  Appl.  Phys.  Lett.,  vol.  36(8),  pp.  672-674.  1980. 
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Publicat Ions ,  cont*d. 

JSEP  Support,  cont'd. 

"A  Study  of  Alignment  Signals  for  Electron  Lithography,"  Y.  C.  Lin, 
I.  Adesida  and  A.  R.  Neureuther.  Presented  at  the  Electrochemical 
Society  Meeting,  St.  Louis,  Mo.,  May  1980. 

"AZ1350  as  a  Negative  Electron  Resist  by  Image  Reversal,"  E.  Hieke 
and  W.  G.  Oldham,  Microelectronics  *80,  Amsterdam,  Oct.  1980. 

"A  High  Resolution  Negative  Electron  Resist  by  Image  Reversal," 

W.  G.  Oldham  and  E.  Hieke,  to  appear  in  Electron  Devices  Letters. 

"Characterizat ion  of  Positive  Photoresist,"  Sanjay  Mehotra,  MS 
thesis  plan  TT,  University  of  California,  June  1980. 

Related  to  JSEP  Work 


"Direct  Observation  of  Voltage  Barriers  in  Zinc  Oxide  Varistors," 

0.  L.  Krivanek,  P.  Williams,  and  Y.  C.  Lin.  37th  Ann.  Proc.  Elect¬ 
ron  Microscopy  Soc*  Amer. ,  San  Antonio,  Texas,  1979.  G.  W.  Bailey  (ed.) 

"Exploration  of  Electron-Beam  Writing  Strategies  and  Resist  Devel¬ 
opment  Effects,"  M.  G.  Rosenfield  and  A.  R.  Neureuther,  Electrochem¬ 
ical  Society  Meeting,  St.  Louis,  May  1980. 


Significant  Accomplishme nts 

A  highly  accurate  Monte  Carlo  algorithm  for  calculating  electron  scattering 
in  solids  has  been  developed.  A  more  physically  correct  model  which  allows 
electron  scattering  to  be  considered  either  as  individual  inelastic  (direct) 
scattering  events  or  through  the  cont inual-sl owing-down- approxiraat ion  (CSDA) 
has  been  simulated.  Excellent  agreement  with  measurements  of  electron  trans¬ 
mission  through  thin  films  was  obtained.  A  significant  extension  of  the 
Monte  Carlo  program  to  allow  nonplanar  substrates  has  made  it  possible  to 
study  alignment  signals  from  tapered  etched  silicon  steps.  Experimental  and 
theoretical  results  have  shown  that  to  a  surprising  extent  simple  empirical  j 

universal  models  can  be  developed  to  predict  electron  beam  alignment  signals.  i 

I 

A  breakthrough  in  the  problem  of  negative  resist  swelling  and  its  inherent  j 

electron  beam  resolution  has  been  made.  We  have  demonstrated  that  a  non¬ 
swelling  positive  photoresist  (AZ1350J)  may  be  used  as  a  negative  electron 
resist  by  means  of  image  reversal.  Resolution  below  1  yin  in  films  1.2  yin 
thick  was  observed  at  a  dose  of  25  yC/cm^  at  20keV.  The  high  etch  resist¬ 
ance  in  plasma  etching  is  also  retained. 

Interaction  with  Other  Work  Units 

We  have  assisted  Professor  Van  Duzer  in  the  development  of  a  pattern  genera¬ 
tion  system  for  the  ETEC  Autoscan  electron  microscope.  We  have  also  inter¬ 
acted  with  Professor  Hess  of  Chemical  Engineering  in  the  testing  of  PMMA- 
related  resist  materials. 
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Scientific  Objective 

The  objective  of  this  work  is  to  determine  the  limitations  on  miniaturi¬ 
zation  of  Josephson  digital  memory  and  logic  circuits. 

State  of  the  Art 


For  the  past  decade,  work  has  been  in  progress  on  developing  the  use  of 
Josephson  junctions  for  digital  applications.  The  driving  motivation  has  been 
the  extremely  small  delay-power  products  that  have  been  demonstrated.  This  is 
a  result  of  the  fact  that  the  voltages  are  in  the  millivolt  range,  the  currents 
in  the  milliampere  range,  and  the  circuit  switching  times  less  than 
10-10  s.  These  low  energy  levels  have  been  demonstrated  with  rather  large 
circuits.  It  is  evident  that  such  circuits  can  be  made  advantageously  in  much 
smaller  sizes  and  that  thermal  effects  may  not  be  the  limiting  factor. 

There  have  been  some  limited  efforts  at  miniaturization  of  Josephson 
circuits.  Josephson  junctions  as  small  as  1  pm^  have  been  reported1-.  Memory 
cells  with  nondestructive  readout  have  been  reported  with  a  size  of  1.4  mil^, 
which  is  comparable  with  small  present-day  semiconductor  memory  cells^.  The 
report  of  that  memory  cell  also  contained  some  theoretical  guidelines  for  cells 
storing  a  substantial  number  of  flux  quanta.  There  have  also  been  reports  of 
single- flux-quantum  memory  cells  with  destructive  readout^-^  an<j  8ome  details 
of  their  operation  have  been  explained.  The  space  taken  by  these  cells  appears 
to  be  somewhat  less  than  1  mil^.  An  IBM  patent  has  appeared  in  1977  that 
treats  miniaturized  memory  cells  where  the  need  for  inductance  is  satisfied 
by  using  very  thin  films  with  the  attendant  kinetic  Inductance  which  results 
from  the  momentum  of  the  electron  pair  fluid.  There  have  been  no  systemmatlc 
miniaturization  studies  presented  of  the  type  made  by  Keyes?  for  semiconductor 
devices. 


Progress  and  Publications  Since  Last  Major  Proposal 

We  have  not  previously  concentrated  on  the  miniaturization  aspect  of 
Josephson  digital  circuits,  but  have  recently  published  several  papers  on 
various  other  aspects  of  Josephson  digital  circuits. 


With  JSEP  Support 

''Resonances  in  Semiconductor-Barrier  Josephson  Junctions,"  T.  Van  Duzer  and 
W.  Y.  Lum,  J.  Appl.  Phys.,  v.  49,  pp,  4560-4563,  August  1978. 

"Switching  Measurements  on  Josephson  Memory  Loops,"  T.  Van  Duzer,  W.  Y.  Lum 
and  H.  W,  K.  Chan,  J.  Appl.  Phys..  v.  49,  pp.  4302-4303,  July  1978, 

Significant  Accomplishments 


An  analysis  was  made  of  the  smallest  possible  superconductive  cells  that 
can  be  made  using  flux  quantization.  Loops  of  about  one  micron  diameter  formed 
with  102nm  thin  film  lines  and  containing  two  Josephson  Junctions  about  0.3  x 
0.3  Um^.  The  current-transfer  time  to  change  states  in  the  loop  from  a  local 
"0"  to  a  logical  "1"  will  be  approximately  one  picosecond.  A  limiting  factor 
will  be  achievable  current  densities  in  Josephson  Junctions. 
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We  have  studied  the  fabrication  problems  associated  with  making  feature 
sizes  on  the  order  of  0,1  piu-  This  has  resulted  in  devices  for  another  pro¬ 
ject  based  on  a  novel  silicon-coupling  arrangement. 

Interaction  with  Other  Work  Units 


This  program  of  studying  miniaturization  forms  a  part  of  a  larger  effort 
on  various  aspects  of  Josephson  digital  circuits.  A  program  has  been  in¬ 
itiated  to  develop  a  high-speed  A/D  converter  and  this  will  use  Josephson 
junctions  and  circuits  of  conventional  size  that  can  be  made  by  photolitho¬ 
graphy.  In  that  work,  we  expect  to  demonstrate  superconductive  integrated 
circuits  containing  some  tens  of  Josephson  junctions.  Professor  A.  M.  Despain 
of  our  Computer  Science  Division  is  cooperating  in  this  project.  We  also  have 
an  NSF-sponsored  project  to  study  the  use  of  coplanar  junctions  (both  elec¬ 
trodes  are  in  the  same  plane)  with  special  emphasis  on  a  structure  we  developed 
here  that  involves  electron-pair  coupling  through  a  highly-doped  surface  of  a 
semiconductor.  These  may  prove  to  be  valuable  for  miniaturization  and  that 
aspect  will  be  considered. 

Our  work  on  Josephson  devices  and  their  frequent  need  for  submicron 
dimensions  has  given  our  group  a  high  degree  of  familiarity  with  electron 
lithography  and  we  have  published  several  papers  in  that  field.  This 
expertise  will  be  supported  by  the  microfabrication  Interests  of  Professors 
A.  R.  Neureuther  and  W.  G.  Oldham,  as  well  as  the  large  body  of  knowledge 
in  our  integrated  circuits  group  (SSD-81-1). 

POD  Interaction 

We  have  frequent  contact  with  ONK  (E.  A.  Edelsack  and  R.  S.  Brandt)  and 
contact  with  the  DOD  supporters  of  the  IBM  project  on  Josephson  digital 
circ lut s . 
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Scientific  Objective 

To  inveatigate  new  methods  of  realizing  low-noise  amplification  in  mono¬ 
lithic  integrated  circuits.  This  will  involve  work  on  noise  sources  and 
characterization,  and  on  new  techniques  of  circuit  design.  Technologies 
employed  will  include  MOS,  bipolar  and  JFET. 

State  of  the  Art 


Low-noise  amplification  is  a  pervasive  need  in  many  types  of  monolithic 
integrated  circuits,  both  MOS  and  bipolar.  In  MOS  Integrated  circuits 
for  example,  recent  work  has  shown  the  possibility  of  including  analog 
functions  such  as  A/D  and  D/A  conversion  on  MOS  chips  containing  complex 
digital  circuits  such  as  microprocessors.  In  another  area,  it  appears  that 
the  use  of  the  charge  conservation  and  transfer  capability  of  MOS  technology 
will  allow  realization  of  analog  filters  of  low  sensitivity  and  high  order 
in  a  small  die  area.  Such  filters  will  find  wide  application  in  bandwidth 
compression,  speech  recognition  and  communication  systems.  Common  to  all 
these  applications  is  the  need  for  low-noise,  wide -dynamic a range  amplifica¬ 
tion  using  minimum  die  area  Cl}.  In  some  cases  the  realization  may  take 
the  form  of  an  operational  amplifier  or  an  instrumentation  amplifier.  In 
other  cases,  a  fixed-gain  amplifier  (either  differential  or  single-ended) 
may  be  the  optimum  solution.  In  all  cases  there  are  interactions  between 
noise  performance  and  circuit  bandwidth,  offset,  power  consumption  and 
die  area. 

Although  analog  bipolar  circuits  are  generally  in  a  more  mature  state  than 
are  MOS  analog  circuits,  there  is  a  continuing  need  for  new  approaches  to 
the  realization  of  low-noise  bipolar  Integrated  circuits.  For  example,  it 
Is  still  common  to  find  discrete  low-noise  pre-amplifiers  used  external  to 
bipolar  monolithic  circuits  in  critical  applications.  The  need  for  low- 
noise  amplification  is  found  both  in  low-frequency  applications  where  1/f 
noise  dominates,  and  also  in  high-frequency  wideband  amplifiers  [  2]  .  As 
in  the  case  of  MOS  circuits,  there  are  interactions  between  noise  perform¬ 
ance  and  other  critical  circuit  parameters. 

Progress 

CHj r  earlier  work  on  the  fabrication,  characterization  and  modeling  of  gen¬ 
eralized  multi-layer  monolithic  structures  allowed  us  to  realize  a  unique 
monolithic  process  incorporating  high  frequency  JFET  structures  with  inde¬ 
pendently  accessible  gates.  Extensive  characterization  on  a  range  of  de¬ 
vices  has  been  undertaken  and  varified  the  advantages  of  these  structures 
for  a  wide  range  of  high-frequency  monolithic  filters  and  real-time  signal 
processing.  Experimental  high-frequency  filters  using  this  process  are 
currently  in  fabrication.  These  filters  incorporate  very-wide-band  opera¬ 
tional  amplifiers  (unity  gain  beyond  100  MHz)  made  possible  by  the  unique 
JFET  structures  described  previously.  This  project  will  continue  with 
further  work  in  combining  advanced  fabrication  techniques  and  new  system 
concepts  with  the  ultimate  goal  being  the  realization  of  unique  methods 
of  high-speed  signal  processing  using  arbitrary  combinations  of  analog 
and  digital  techniques. 
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A  second  aspect  of  our  work  has  involved  new  techniques  to  achieve  max¬ 
imum  possible  dynamic  range  in  bipolar  and  MOS  analog  amplifiers  and  sig¬ 
nal-processing  circuits.  Our  work  in  this  area  has  focused  initially  on 
methods  of  achieving  maximum  dynamic  range  in  swi tched-capacitor  filters. 
In  such  devices,  the  fundamental  limit  on  dynamic  range  results  from 
thermal  noise  in  the  MOS  switch  transistors  used  in  the  filter  structure, 
hut  monolithic  switched-capacit or  filters  produced  to  date  have  not  ap¬ 
proached  this  limit  because  of  the  dominant  effects  of  flicker  noise  in 
the  MOS  transistors  making  up  the  operational  amplifiers  [1,2].  Our  re¬ 
search  has  resulted  in  a  new  differential  chopper-stabilized  technique 
which  removes  this  1/f  noise,  resulting  in  experimental  fifth-order  filt¬ 
er  prototypes  with  dynamic  range  in  excess  of  105  dB.  This  is  more  than 
L5  dB  higher  than  previously  reported  performance  obtained  from  such 
filters  [3].  We  expect  to  continue  this  work  with  the  objective  of  ex¬ 
ploring  the  maximum  theoretical  dynamic  ranges  achievable,  and  formul- 
nf lug  detailed  theories  of  the  noise  performance  of  such  filters. 

A  c  c  ( )  mp 1  I shme nts 


The  paper  entitled,  "Design  Techniques  for  MOS  Switched  Capacitor  Ladder 
Kilters, M  by  Gordon  M.  Jacobs,  David  J.  Allstot,  Paul  Gray  and  Robert  W. 
Brodersen  received  the  W.  R.  G.  Baker  Prize  Paper  Award.  The  research 
reported  in  this  paper  is  based  in  part  on  this  work  unit.  An  abstract 
of  the  paper  follows: 

Design  techniquest  for  monolithic,  high-precision,  MOS 
sampled-data  active-ladder  filters  are  described.  Switched 
capacitor  integrators  are  used  to  implement  the  "leapfrog" 
configuration  for  simulating  doubly  terminated  LC  ladder 
networks.  Techniques  are  presented  for  designing  all-pole 
low- pass  filters,  as  well  as  r-athods  for  including  trans¬ 
mission  zeros.  An  approach  for  implementing  bandpass  filt¬ 
ers  is  described  which  is  derived  from  the  conventional 
1 ow-pass-to-bnndpass  transportation.  Monolithic  realiza¬ 
tions  for  two  different  low-pass  filters  are  briefly  de¬ 
scribed  which  show  excellent  agreement  with  theory. 

Other  accomplishment s  include  the  following: 

1.  A  new  high-frequency  bipolar  monolithic  process  has  been  devised  and 
extensively  characterized  for  both  low-frequency  and  high-frequency 
parameters .  This  research  has  led  to  new  circuit  configurations  for 
high-frequency  monolithic  filters,  switches  and  amplifiers  which  are 
currently  being  fabricated. 

2.  A  new  differential  chopper-stabilized  swit ched-capacit or  filtering 
technique  was  developed,  and  evaluated  by  fabricating  a  monolithic 
prototype.  This  fifth-order  lowpass  3.5  kHz  filter  displayed  a  dyn¬ 
amic  range  in  excess  of  105  dB. 

Pub  1 Icat ions 

A.  Milder  .ISKP 


l.  C.  M.  Jacobs,  D.  J.  Allstot,  R.  W.  Brodersen  and  P.  R.  Gray, 


Publications,  cont'd. 


A.  Under  JSEP,  cont'd. 

1.  cont'd. 

’’Design  Techniques  for  MOS  Switched  Capacitor  Ladder  Filters/' 
IEEE  Trans,  on  Circuits  and  Systems,  Vol.  CAS-25,  No.  12,  Dec. 
1978,  pp.  1014-1021.  Received  the  W.R.G.  Baker  Prize  Paper 
Award  from  IEEE. 

2 

2.  S.  K.  Lui  and  R.  G.  Meyer,  "High  Frequency  Bipolar- JFET-I  L 
Process,”  submitted  to  International  Electron  Devices  Confer¬ 
ence,  1980. 

3.  K.  C.  Hsieh  and  P.  R.  Gray,  "A  Wide  Dynamic  Range  Differential 
Chopper-Stabilized  Switched-Capacitor  Filter,"  submitted  to 
1981  International  Solid-State  Circuits  Conference. 

B.  Under  Other  Sponsorship 

1.  S.  K.  Lium^R.  G.  Meyer  and  N.  Kwan,  "An  Ion-Implanted  Sub- 
Surface  Monolithic  Zener  Diode,"  IEEE  Journal  Solid-State 
Circuits,  Vol.  SC-14,  No.  4,  August  1979,  pp.  782-784. 

2.  D.  J.  Allstot,  R.  W.  Brodersen  and  P.  R.  Gray,  "An  Electric¬ 
ally  Programmable  Switched-Capacitor  Filter,"  IEEE  Journal 
Solid-State  Circuits,  Vol.  SC-14,  No.  6,  Dec.  1979,  pp-  1034- 
1041. 
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Scientific  Objectives 

We  propose  to  investigate  the  technology  for,  and  the  physical  behavior 
'*  o’ect  ronic  systems  that  are  composed  of  polymer-film  materials  and  inte- 
gi a t ed-cl rcult  devices.  The  objectives  are  to  define  the  materials  and 
interlace  properties  and  to  evaluate  the  technologies  that  may  ultimately 
find  application  to  new  and  improved  transducer  devices. 

Major  task: 

To  investigate  the  control  of  piezoelectric  properties  in  poly-viny lidene 
fluoride  and  poly-acrylonitrile  that  can  be  achieved  through  plasma-initiated 
polymerization.  The  purpose  of  this  effort  will  be  to  develop  a  novel  method 
for  the  production  for  piezoelectric  polymer  thin  films  directly  from  monomer 
sources,  thereby  eliminating  many  of  the  constraints  inherent  in  incorporating 
such  i 1 1ms  into  solid  state  devices.  This  study  is  based  upon  current  work 
in  our  laboratory  which  demonstrates  the  feasibility  of  preparing  films  of 
vinyl  polymers  from  monomer  sources. 

State  of  the  Art 

The  above  task  is  proposed  in  order  to  apply  newly  developed  polymer 
format  ion  technology  (1,  2)  to  poly-viny lidene  fluoride  (and  possibly  poly- 
acr\ loni t ri le) ,  a  polymer  material  that  has  very  interesting  and  useful 
piezoelectric  properties  (3,  4).  The  technique  of  plasma-initiated  polymer¬ 

ization  (1,  2)  has  been  demonstrated  for  certain  vinyl  monomers  in  bulk 
condensed  phases,  but  has  not  been  attempted  for  thin  film  formation.  Our 
experience  with  piezoelectric  films,  particularly  with  ZnO,  has  given  us  a 
background  in  exploitation  of  this  property  in  novel  transducers.  We  expect 
to  show  that  the  plasma-ini t ia ted  polymerization  process  can  permit  tailoring 
cf  the  piezoelectric  activity  of  the  film  in  order  to  optimize  certain  trans¬ 
duce'  behavior. 

Progress  to  Date. 

A  parallel  plate  plasma  initiated  polymerization  reactor  has  been  designed 
and  fabricated.  The  reactor  consists  of  two  three-inch  diameter  parallel  plate 
electrodes  mounted  in  a  vacuum  chamber.  Monomer  vapor  and  inert  gases  enter 
the  reactor  through  the  upper  electrode,  and  the  lower  electrode  supports 
the  substrate  to  be  coated  with  polymer.  A  13.5  mHz  radio  frequency  generator 
is  attached  to  the  system,  and  a  plasma  can  be  generated  between  the  electrodes. 

Sign tficant  Accomplishments 

The  above  equipment  will  be  used  to  test  the  feasibility  of  directly 
forming  piezoelectric  films  from  monomer  vapor.  This  technique  offers  the 
potential  for  eliminating  the  normally  separate  stretching  and  poling  process 

steps . 

Update  of  Current  Period 


A  glass  vacuum  system  has  been  designed  and  constructed.  Tills  system 
will  allow  plasma  Initiated  polymers  to  he  formed  In  the  bulk  phase  (1,  ?). 
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Subsequently,  these  polymers  can  be  applied  to  substrates  during  polymer¬ 
ization  using  spin-on  techniques.  The  bulk  polymerization  represents  an 
important  complimentary  process  to  the  direct  film  formation  technique 
described  above.  In  combination,  the  two  techniques  should  allow  the 
formation  of  polymer  films  that  span  a  wide  range  of  physical  properties 
such  as  molecular  weight  and  crystallinity. 

During  the  coming  year,  preliminary  studies  will  examine  the  manner  in 
which  processing  conditions  affect  the  physical  properties  of  plasma  initiated 
polyvinylidene  fluoride  and  polyacrylonitrile  films.  The  physical  properties 
of  these  films  will  be  correlated  with  the  relative  piezoelectric  strength 
of  these  films.  The  relative  piezoelectric  strength  of  the  films  will  be 
initially  studied  in  our  laboratory  by  measuring  the  voltage  generated 
across  a  film  when  it  is  mechanically  strained  by  another  transducer  (5). 
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About  twenty  years  ago,  the  unity  of  the  basic  concepts  and  techniques  under¬ 
pinning  the  fundamental  study  of  circuits,  control  systems  and  communicat ion 
systems  became  apparent  and  became  known  as  system  theory.  With  rapidly  in¬ 
creasing  power  and  decreasing  cost  of  computing,  these  concepts  and  techniques 
were  applied  to  larger  and  larger  systems.  This  commonality  of  basic  theoret¬ 
ical  tools  and  approaches  still  exists  today  and  as  the  size  the  system  under 
study  increases,  the  nature  of  the  problem  changes. 

These  features  are  particularly  apparent  in  the  first  pair  of  topics  which 
cover  respectively  large  circuits  and  large  control  systems.  Tn  each  of 
these  topics  we  see  a  fundamental  interweaving  of  decomposition  techniques  — 
such  as  diakoptic  analysis,  graphical  decompose  t  ion  by  strongly  connected 
components  and  minimum  essential  sets,  and  decent ralized  control  by  creation 
of  a  hierarchy  — ,  modeling  techniques  —  piecewise  linear  and  apllne  approx¬ 
imations  —  and  analytical  techniques. 

The  third  topic  is  a  bold  project  aimed  at  stepping  out  of  the  by-now-tradi¬ 
tional  use  of  simulation  for  computer-aided  electronic  circuit  design  and 
reaching  for  a  more  sophisticated  level  where  the  simulation  is  coupled  with 
up-to-date  optimization  algorithms  to  direct  the  design  toward  a  minimum  of 
a  previously  chosen  objective  function  or  towards  the  display  of  a  trade-off 
curve  (or  surface),  so  that  the  engineer,  using  additional  data  not  present 
in  the  mathematical  formulation,  can  choose  what  he  deems  to  be  optimum.  A 
significant  feature  of  the  optimization  algorithm  is  that  inequality  con¬ 
straints  may  easily  be  imposed,  a  fact  which  we  consider  crucial  in  engineer¬ 
ing  design. 

The  fourth  topic  addresses  itself  to  signal  processing;  in  particular,  two  di¬ 
mensional  adaptive  signal  processing.  It  is  proposed  to  stimulate  the  syner¬ 
gism  of  competence  in  information  theory,  in  algorithms  and  discrete  trans¬ 
forms,  and  in  LSI  digital  design  to  pitch  tracking  transform  coding  of  speech 
—  by  applying  to  speech  two-dimensional  signal  processing, —  to  the  implemen¬ 
tation  of  the  Winograd  Fourier  transform  algorithm  to  two-dimensional  digital 
filter  design,  to  the  efficient  representation  of  objects  within  two-dimensional 
images,  and  finally  to  the  adaptive  recognition  of  signals. 


To  anyone  who  studies  these  four  proposals  it  is  clear  that  in  the  course  of 
the  work  there  will  be  ample  opportunities  for  interaction  and  cross- fert i 1 1- 
zation,  both  through  normal  informal  contacts  and  through  regular  seminars. 
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Our  proposed  program  in  .‘.uact  electronics  and  optics  emphasizes  the  funda¬ 
mental  problems  of  cb-v:;  o  tor  generation  and  detection  of  radiation  at  wave¬ 
lengths  shorter  than  microwave  .  Much  of  the  work  deals  with  extremely  novel 
devices  whose  potential';  -r-  >ofng  explored:  these  include  the  Josephson 
and  super-Schottky  d  1*  »  :  •*•.•:!  ■;  k  -e  metal-barrier-metal  (MBM)  junction.  Other 

devices  in  the  proposed  *dy ,  such  as  (GaAl)As  lasers  and  high-frequency 
Sc.hottky  diodes,  are  M>a,!  k;.v.vn,  and  the  emphasis  is  on  better  understand- 
ing  of  their  fundament !  :  red  >n  possible  new  applications. 

There  is  a  need  for  bet*,  w  diodes  at  wavelengths  between  10  microns  and 

a  few  millimeters.  Ja  ?'  *.  !\l.imetor  and  submillimeter  ranges,  cryo-devices 

(Josephson  and  super  •  Pola' "  ‘  i-v  junctions)  have  promise  because  of  their  high 
degree  of  nonlinear,  .  y,  '»  «T  'evious  coupling  problems  remain.  At  somewhat 
shel  ter  wavelengths  t.pe bp  •:  d-SOO  nitrons)  Schottky  diodes  appear  promising 
but  as  yet  are  many  orhvs  ot-  magnitude  away  from  their  theoretical  potential. 

At  still  shorter  w.veic, ;,:bo.  met al* barrier-metal  junctions  probably  are  the 
best  potential  rdxi-r  >,  but  ~*ir  promise  is  also  as  yet  largely  unexplored. 

We  feel  that  an  at  :  ar!  o;,  *  fundamental  problem  of  radiation  ln-coupling  is 
in  order,  based  on  the  us  ■  .?•  .n  pumas,  novel  waveguides,  and  mixer  arrays.  The 
devices  mentioned  .ih-'vc  i  *>  lend  themselves  to  these  techniques,  and  if 

better  coupling  car.  h.  =.  :p  1 ished  the  potential  of  the  devices  can  be  re¬ 
al  .1  r  e  d . 

The  metal-barrier  r  :  interesting  because  it  probably  has  the 

highest  frequency  «  ;  '  •  .  ur  v  known  lumped  circuit  element.  It  is  be¬ 
lieved  that  it  operate!.  v  •  -*  ;io  of  nc  currents  at  frequency  as  high  as  t  he 

visible.  This  env  i  :  look  for  novel  ways  to  exploit  such  properties. 

One  part,  of  the  pi  -'r-  .  -  :*t  has  to  do  with  the  addition  of  a  second  ’’collec¬ 
tor"  junction  so  as  1 o  mlm-  whether  a  t  rans  is  tor-  like  throe-element  de¬ 

vice  is  possible.  Th<s  vi'iv  sb-'  involves  fundamental  studies  of  the  MBM 
J  imct  ion  and  its  r^-chan  tar: 

Tim  sum  i  conductor  l.m  !  *  air.se  a  well-known  device  of  great  importance, 

p  ingle  -mode  op^ra;  .  <;  :  .  .•  ;*  still  more  useful,  and  much  lias  been 

accomplished  her:*  ■  •!  but  the  fundamental  physical  mechanisms 

of  mod-.*  .n  ■■  ?*  .m-:  clear,  and  studies  of  them  will  continue, 

with  emphasis  on  carr  b  ;  M  t  .  f  n.  The  study  of  coupled  oscillators, 
f.  r  use  as  phased  .mrm  radiation,  will  continue. 

The  research  program  ;»•  i  s.-d  is  c  lose!  y  related  with  work  under  other 

.sponsorship  and  tm  t b  •  v.  •.*  iiei  units.  (1)  The  mixer-array  studios 
will  involve  u*o  el  -  ■  lubricated  under  other  sponsorship.  Pro¬ 

gress  with,  the  ca.M-'iirr  ,  i  •,  v  will  provide  guidance  tor  the  fabrication 
work.  Thorp  ;«  ,  *t..-i;hip  with  the  e  1  ec  t  roinagnot  i  c  theory  unit 
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of  this  proposal;  certain  fundamental  problems  of  electromagnetics  have  been 
raised  by  the  problem  of  coupling  to  mixers ,  and  their  solution  will  contri¬ 
bute  to  achievement  of  the  goals  proposed  here.  There  is  also  considerable 
interest  in  millimeter/infrared  heterodyne  detection  in  the  Departments  of 
Physics  and  Astronomy  at  this  university.  We  have  been  maintaining  contact 
and  plan  to  continue,  as  this  provides  guidance  as  to  the  state  of  the  art 
and  as  to  what  improvements  are  required.  Considerable  work  on  development 
and  fabrication  of  semiconductor  lasers  is  being  carried  on  in  our  laboratory 
under  other  support,  and  techniques  made  available  by  that  work  will  facili¬ 
tate  the  more  fundamental  studies  and  experimental  in  this  proposal.  There  is 
also  close  cooperation  between  the  two  projects  in  Quantum  Electronics  and 
Optics . 
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Scientific  Ojectlves 

We  propose  to  explore  and  to  study  novel  and  complex 
semiconductor-laser  structures  for  tunability  and  traveling-wave  opera¬ 
tion,  Theoretical  analysis  based  on  scattering-matrix  formalism  will  be 
carried  out  to  gain  insight  on  possible  modes  of  oscillation  in  complex 
laser  structures  and  to  guide  us  in  the  search  for  new  laser  structures. 
Experimental  techniques  using  LPE  and  MBE  will  be  developed  to  fabricate 
the  various  structures.  Our  objectives  are  (1)  to  advance  the 
semiconductor-laser  technology  beyond  its  present  state  as  an  individual 
coherent  optical  source,  (2)  to  explore  new  modes  of  operation  for  a 
semiconductor  laser  other  than  the  simple  Fabry-Perot  mode,  and  (3)  to 
study  laser  performance  when  interacting  with  other  lasers  or  an 
integrated-optical  network . 

State  of  the  Art 

In  recent  years  important  advances  in  semiconductor  lasers  have 
been  reported  in  the  attainment  of  linear  light-output  versus  input- 
current  relation,  stable  transverse  and  lateral  mode  pattern,  and  single 
longitudinal  mode  operation.  Equally  impressive  are  achievements  in 
guided-wave  optical  circuit  elements  such  as  directional  couplers 
switches  and  modulators.  However,  except  for  some  interesting  work  on 
arrays  [1,2]  most  of  the  research  efforts  on  semiconductor  lasers  and 
waveguide  optical  elements  have  been  limited  to  the  performance  and 
characteristics  of  individual  devices.  To  expand  the  usefulness  of 
these  devices  and  to  explore  the  potential  as  well  as  limitations  of 
integrated  optics,  these  devices  must  be  studied  in  the  frame  of  an  opt¬ 
ical  network.  We  believe  that  it  is  opportune  to  apply  the  knowledge 
accumulated  and  the  techniques  developed  for  the  microwaves  to  the 
advancement  of  guided-wave  optical  systems  involving  sources  and  other 
optical  elements. 

One  major  advance  in  semiconductor-laser  technology  is  the  incor¬ 
poration  of  two-dimensional  waveguiding  in  the  transverse  plane  of  the 
laser  structure.  During  the  past  few  years #  many  lasers  exhibiting 
stable  transverse  and  lateral  mode  pattern  and  single  longitudinal  mode 
operation  have  been  developed  at  several  laboratories.  These  include 
the  buried-heterostructure  (BH)  laser  [ 3 *  ^ »  the  channeled-substrate- 
planar  (CSP)  laser  [5,6],  and  the  laser  utilizing  thickness  variation  in 
the  grown  layer  to  provide  waveguiding  in  the  lateral  direction  [7-10], 
Among  the  lasers  referenced  here  are  the  1RW  laser  [1],  the  CJSF  laser 
[7]  and  the  narrow  DCC-CSP  laser  [6]  being  developed  and  studied  in  our 
laboratory.  All  the  lasers  mentioned  above  oscillate  in  well  defined 
cavity  modes  in  all  three  directions,  transvers,  lateral,  and  longitudi¬ 
nal.  This  fact  facilitates  the  application  of  waveguide-network  theory 
to  semiconductor  lasers  and  makes  possible  an  analysis  of  laser  struc¬ 
tures  more  complex  than  the  simple  Fabry-Perot  structure  which  has  been 
almost  exclusively  employed  in  the  semiconductor  lasers  thus  far. 

In  the  meantime,  remarkable  progress  has  also  been  made  in  the 
fabrication  and  control  of  modes  in  optical  waveguides.  Low-loss, 
single-mode  waveguides  have  been  reported  in  GaAs  [11]  and  LiNbO^  [12], 
and  active  waveguide  devices  using  coupled  guides  and  Y  junctions  either 
alone  or  in  the  Mach-Zehnder  interferometer  conf iguration  have  been 


developed.  These  advances  further  support  our  view  that  single-mode 
guided-wave  technology  is  ready  for  the  exploration  of  new  and  complex 
semiconductor-laser  structures  and  the  development  of  simple  Integrated 
optical  systems  involving  the  source  and  other  optical  elements.  On  the 
other  hand,  when  semiconductor  lasers  are  used  in  practical  systems  such 
as  in  conjunction  with  optical  fibers,  practical  problems,  such  as 
induced  pulsation  in  the  laser  output,  do  arise.  The  problems  caused  by 
feedback  are  well  known  in  network  theory.  Therefore,  it  seems  timely 
that  we  utilize  the  vast  amount  of  knowledge  in  waveguide  network  theory 
in  advancing  the  state  of  semiconductor  lasers  beyond  its  present  state 
as  individual  sources  and  in  analyzing  and  overcoming  problems  encoun¬ 
tered  in  practical  systems. 

Much  of  the  theory  of  microwave  networks  [13]  is  directly  applica¬ 
ble,  since  both  microwave  and  optical  networks  are  distributed  elec¬ 
tromagnetic  systems  with  waveguide  modes  acting  as  input  and  output 
ports.  The  known  limitations  on  scattering-matrix  parameters  [131  must 
be  considered  as  constraints  on  the  system.  Optical  parallel  to  most  of 

the  passive  microwave-network  elements  -  hybrids,  filters,  isolators, 

power  dividers  -  have  been  worked  out  in  integrated  optic  form,  but 

not  all  in  a  form  directly  integrable  with  the  semiconductor  lasers. 
Moreover,  new  fabrication  techniques  such  as  MBE  will  permit  still  dif¬ 
ferent  configurations  of  these  elements.  Elegant  theoretical  analyses 
of  the  behavior  of  certain  specific  junctions  have  been  made  [1*4]*  but 
again  not  for  all  that  are  appropriate  to  the  specific  conf igurations  or 
fabrication  techniques.  Thus  much  remains  to  be  done  on  the  network 
aspect  of  this  problem. 

Progress  to  Date 

During  the  period  October  1978  -  September  1980,  a  major  effort  was 
devoted  to  a  study  of  the  mode-stability  problem  in  semiconductor  lasers 
and  the  development  of  a  mode-stabilized  laser,  the  narrow  DCC-CSP 
laser.  The  idea  behind  the  laser  can  be  summarized  as  follows.  If  the 
width  of  active  region  is  made  small  compared  to  the  carrier  diffusion- 
length,  the  effect  of  hole  burning  (which  is  a  major  cause  for  lateral 
mode  instability)  caused  by  stimulated  omission  can  be  greatly  reduced 
through  carrier  diffusion.  The  approaches  taken  in  implementing  the 
idea  are  (1)  the  use  of  a  narrow  channel  (?  pm)  smaller  than  the  carrier 
diffusion  length  and  (?)  the  use  of  a  double-current-confinement  (DCC) 
scheme  to  confine  the  injection  current  to  the  narrow  channel  region.  A 
series  of  publications  have  resulted  from  this  investigation  which  are 
listed  under  Publications.  The  outstanding  features  of  the  laser  are 
presented  under  Significant  Accomplishments.  Now  that  we  have  acquired 
an  understand ing  of  the  mode-instability  problem  and  have  developed  a 
laser  with  excellent  mode  stability  which  oscillates  in  a  well-defined 
cavity  mode,  we  can  move  to  a  new  and  exciting  phase  of  our  research  on 
combining  semi  conductor-] aser  technology  and  waveguide-network  tech¬ 
niques  (developed  for  microwave)  for  exploration  of  this  relatively 
unexplored  area. 

A  theoretical  study  has  been  initiated  on  the  proper  modes  of 
oscillation  for  two  complex  laser  structures:  one  with  a  branch 
waveguide,  and  the  other  with  a  ring  waveguide,  attached  to  a  straight 
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Fabry-Perot  cavity.  The  former  structure  is  to  provide  longitudinal¬ 
mode  discrimination  and  hence  wavelength  selectivity  for  the  laser  while 
the  latter  structure  is  to  support  traveling-wave  lasing  modes  in  the 
laser  in  contrast  to  standing-wave  modes  in  the  conventional  Fabry-Perot 
laser.  Scattering-matrix  formalism  is  used  to  relate  the  inputs  to  the 
outputs  of  the  laser  structures.  The  proper  longitudinal  mode  of  oscil¬ 
lation  can  be  found  by  setting  the  determinant  of  the  scattering  matrix 
to  zero  (meaning  nonzero  outputs  with  zero  inputs). 

Analytical  results  on  the  wavelength-selective  laser  are  very 
encouraging.  The  oscillation  condition  for  such  a  laser  is  given  by 


?  2  2  2 
R4"P i  (p 2T2  +  P3T3)  = 


(1) 


where  R  is  reflection  coefficient  at  the  three  exit  facets  (two  for  the 
Fabry-Perot  cavity  and  one  for  the  branch  waveguide),  P  is  the  appropri¬ 
ate  propagation  factor  in  each  segment  of  the  composite  laser  structure, 
and  T2  and  T 3  represent  the  transmission  coefficient  for  each  branch  at 
the  waveguide  junction.  The  above  equation  can  be  re-arranged  into  the 
familiar  form  PjP*j  =  1  with  an  effective  reflection  coefficient 


R2(T2  +  t3p2/p2) 


(2) 


Note  that  the  propagation  factors  P3  and  p2  contain  both  magnitude 
and  phase.  The  quantity  Reff  reaches  a  maximum  when  the  two  phase  fac¬ 
tors  are  in  phase  and  is  a  minimum  when  the  two  phase  factors  are  out  of 
phase.  Therefore,  the  interference  of  the  two  fields  at  the  waveguide 
junction  effectively  modulates  the  threshold  condition  of  the  laser,  and 
thus  provides  wavelength  selection  for  the  laser.  We  should  point  out 
that  wavelength  selectivity  is  lost  if  P3  =  P p,  For  wavelength  selec¬ 
tivity  to  be  most  effective,  the  laser  structure  should  be  so  designed 
that  p2  and  P3  are  in  phase  only  at  one  wavelength  over  the  spectral 
width  of  the  gain  profile. 

Experimental  work  on  the  wavelength-selective  laser  is  under 
preparation.  The  first  structure  to  be  tested  will  be  one  with  the 
branch  waveguide  made  of  a  quarter-ring  waveguide.  The  circular 
waveguide  is  chosen  in  order  to  make  the  two  parts  asymmetrical  so  that 
a  phase  P2  not  identical  to  phase  P3  is  guaranteed.  The  analytical  work 
on  the  traveling-wave  laser  is  somewhat  more  complicated.  Preliminary 
results  indicate  that  the  laser  output  can  can  be  made  to  contain 
predominantly  one  traveling-wave  component  of  the  ring  cavity,  either 
clockwise  or  counterclockwise .  More  work  is  needed  to  assure  the 
correctness  of  the  results  and  to  put  the  results  in  a  form  most  easily 
interpretable . 
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Progress  to  Date,  c o n t 1 d 

During  the  past  year,  we  have  also  investigated  the  metal-oxide  barrier 
in  its  potential  for  an  optical  amplifier.  This  has  included  analysis  of  the 
MOMOM  structure,  studies  of  coupling  into  the  tunneling  region  through  the 
surface  wave  of  the  metallic  portion,  and  development  of  an  "edge-diode" 
configuration  which  gave  an  Interaction  area  or  the  same  order  as  a  point 
contact  diode  (10“-^  cnr)  but  with  the  stability  of  integrated-optics  con¬ 
figurations.  The  conf igurat ion  has  since  been  used  to  advantage  in  small 
superconductivity  configurations.  Pinholes  were  a  problem  in  the  thin  films 
needed  for  the  MOMOM  dev  let1,  bet  it  is  believed  that  the  new  molecular-beam 
epitaxy  unit  will  minimize  these.  There  are  many  more  fundamental  problems 
before  it  can  be  determined  it  optical  amplification  in  this  device  is  a 
possibility,  but  as  these  are  also  applicable  to  the  detection  and  mixing 
problem,  they  will  be  carried  out:  in  the  next  work  period  in  the  work  re¬ 
lated  to  met al-bar rier-met al  detectors  and  mixers. 

Publicati ons 

A,  Under  JSEP 

L.  Figueroa,  ”GaAs  Double-He terostructure  Injection  Lasers  with 
Lateral  Optical  and  Current  Confinement,”  Ph.D.  thesis,  University 
of  California,  Berkeley,  November  1978. 

C.  Y.  Chen,  "Control  of  Lateral  and  Longitudinal  Modes  of  the 
GaAs/Ga i_xAlx As  Semiconductor  Laser,”  Ph.D.  thesis.  University  of 
California,  Berkeley,  September  1980. 

B.  Under  Combined  Sponsorship 

C.  Y.  Chen  and  S.  Wang,  ’’Narrow  Double-Current-Confinement 
Channeled-Substr ate-Pl  ana’"  Laser  Fabricated  by  Double  Etching  Tech¬ 
nique,”  paper  16.6  Conference  Proceedings,  Post  deadline  papers, 
Second  International  Conference  on  Integrated  Optics  and  Optical 
Fiber  Commun icaci on ,  Amsterdam,  September  17-19,  1979. 

C.  Y.  Chen  and  G.  Wang,  "Narrow-Channel  Double-Current-Confinement 
CSP  Laser  and  Effect  of  the  Current  Distribution  on  Its  Perfor¬ 
mance,”  paper  MCI,  Technical  Digest,  Topical  Meeting  on  Integrated 
and  Guided-Wave  Optics,  Incline  Village,  Nevada,  January  28-30, 

1980. 

C.  Y.  Chen  and  S.  Wang,  "Double  Ftching  Technique  for  the  Fabrica¬ 
tion  of  Submicron  Charnels  on  a  GaAs  Wafer  and  Its  Application  to 
Laser  Fabrication,”  Jour.  Appl .  Phys.,  vol.  SI,  pp.  1802-1808, 

1980. 

C.  Y.  Chen  and  S.  Win*.  "Narrow  Double-Cur ren t-Con f inement 
Channeled-Substrate-Pl anar  laser  Fabricated  by  Double  Etching  Tech¬ 
nique,”  Appl  .  Phys.  Lett.,  vM  .  96,  pp .  69*1-696,  1980. 
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Publications,  cont'd 

B.  Under  Combined  Sponsorship,  cont'd. 


C.  Y.  Chen  and  S.  Wang,  "Near-Field  and  Beam-Waist  Position  of  the 
Semiconductor  Laser  with  a  Channeled-Substrate  Planar  Structure," 
Appl .  Phys.  Lett. ,  vol .  37,  pp.  257-260,  1980. 

Invited  paper,  S.  Wang,  C.  Y.  Chen,  A.  S.  H.  Liao,  and  L.  Figueroa, 
"Control  of  Modal  Behavior  in  Semiconductor  Lasers,"  under  prepara¬ 
tion  for  Jour .  Quantum  Electronics . 

Significant  Accomplishments 

We  have  developed  a  double-etching  technique  (DET)  to  make  grooves 
of  micron  and  submicron  dimensions  using  conventional  photolithography. 
Grooves  having  channel  width  as  small  as  0.8  jjm  have  been  made  with  uni¬ 
form  quality  and  the  results  are  reproducible.  Using  the  technique,  we 
have  developed  a  new  laser  structure,  the  narrow  DCC-CSP  laser  in  which 
the  active  region  is  confined  to  a  half-width  comparable  to  or  smaller 
than  the  diffusion  length.  By  so  doing  we  minimize  the  effect  of  hole 
burning  and  by  utilizing  the  built-in  index  difference  provided  by  the 
differential  loss  in  the  tsr  laser,  we  ntnre  ucen  able  to  achieve  excel¬ 
lent  lateral-mode  stability.  The  narrow  DCC-CSP  laser  shows  linear 
light-current  relation  and  stable  fundamental  mode  in  both  transverse 
and  lateral  directions  up  to  51th  of  catastrophic  damage)  and 

single  longitudinal  mode  up  to  2.21^  (the  highest  reported  in  the 
literature).  An  experiment  on  a  new  laser  structure,  the  TEPOSE  laser 
with  two  electrodes  for  pumping  and  one  electrode  for  stimulated  emis¬ 
sion  is  devised  to  study  the  effect  of  nonuniform  and  asymmetric  current 
distribution  as  a  result  of  unequal  excitation  of  the  electrodes.  It  is 
found  that  asymmetric  current  distribution  can  degrade  the  fundamental 
lateral-mode  behavior  and  nonuniform  current  distribution  can  lead  to 
multiple  longitudinal-mode  operation.  This  study  demonstrates  the 
importance  of  uniform  and  symmetric  excitation  as  a  controlling  element 
to  mode  stability  in  addition  to  built-in  index  difference.  An  experi¬ 
ment  is  also  established  for  measuring  the  beam-waist  width  and  the 
beam-waist  position  to  ascertain  the  guiding  mechanism  in  a  semiconduc¬ 
tor  laser  and  to  observe  the  stability  of  the  laser  characteristics  as  a 
function  of  the  pumping  current.  We  have  been  asked  by  the  Editor  of 
Jour.  Quantum  Electronics  to  write  an  invited  paper  on  the  subject  of 
control  of  modal  behavior  in  semiconductor  lasers. 

The  edge-diode  configuration  developed  for  the  MOM  devices  has  been 
widely  used  in  these  and  in  Josephson  junction  devices. 

Interaction  with  Other  Work  Unit 

Frequent  consultation  is  anticipated  with  the  work  unit  on  elec¬ 
tromagnetics  concerning  wave  propagation  in  dielectric  waveguides. 
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Scientific  Objective 

In  this  project  we  propose  to  study  radiation  detection  systems 
involving  direct  detection  and  superheterodyne  mixing  techniques  using 
various  nonlinear  devices  for  mixers  and  detectors.  The  wavelength 
range  of  interest  is  the  entire  band  of  frequencies  from  millimeter 
through  far  infrared,  that  is  10  mm  to  10  microns  wavelength.  It  is  of 

particular  interest  to  understand  which  techniques  are  favored  -  in 

the  framework  of  realistic  receiver  systems  -  for  application  in  vari¬ 
ous  parts  of  the  frequency  range.  The  investigation  will  involve  cou¬ 
pling  of  radiation  to  devices  with  impedance  matching,  selection  of  the 
best  devices  and  obtaining  their  characteristic  parameters,  and  a  study 
of  optimal  coupling  to  intermediate-frequency  amplifiers.  Objectives 
include  low-noise  receiver  performance,  plus  other  advantages  such  as 
wideband  response  and,  for  some  devices,  room-temperature  operation. 
The  novel  devices  in  question  (being  developed  partially  under  other 
sponsorship)  include  planar  Schottky  diodes,  super-Schottky  diodes, 
cryogenic  SIS  and  Josephson  junctions,  metal-barrier-metal  diodes,  and 
novel  configurations  of  other,  more  conventional  detectors.  Attention 
will  be  given  to  principles  of  coupling  to  multiple-device  arrays. 

State  of  the  Art 

The  noise  temperatures  of  mixer-IF  amplifier  combinations  have  been 
surveyed  for  the  frequency  range  1  GHz  (wavelengths  of  30  cm  to  0.3  cm) 

[1] .  Room  temperature  mixers  are  usually  Schottky  diodes  in  this  range 
and  their  best  noise  properties  are  typified  by  a  measurement  by  Kerr 

[2]  at  100  GHz;  assuming  an  IF  amplifier  with  a  15  K  noise  temperature, 
the  single-sideband  noise  temperature  for  the  mixer-IF  amplifier  combi¬ 
nation  is  about  550  K.  Equivalently,  this  is  a  minimum  detectable  power 
of  0.76  x  10~20  w/Hz.  That  value  was  then  further  reduced  by  nearly  a 
factor  of  two  by  cooling  the  diode  (essentially  the  same  results  are 
achieved  at  T  r  18  K  or  T  =  77  K). 

To  achieve  large  reductions  of  noise  temperature  for  frequencies  in 
the  several-millimeter  wavelength  range,  superconductive  devices  must  be 
used.  There  are  two  general  categories.  In  one,  the  Josephson  effects 
are  employed  and  in  the  other,  the  nonlinearity  that  comes  with  single¬ 
particle  tunneling  in  superconductors.  The  best  results  obtained  with 
the  Josephson  devices  have  mixer  noise  temperatures  in  the  range  of  2 0- 
50  times  the  liquid-helium  bath  temperature  with  conversion  efficiencies 
in  the  range  of  0.5-1. 4.  Mixer  noise  temperature  of  54  K  with  conver¬ 
sion  efficiency  of  1.4  at  36  GHz  [31,  and  140  K  with  conversion  gain  of 
2.4  at  140  GHz  [4]  have  been  achieved.  Opinion  now  is  that  the  Joseph¬ 
son  devices  have  too  much  noise  down-converted  into  the  IF  range  by 
beating  of  the  RF  noise  with  the  many  harmonics  of  the  L0  signal  and  the 
Josephson  oscillations.  Still,  Taur  estimates  from  detailed  computer 
calculations  that  a  receiver  noise  temperature  of  70  K  should  be  achiev¬ 
able  for  a  SS8  receiver  at  100  GHz,  assuming  an  IF  amplifier  with  a  20  K 
noise  temperature  [5]. 

Considerable  emphasis  is  now  being  placed  on  the  "quasiparticle” 
(single  particle)  tunneling  devices,  also  known  as  NIS  and  SIS  junc¬ 
tions.  This  has  grown  out  of  the  realization  that  very  strong  non- 
linearities  of  the  I-V  characteristic  are  available  for  classical 
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resistive  mixing  through  the  use  of  devices  that  function  at  millivolt 
levels  where  a  sharp  change  of  slope  results  from  the  existence  of  the 
superconducting  energy  gap.  The  experimental  work  has  been  further 
inspired  by  the  quantum-mechanical  analysis  of  Tucker  [6]  showing  that  a 
resistive  mixer  can  be  expected  to  give  quantum-limited  detection  sensi¬ 
tivities  at  high  frequencies  (>70  GHz).  It  has  also  shown  that  conver¬ 
sion  gain  can  be  expected  from  these  devices.  The  earliest  of  the 
quasiparticle  devices  was  the  super-Schottky  diode,  which  is  a  Schottky 
diode  with  the  metal  electrode  replaced  by  a  superconductor.  These 
showed  the  best  diode  noise  temperature  ever  seen  at  10  GHz  [8].  In  the 
usual  configuration,  the  series  resistance  leads  to  a  parasitic  conver¬ 
sion  loss  which  increases  strongly  with  frequency  so  that  usefulness  of 
the  device  is  limited  to  wavelengths  longer  than  several  millimeters. 
In  the  normal  metal-insul ator-superconductor  (NIS)  tunneling  device,  the 
Schottky  barrier  is  replaced  by  that  of  the  insulator,  and  the  series 
resistance  is  largely  eliminated.  The  same  applies  in  the  case  of  the 
superconductor-insulator-superconductor  (SIS)  junction;  in  this  device, 
the  nonlinearity  is  stronger  but  it  also  has  Josephson  current  which  can 
add  noise.  Results  to  date  on  the  SIS  and  NIS  junctions  have  been 
excellent  [9].  Experiments  with  SIS  junctions  at  36  GHz  have  demon¬ 
strated  the  predicted  conversion  efficiency.  The  single-sideband  value 
of  0.7  (or  1.4  in  DSP)  is  greater  than  can  be  expected  from  classical 
mixing  theory.  The  mixer  noise  temperature  of  3  K  (SSB)  is  very  close 
to  the  photon  noise  limit  hv/kp  =  1.7.  These  experiments  were  done  with 
lead-alloy  oxide-barrier- junctions .  It  would  be  advantageous  to  use  a 
more  rugged  system  such  as  niobium  on  thin  silicon  membranes. 

The  highest  frequency  at  which  the  quasi  particle  mixers  have  been 
reported  is  115  GHz  where  a  conversion  efficiency  of  0.2  and  a  mixer 
noise  temperature  of  less  than  about  80  K  were  deduced  from  a  receiver 
noise  temperature  of  470  K.  One  part  of  our  work  is  aimed  at  making  an 
antenna-coupled  mixer  at  one  millimeter  wavelength;  the  first  measure¬ 
ments  shall  be  at  70  GHz.  Carrying  the  advantages  of  these  devices  to 
the  submi 1 1 imeter  wavelength  range  is  a  large  and  important  task. 

At  wavelengths  near  10  microns,  extrinsic  and  intrinsic  photoeon- 
due ti vo  and  photovoltaic  detectors  can  be  used  as  mixers.  Minimum 
detectable  power  near  the  quantum  limit  ( about.  10-00  Watt  for  1  Hz 
bandwidth)  has  been  demonstrated  [ 9 . 1 ^ 1 .  Although  results  have  been 
impressive,  it  is  difficult  to  obtain  IF  bandwidths  much  larger  than  1 
GHz  with,  such  mixers.  Moreover,  they  all  require  refrigeration,  at 
least  to  77°  K.  As  one  goes  toward  longer  wavelengths,  mixers  are  more 
difficult  to  make.  A  detectable  power  of  3  x  10~^  W  at  119  i-J  using  an 
extrinsic  Si  mixer  and  ?°  K  has  bbeen  reported  [11];  in  another  experi¬ 
ment  with  an  extrinsic  Ga As  mixer  at  4°  K,  a  detectable  power  of  10“^4 
watt  was  reported  at  3a?  j)  [12]. 

In  principle,  one  should  be  able  to  do  much  better  t.han  this  at 
wavelengths  around  0.1  mm,  using  resistive  diode  mixers  at  room  tempera¬ 
ture.  Minimum  deter  table  powers  on  the  order  of*  kTH,  or  4  x  10“^ 
watt/Hz,  are  theoretically  possible.  Rut  results  achieved  to  date  fall 
short  of  theoretical  limits:  with  Schott ky  diodes,  for  example,  room- 
temperature  detectable  powers  of  10”°,  and  r>  x  10"^'  W/Hz  at  118, 
170  and  83'7  microns,  respectively  [13,141.  One  experiment  has  also  been 
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done  to  test  the  metal-barrier-metal  diode  as  a  mixer  [1*4]  at  10 
microns;  the  detectable  power  was  w/Hz  at  room  temperature.  This 
is  the  lowest  value  that  has  been  obtained  for  10-micron  t oom- 
temperature  detection,  but  it  also  far  above  theoretical  limits. 

With  regard  to  room-temperature  detection  and  mixing  in  the  1-3  mm 
range,  conventional  Schottky  diodes  have  given  good  performance.  How¬ 
ever,  there  are  difficulties  with  stability,  ruggedness,  and  cost. 
Presently  a  new  technology  seems  to  be  possible,  based  on  the  integra¬ 
tion  of  dielectric  waveguides  and  planar  (or  "surface  or iented,,M ) 
Schottky  diodes,  and  perhaps  other  components  into  man-wave  integrated 
circuits.  Under  other  support  we  are  presently  studying  devices  for 
this  purpose.  Our  work  under  JSF.P  relates  to  the  principles  and  poten¬ 
tial  capabilities  of  such  systems,  especially  as  compared  with  the  other 
types  of  receivers  already  mentioned. 

In  summary,  progress  in  infrared/man  detection  and  mixing  in  the 
coming  years  will  involve  (a)  extension  of  the  nearly  quantum-limited 
sensitivities  of  superconductive  devices  as  far  as  possible  into  the 
range  above  100  GHz;  (b)  extension  of  the  capabilities  of  room- 
temperature  devices  to  wavelengths  below  one  millimeter;  and  (c) 
development  of  efficient  coupling  schemes  and  rugged  devices. 

Progress  and  Publications  Since  Last  Major  Proposal 

Work  has  continued  on  silicon  membrane,  sandwich-type  superconduc¬ 
tive  devices  which  can  be  made  either  as  Josephson  junctions  or  super- 
Schottky  diodes.  These  are  being  made  with  a  two-step  etching  pro¬ 
cedure.  First,  a  membrane  of  about  50  x  50  and  3  pm  thick  is  made 
by  preferential  etching.  The  boron  dopant  used  to  stop  the  etching  is 
then  removed.  A  new,  shallow  doping  is  then  done,  and  a  small  etch  pit 
is  made  in  the  3  JJm  thick  layer.  The  final  membrane  is  to  be  about  1 
pn^  with  a  thickness  of  50  -  100  nm .  The  reason  for  using  the  two-step 
process  is  that  it  is  necessary  to  know  the  thickness  with  greater  pre¬ 
cision  than  the  allowed  errors  in  the  dimensions  of  the  final  membrane. 
We  can  measure  the  =3  pm-thick  layer  with  some  precision.  We  have  had 
difficulty  with  spurious  etching  during  the  out-diffusion  step  and  con¬ 
siderable  effort  has  been  devoted  to  developing  an  understanding  of  the 
process.  Some  improvements  have  been  made  in  the  yield.  We  have  also 
taken  the  approach  of  doing  the  first  etching  step  without  a  prior  dop¬ 
ing;  the  etching  control  is  by  timing  and  visual  observation.  The  mem¬ 
branes  thus  formed  do  not  have  the  damage  which  is  a  natural  ooncommi- 
tant  of  the  heavy  doping  for  the  thick  membrane;  we  will  seek  to  deter¬ 
mine  whether  that  changes  the  performance  of  the  device.  We  have  also 
developed  a  procedure  for  nondestructive  measurement  of  the  membrane 
thickness  by  using  the  penetration  of  the  electron  beam  in  a  scanning 
electron  microscope.  We  have  started  to  develop  a  model  for  the 
semiconductor-bar r ier  junction.  A  00  GHz  microwave  mixing  system  has 
been  prepared  for  testing  the  devices  both  as  Josephson  junctions  and  as 
super-Schot tky  diodes.  This  work  is  continuing. 

Quasipartiele  SIS  mixers  for  the  millimeter  range  are  also  being 
developed.  Here  our  approach  havS  centered  on  unifying  the  STS  junction 
with  an  antenna  structure,  the  latter  to  provide  efficient  coupling.  A 
typical  design  consists  of  a  sandwloh-t  yp«*  vee  antenna  with  four  SIT. 
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junctions  in  series  at  its  termin'' 1  , .  Performance  approaching  theoretic¬ 
al  quantum  limits  is  exp**'*  1  ,iear  70  GHz.  At  higher  frequencies,  the  cap¬ 

acitance  of  the  junctions  becomes  a  problem,  and  a  novel  method  has  been 
developed  to  tune  reactance  out  This  method  involves  the  use  of  an  elect¬ 
rically  adjustable  reactance  used  to  cancel  this  SIS  junction's  capacitance. 
The  adjustable  reactance  is  a  large  Josephson  junction,  dc  biased  below  the 
critical  current  and  connected  across  the  SIS  junctions  with  a  length  of 
superconducting  transmission  line.  The  latter  is  used  to  transform  the  ad¬ 
justable  reactance  into  tin*  proper  magnitude  and  phase  for  cancellation. 

With  this  scheme  it  should  he  ossible  to  use  the  antenna-coupled  SIS  at 
frequencies  up  to  300  GHz. 


In  order  to  te'f  •  higher  frequency,  a  23 0  GHz 
isotopic  methyl  fluor:  :•  *•••'*'■  ('instructed.  The  COp  pump  laser 
for  it  is  now  working  w  *  •*  1-ser  itself  Is  being  perfected. 


This  JSEP  grant 
planar  Schottky  d  o- 
and  analyzed.  Consider 
to  have  NEP  higher  vV  : 
ten.  This  drawback  i 
integrabil ity  of  the  d 
carried  on  further  • 


'  'ploratory  research  with 
*  diodes  has  been  developed 

’t  ns,  these  diodes  are  expected 
:  »n  diodes  by  a  factor  of  about 
the  greater  ruggedness  and 
■Motion  1  development  is  now  being 


We  have  shown  t.ha*  sling  of  radiation  at  3b  GHz  into 
thin  film  metal-barr : er-me*  i  -t  :•»!'.  s  can  occur.  Sensitivities  within 
a  factor  of  about  three  *  o  J  ;  of  :>hottky  harrier  junctions  can  be 
obtained  for  direct  de*  >  f  »•  >t  ^  GHz.  The  noisf?  was  1/f  limited  up 
to  a  few  tens  of  megahen1  ?  whe-.  a  Mas  was  applied.  The  characteristics 
of  this  noise  spectrum  ar  j  ■  ’ tr i d stub! e  to  impurity  noise  and  can  possi¬ 
bly  be  alleviated  by  bet  Mr  : abr  icat i on  of  the  barrier  material.  In  the 
absence  of  the  1/f  spec  ire*  'Mgh  frequency  or  zero  bias),  ihe  noise  was 
Johnston  or  shot  in  c-tcr>. 


Heterodyne  det*  ;  \ *(  and  harmonic  mixing  using  a  72  GHz 
signal  were  also  serf  » -i“d .  G  iversion  losses  were  respectively  about 
ten  to  twenty  dh  greate'"  ^  reported  for  commercially  available 
detectors.  Hov;ever,  .  ’  : design  for  the  harmonic  mixing  experi¬ 
ment  was  not  used. 

These  experiments,  *;.«■  ’-e.  alts  of  which  are  described  in  detail 
elsewhere  [IS],  prov i  rfr,;,.-  evidence  that  the  metal-barr i er-metal 
structure  would  be  :•!;  *-r  i  :  to  G-nottky  barriers  at  higher  frequencies: 
frequencies  above  w:  i  do  Mo  spread ing  resistance  results  in  inferior 
performance  for  the  So  hot  -  gy  devices,  We  have  thus  initiated  a  study  of 
Me  fundamental  character ; of ;  s  of  motal-harr ier-metal  junction  struc¬ 
tures  at  shorter  wove  I  . 


An  initial  e x r  i  m  •  >. 
radiation  to  an  *‘d  t-* 

center  of  an  e 1  over  w.j  .  '  *  .• 

detection  was  nut.  h ; r ’•  . 
Hwang,  et  al  ,  [1?1  were  oy* 
structure  for  the  pin*  r  ; 
parable  to  a  p'>  i  nt,  '•  : 

corresponded  to  a  r:vqt-' 
detection  experiment  I  : 

;-;s*»i:.le  and  thot.  t  'usd  >s * 
r  n x per  i non t  s  <i: 

*  *  ^  ' '  n  t  1  n  i  J #  *  \  , 


if  demonstrated  the  coupling  of  118  pim 
i er-metal  structure  situated  at  the 
ib'rinn.  Although  the  sensitivity  of 
opl  i  ng  principles  first  suggested  by 
died  for  M-H— M  junctions.  The  edge 
:  al  an  area  of  ■  10“^r5  cm'0  w^jch  is  com- 
Thus  an  Pf  time  constant 
"■neb  shorter  than  118  pm.  This  initial 
Miat.  coupling  to  planar  junctions  is 
.vest  i  gabions  of  a  variety  of  nixing  and 
,,v*v'- 1  engths  are  possible.  This  work  is 
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(1)  "Characteristics  of  Integrated  MOM  Junctions  at  D.C.  and  at  Optical 
Frequencies 1 11  IEEE  Jour.  Quantum  Elec.,  QE-14,  159  (1978). 

(?)  "Electronic  Resonance  of  CARS,"  Phys .  Rev.,  18  (  1978). 

(3)  "Photo-Induced  Currents  in  Metal-Barrier-Metal  Junctions,"  Radia¬ 
tion  Energy  Conversion  in  Space ,  ed .  by  K,  W.  Billman,  1978.  ( Vol . 

61  of  Progress  in  Astronautics  and  Aeronautics  Series.) 

(4)  "Coupling  to  an  Edge  1  Metal-Ox ide-Metal 1  Junction  via  an  Evaporated 
Long  Antenna,"  Appl .  Phy3 .  Lett . ,  34,  823  (1979). 

(5)  "Small  Area  Metal-Oxide-Metal  Junctions  as  Picosecond  Photoemissive 
Detectors,"  in  Picosecond  Phenomena ,  Springer  Series  in  Chemical 
Physics,  vol.  4  (edited  by  C.  V.  Shank,  E.  P.  Ippen,  and  S.  L. 
Shapiro) . 

(6)  "Moving  Toward  NMM  Wave  Integrated  Circuits,"  Microwave  Journal , 
47-52  (June  1980) . 

(7)  "Progress  in  NMM  Wave  Integrated  Circuits,"  to  be  presented  at  the 
Fifth  International  Conference  on  Infrared  and  Millimeter  Waves, 
October  1980. 

The  following  doctoral  theses  were  completed  by  students  receiving 
support  from  JSEP: 

"Character istics  of  Metal-Oxide-Metal  Devices,"  Mordehai  Heiblum. 
The  abstract  is  as  follows: 

In  the  past  ten  years,  metal-oxide-metal  (MOM)  devices,  in  their 
point-contact  configuration,  have  been  used  for  harmonic-,  sum- 
and  difference-generation  at  frequencies  up  to  the  near  in¬ 
frared,  as  well  as  for  detection  of  visible  radiation. 

Because  of  the  lack  of  mechanical  stability  of  the  point-contact 
configuration,  planar  MOM  devices  are  desirable.  Earlier  ef¬ 
forts,  however,  have  resulted  in  devices  with  undesirably  large 
area.  In  this  work  I  report  a  new  integrated  device  with  area 
comparable  to  that  of  the  point-contact  configuration  (-10“ 
^cm^),  which  we  name  "Edge  MOM."  Fabrication  methods  and 
detection  characteristics  at  DC  and  at  10.6,  3.39,  and  0.6328  pm 
are  described.  A  detailed  description  of  the  Edge  MOM  junction 
and  its  behavior  under  the  influence  of  impinging  radiation  is 
given . 

An  extension  of  the  MOM  to  an  MOMOM  is  suggested  in  part  two  of 
the  work.  The  suggested  MOMOM  has  the  theoretical  capability  of 
amplifying  optical  signals. 


In  order  to  find  the  limitations  of  the  MOMOM,  hot-electron  dev¬ 
ices  such  as:  MOM-vaeuum-motal ,  MOM-semtoonduotor ,  and  MOMOM 
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are  theoretically  investigated.  As  a  result  a  proposed  struc¬ 
ture  with  effective  interaction  area  ( and  oxide 

o 

thickness  on  the  order  of  10A  is  proposed,  which  might  have  a 
transfer  ratio  of  0.3  -  0.5. 

Different  coupling  mechanisms  such  as:  wave-guiding,  antenna, 

and  sur face-plasmon  coupling  are  looked  into  and  frequency  limi¬ 
tations  are  estimated. 


”Metal-Barr ier-Metal  Junctions  for  Millimeter  Wave  Mixing  and 
Detection,”  Charles  William  Slayman.  This  thesis  has  the  following 
abstract : 

Metal-barrier-metal  (MBM)  junctions  were  used  to  mix  and  detect 
millimeter  wave  radiation.  The  MRM  junctions  used  in  this  work 
were  of  the  stable  thin  film  type  with  areas  of  roughly  1  jjm  x  1 
pm  and  fabricated  using  conventional  photol i thographic  tech¬ 
niques.  The  devices  tested  were  predominantly  of  the  Ni-NiO-Ni 
variety . 

A  brief  review  of  1)  the  theory  of  the  nonlinear  current  vs, 
voltage  (I-V)  characteristic  of  the  MBM,  2)  the  equivalent  cir¬ 
cuit  of  the  device,  and  3)  the  processing  steps  used  to  fabri¬ 
cate  the  Ni-NiO-Ni  junctions  is  presented.  Devices  of  different 
impedances  were  produced  by  varying  the  processing  procedure. 

The  I-V  character i st ic  ,  dynamic  resistance  Rcj yn  ,  and  nonlineari¬ 
ty  I”/If  were  measured  for  the  different  devices.  It  was  found 
that  the  nonlinearity  was  not  dependent  on  the  processing  param¬ 
eters.  The  low  temperature  (77K)  charac ter i st i cs  of  an  Al- 
AlpOo-Ni  and  a  Ni-NiO-Ni  were  measured.  Between  4 . 2K  and  77K, 
the  I-V  character i st ic  of  the  Ni-NiO-Ni  junction  did  not  change. 

A  brief  introduction  to  the  possible  noise  mechanisms  in  MBM 
junctions  is  presented.  Any  two  terminal  device  (both  of  whose 
electrodes  obey  Ferni-Dirac  statistics)  that  displays  shot  noise 
under  bias  will  display  shot  noise  at  zero  bias  equivalent  to 
the  thermal  noise  of  a  resistor.  This  fact  was  experimentally 
observed  in  a  9407,  Ni-NiO-Ni  junction.  Under  bias,  however, 
the  device  displayed  a  1/f^*1  noise  spectrum  up  to  100  kHz  and 
the  RMS  noise  current  was  found  to  be  proportional  to  the  bias 
current  and  not  the  square  root  of  the  bia3  current. 

Detection  experiments  were  performed  at  36  GHz  with  Ni-NiO-Ni 
junctions  ranging  from  180.:  to  17. 7  .  The  detected  signal  sa¬ 
turated  at  power  levels  above  10  pW.  The  millimeter  wave 
response  was  studied  under  short  pulse  (60  ns)  and  long  pulse  (6 
iis)  conditions  for  a  normal  and  a  burned  Ni-NiO-Ni  device.  The 
sensitivity  of  the  junction  was  found  to  increase  threefold  at 
77K.  A  comparison  of  the  performance  of  the  various  Ni-NiO-Ni 
junctions  to  point  contact  MBMs  and  commercial  Gohottky  detec¬ 
tors  is  made. 
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Fundamental  mixing  between  two  36  GHz  klystrons  and  harmonic 
mixing  between  a  36  GHz  klystron  and  a  73  GHz  klystron  was  per¬ 
formed  using  various  Ni-NiO-Ni  junctions*  Mixer  conversion  loss 
and  noise  temperature  were  studied  as  a  function  of  bias  and  lo¬ 
cal  oscillator  (LO)  drive.  Though  the  Ni-NiO-Ni  diodes  are  not 
as  sensitive  as  Schottky  mixers  and  detectors,  different  MBM 
systems  might  ultimately  prove  to  be  more  nonlinear. 

Significant  Accomplishments 

We  have  gained  an  improved  understanding  of  the  processes  involved 
in  the  fabrication,  by  the  two-step  etching  procedure,  of  silicon  mem¬ 
brane  Josephson  devices.  Successful  devices  have  been  made  and  tested. 

A  new  approach  to  tunable  quasi particl e  mixers  has  been  devised. 
The  novel  features  of  this  design  are  (a)  use  of  an  integral  planar 
antenna  to  improve  coupling  and  (h)  use  of  an  elec trical ly-tunable  vari¬ 
able  reactance  to  cancel  device  capacitance. 

The  usefulness  of  metal-barr ier-metal  junctions  as  millimeter-wave 
mixers  has  been  demonstrated  experimentally.  Results  indicate  that 
their  mixer  noise  temperatures  are  comparable  to  those  of  Schottky 
diodes,  although  conversion  is  less. 
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Scientific  Objective 


Our  objective  is  to  continue  our  research  efforts  on  large  scale 
and  nonlinear  circuits  which  were  initiated  several  years  ago.  The  many 
tools  and  results  that  we  have  developed  so  far  continue  to  serve  as  a 
foundation  for  our  future  research  and  will  he  used  to  attack  several 
interrelated  fundamental  problems  in  this  relatively  unexplored  but 
important  area. 

State  of  the  Art 

Network  theory  played  a  major  role  from  the  early  '10 *s  to  the  60*s 
in  the  development  of  the  present  day  communication  systems.  The  main 
thrust  of  network  theory  then  was  directed  at  the  analysis  and  synthesis 
of  linear  circuits  which  can  be  characterized  by  a  relatively  small 
number  of  differential  equations.  Although  in  recent  years  many  authors 
have  published  papers  on  nonlinear  circuits  and  large-scale  networks, 
the  extent  of  the  development  of  large-scale  nonlinear  circuits,  both  in 
theory  and  applications,  is  nowhere  near  the  kind  of  maturity  in  classi¬ 
cal  network  theory.  The  advances  during  the  past  decade  in  elecronic 
devices  and  integrated  circuits,  coupled  with  the  immense  impact  of  com¬ 
puters  on  circuit  analysis  and  design,  have  drastically  changed  the 
nature  of  the  problems  faced  by  circuit  theorists.  While  anlysis  of 
circuits  is  still  the  key  problem,  it  is  more  difficult  in  view  of  the 
size  and  complexity  of  the  circuit  and  the  inherent  nonlinearity  which 
exists  in  most  electronic  devices.  Furthermore,  the  scope  and  signifi¬ 
cance  of  "circuit  analysis"  has  been  much  broadened.  Theory  and  tech¬ 
niques  developed  from  large-scale  networks  and  nonlinear  circuit 
analysis  can  often  be  used  in  other  types  of  large-scale  network  and 
system  problems;  for  example,  computer  networks,  power  networks  and  com¬ 
plex  control  systems. 

It  should  be  noted  that  because  of  the  size  and  complexity  of  LSI 
and  VLSI,  the  objectives  and  approaches  in  analysis  are  considerably 
different  from  those  of  the  classical  network  theory.  First,  computa¬ 
tion  becomes  a  key  element  in  the  study  of  any  large  circuit  or  system. 
The  development  of  efficient  algorithms  is  crucial.  While  theoretical 
results  are  essential,  often  one  must  also  depend  on  heuristics. 
Second,  the  study  of  the  structure  of  a  large  circuit,  and  the  property 
of  interconnection  is  important.  Frequently,  large  problems  need  to  be 
divided  into  smaller  sub-problems  in  order  to  facilitate  analysis,  data 
storage  and  computation.  Third,  design  is  usually  accomplished  by 
experience  as  a  first  step  and  improved  on  by  repeated  analysis. 
Fourth,  such  problems  as  reliability  study,  tolerance  analysis,  diag¬ 
nosis  and  testing  are  important  and  should  he  kept  In  mind  in  the 
overall  design  process.  The  last  item  constitutes  a  new  area  of 
research  which  was  not  a  part  of  classical  network  theory.  All  these 
indicate  that  there  exist  n  large  variety  of  problems  of  a  eircuit- 
theoretlc  nature  which  need  to  he  studied* 

Progress  and  Publications  Since  Last  Major  Proposal 
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The  following  is  a  summary  of  research  progress  in  the  area  of  non¬ 
linear  and  large-scale  networks  which  were  published  since  the  last 
major  proposal . 

JK  Research  Progress  on  Piecewise -Linear 
Analysis  and  Piecewlso-L inear  Algebra 

Piecewi se-1 inear  techniques  have  been  used  extensively  in  circuits 
and  systems  theory  to  model  dc  non] inear  character i sties  of  electronic 
devices,  and  to  study  a  large  class  of  nonlinear  resistive  networks.  In 
1977,  Chua  and  Kang  introduced  new  analytical  representations  for  one¬ 
dimensional  piecewi se-1 i nrar  functions^.  This  piecewi  se-1  i  near 
representation  was  subsequently  generalized  to  the  multi-dimensional 
case:  it  is  applicable  for  any  m-dimensional  piecewi se-linear  function 
which  is  affine  over  convex  polyhedral  regions  bounded  by  linear  parti- 
tions^.  This  general i zat i on  results  in  a  tremendous  amount  of  data 
compression,  thereby  allowing  large-scale  piecewi se-1 inear  equations  to 
he  stored  in  significantly  less  computer  memory  space.  Application  of 
this  result  to  the  modeling  and  analysis  of  large-scale  networks  is 
presently  under  investigation. 

Although  motivated  originally  by  analysis  problems,  the  piecewi se- 
linear  representation  in  reference  2  led  naturally  to  important  applica¬ 
tions  in  nonlinear  circuit  synthesis .  In  particular,  systematic  methods 
for  synthesizing  nonlinear  networks  having  a  prescribed  scalar  or  mul¬ 
tidimensional  piecewi ue- 1  inear  function  have  been  developed^.  jn  the 
scalar  case,  precision  active-circuit  building  blocks  using  operational 
amplifiers,  transistors,  diodes,  and  resistors  have  been  derived  for 
realizing  piecewi so-* ’re  u-  driving  point  and  transfer  characteristic 
plots.  In  the  mul tid imet’s i onal  case,  methods  have  also  been  obtained 
for  realizing  a  nul  ti  t.ermi  nal  nonlinear  network  having  a  multidimen¬ 
sional  i  ecewi  se- 1  i  near  *.  rnnsfer  function.  Moreover,  these  methods  have 
been  generalized  by  synthesizing  nonlinear  n-ports  having  a  prescribed 
multidimensional  piecewi se- i i near  driving-point  function.  Although  most 
of  the  basic  buildinr  Mocks  are  grounded  active  networks,  converter 
networks  have  been  developed  for  transforming  such  grounded  networks 
into  floating  networks  coving  the  same  properties.  By  slight  modifica¬ 
tions  of  these  converter  networks,  other  useful  conversion  properties 
have  also  been  developed. 

More  recently,  this  multidimensional  piecewise-1 inenr  representa¬ 
tion  has  led  to  the  partial  resolution  of  a  heretofore  unsolved  conjec¬ 
ture  first  posed  in  In  a  recent  paper^,*  we  have  proved  that 
every  reciprocal  n-port  resistor  represented  by  our  n-dimensional 
piecewi se-1 inear  funct:  n  is  realizable  by  a  nonliner  circuit  containing 
only  /-terminal  piecewi inear  i  resistors  and  a  reciprocal  linear 
e  Jem  ,  such  as  '  p  +  q ') -p  or  i  transformer.  This  result  thereby  resolves 
one  of  the  many  fundament  u  questions  raised  in  reference  6.  In  [(artic¬ 
ular,  we  have  pre^d  the  ''-allowing: 

*The  second  author  (;>a;id  I.  Curtin)  was  supported  by  the  JSFP  Con¬ 
tract  F-^i4b?(}-?l-b0g7  # 
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Theorem  n-port  Reciprocity  Criteria 

A  piecewise-1 inear  n-port  described  by 


1,  -  a,  +  b,v,  +  ...  *  b,nv„  ♦  £  9U 
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Theorem  2 .  Reciprocal  n-port  Synthesis  Criteria 

A  piecewise-lmear  n-port  satisfying  the  above  reciprocity  criteria 
can  be  synthesized  using  only  passive  piecewi se-1 inear  7-terminal  resis¬ 
tors  and  a  (p+q)-port  transformer  if  and  only  if  the  set  of  algebra i c 
conditions  given  in  reference  5  satisfied. 

Both  theorems  1  and  ?.  could  not  be  proved  without  using  our  mul¬ 
tidimensional  piecewi se-1 inear  representation  whose  explicit  analytical 
represen tation  was  crucial  in  constructing  the  piecewi se-1 inear  alge¬ 
braic  proof. 

Research  Progress  on  Qua I i tat i ve  Nonlinear  Network  Analysis 

Among  the  many  important  unsolved  problems  posed  in  the  last  major 
proposal  concerning  the  f oundat i on  of  nonlinear  circuit  theory,  the  fol¬ 
lowing  are  some  of  the  recent  results  published  on  this  subject. 

a  *  Implicate ons  of  Capaci tor-only  Cutsets  and 

Inductor-only  Lo^ps  iji  Nonl  i  near  Networks.! 

Let  N  be  an  autonomous  dynamic  nonlinear  network.  Let  Np^  be  the 
associated  resistive  subnetwork  obtained  by  open  circuiting  all  capaci¬ 
tors  and  short  circuiting  all  inductors.  We  have  proved  the  following 
results  in  reference  7. 

(1)  Suppose  that  Np;-  has  only  isolated  operating  points.  Then  N  has 
only  isolated  equilibria,  if  and  only  if,  there  are  no  capacitor- 
only  cutsets  and  inductor-only  loops. 

(7)  If  the  above  condition  is  violated,  then  there  are  a  continuum  of 
equilibria  even  if  the  operating  points  are  isolated. 

(j)  Let  M  be  the  set  of  equilibria.  Then  each  trajectory  is  con¬ 
strained  to  lie  on  an  affine  submanifold  M* f  which  depends  on  the 

g 

initial  state,  3ueh  that  M  M  nas  only  isolated  points.  Hence, 
each  trajectory  behaves  as  if  it  has  only  isolated  equilibria.  The 
space  M  ,  because  of  its  nature,  can  be  considered  as  the  minimal 
dynamic  space  of  the  network. 

The  above  results  are  of  special  relevance  to  the  computer  simula¬ 
tion  of  nonlinear  electronic  circuits^.  The  presence  of  inductor-only 
loops  and  capacitor-only  cutsets  could  lead  to  serious  numerical  prob¬ 
lems  during  integration.  These  issues  have  been  the  suhject  of  several 
heated  debates^* ^  almost  70  years  ago.  The  results  presented  in  ?  pro¬ 
vide  a  definitive  and  rigorous  resolution  of  the  fundamental  issues 
raised  in  references  9  and  10. 

b.  Pynamics  of  Joseph  son- junction  Circuits.!! 

Although  seemingly  esoteric,  the  electromagnetic  properties  of 
vlosephson  junctions  have  been  used  in  applications  ranging  from  the 
measurement  of  mine-te  computers  with  pie<>  second  switching  times. 

A 1 t  hot  vh  nil*/  has  torn  published  concerning  Joseph son- junction  cir¬ 
cuits,  these  articles  are  either  quantum-mechanical  analyses  of  the 
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Josephson  effect  or  sundry  reports  of  experiments  which  reveal  some  very 
remarkable  phenomena  associated  with  circuits  containing  Josephson  Junc¬ 
tions,  Those  who  have  ventured  to  solve  the  differential  equations  of 
Josephson-junction  circuits  have  necessarily  restricted  themselves  to 
approximate  methods,  or  the  use  of  approximate  analog  models.  The 
situation  remains  somewhat  unsatisfactory  from  the  circuit  theorist's 
point  of  view,  as  only  a  very  limited  insight  is  available  into  the  gen¬ 
eral  circuit  behavior  of  the  devices. 

Our  research  on  Josephson-junction  circuits  was  conerned  mainly 
with  developing  a  unified  qualitative  theory  of  the  phenomena  associated 
with  Josephson-junction  circuits.  Some  recent  progress  has  been 
reported^.  In  particular,  the  features  of  the  d.c.  I-V  characteristic 
of  a  Josephson  junction  have  been  explained  rigorously  in  terms  of  the 
flows  on  a  cylindrical  phase  space.  The  junction  phase  difference  6 
attains  an  equilibrium  point  in  the  supercurrent  regime,  and  acts  like  a 
current-controlled  oscillator  in  the  finite  voltage  regime.  The  hys¬ 
teresis  in  these  characteristics  is  due  to  the  coexistence  of  an  equili¬ 
brium  point  and  a  periodic  solution. 

For  the  a.c.  excited  Josephson  junction,  the  dynamics  are  described 
for  a  simpler  circuit  model  as  trajectories  on  the  surface  of  a  torus. 
Using  the  concept  of  a  turning  point  to  define  the  various  possibilities 
of  periodic  flows  on  the  torus,,  the  character  of  harmonic  and  subhar¬ 
monic  oscillatory  waveform  has  been  derived.  Furthermore,  the  existence 
of  almost  periodic  waveforms  has  been  established.  The  presence  of 
constant-voltage  steps  is  related  to  the  structural  stability  of  the 
turning  point,  which  roughly  means  that  the  character  of  a  periodic 
oscillation  is  not  affected  by  smaM  enough  perturbations  in  the  excita¬ 
tion.  An  algorithm  is  given  in  reference  11  to  numerically  determine 
the  heights  of  these  constant  voltage  steps. 

A  precise  interpretation  is  also  given  in  reference  11  showing  how 
an  a.c.  excited  junction  is  subject  to  synchronization  phenomena,  and 
that  the  step  height  is  merely  the  entrainment  range,  or  locking  range, 
of  each  synchroni zation  event.  It  is  also  shown  how  synehroni zation ,  in 
general,  can  be  geometrically  interpreted  as  closed  trajectories  on  an 
integral  manifold  which  can  be  smoothly  transformed  into  an  n- 
dimensional  "torus 

Theory  <3  f  Symmetry  in  Non  linear  Circuits  12,13 

Simple  group-theoretic  concepts  have  been  used  successfully  in 
references  12  and  13  to  develop  a  rigorous  and  comprehensive  theory  of 
symmetry  for  nonlinear  elements  and  circuits.  This  theory  does  not  rely 
on  geometrical  arguments  or  other  ad  hoc  techniques  normally  invoked  in 
such  studies.  This  theory  also  unifies  all  forms  of  symmetry,  including 
rotation,  reflection  and  complementary  symmetry  into  a  3ingle  framework. 
It  also  includes  all  known  nonlinear  symmetry  principles  as  special 
cases.  Moreover,  a  general  method  for  identifying  all  symmetry  charac¬ 
teristics  possessed  by  a  nonlinear  multiterminal  element  or  circuit  has 
been  obtained.  Some  of  the  results  that  have  been  obtained  concerning 
symmetry  in  a  nonlinear  element  are: 
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(1)  Several  algorithms  for  synthesizing  a  nonlinear  multiport  or 
multiterminal  element  having  any  prescribed  form  of  symmetry  have  been 
derived.  In  particular,  various  examples  have  been  used  to  illustrate 
how  these  algorithms  can  he  used  to  derive  well-known  symmetrical  non¬ 
linear  circuit  modules  such  as  push-pull  amplifiers,  complementary- 
symmetric  amplifiers,  rectifiers,  modulators,  etc. 

(?)  A  reduction  algorithm  has  been  developed  which  allows  a  com¬ 
plicated  symmetric  element  to  be  analyzed  by  a  much  simpler  reduced  ele¬ 
ment  . 


(9)  A  general  principle  has  been  derived  for  applying  symmetry  to 
achieve  frequency  separation  in  nonlinear  communication  circuits  where 
the  even  harmonic  components  are  separated  from  the  odd  harmonic  eom- 
ponen ts . 
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Interaction  With  Other  Work  U n 1 1 s 

We  anticipate  interacting  with  the  CAD  efforts  of  Professors  Polak 
and  Sangiovanni-Vincentel 1 1 ,  as  we  1  i  as  the  large-scale  system  research 
of  Professors  C.  A.  Desoer ,  P.P.  Varaiyn,  R.  W.  Prodersen  et  al  .  (Work 
units  TSS-8 1 -2  and  5SD-R1-U  respectively^. 
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Scientific  Objective 

The  overall  long  term  objective  Is  to  design  control  structures  for 
unknown  systems  subject  to  random  disturbances.  The  class  of  structures 
being  investigated  is  the  so-called  "adaptive  feedback  control"  laws, 
the  unknown  system  is  any  countable  state  Markov  chain  whose  unknown 
transition  probabilities  are  parameterized .  Also  to  be  studied  is  the 
practically  important  case  where  the  "real"  system  does  not  belong  to 
the  class  of  a  priori  selected  models.  In  conjunction  with  this  work, 
deterministic  system  design  methods  and  evaluations  are  also  investi¬ 
gated  . 

State  of  the  Art 

Over  the  past  twenty-five  years  engineers  have  invented  a  variety 
of  procedures  for  designing  control  strategies.  The  sophistication  of 
these  strategies  has  grown  with  the  increase  in  the  use  of  computers 
which  enables  the  implementation  of  complicated  strategies  in  real  time 
(feedback  control).  But  all  these  procedures  rest  on  the  assumption 
that  the  information  available  about  the  system  and  the  calculation  of 
the  choice  of  control  which  is  based  on  this  information  is  centralized , 
that  is  it  takes  place  in  a  single  system. 

As  engineers  began  to  consider  the  control  of  large  systems*  (e.g., 
power  systems,  communication  networks)  in  a  scientific  manner,  it  became 
clear  that  "classical"  procedures  of  design  are  inapplicable  to  the 
problem  of  Combinatorial  Complexity.  As  the  number  of  variables 
describing  the  system  increases,  the  amount  of  computation  necessary  to 
apply  classical  control  strategies  grows  even  faster,  soon  exhausting 
the  capability  of  existing  and  foreseeable  computing  capability.  It  was 
soon  recognized  that  a  way  of  dealing  with  this  complexity  is  to  insist 
on  decentralized  control  strategies ^ .  After  this  recognition, 
research  has  proceeded  in  three  directions.  First ,  there  has  been  the 
attempt  to  derive  design  principles  applicable  to  a  fairly  general  class 
of  decentralized  controllers^. 6.  But  this  attempt  was  soon  frustrated 
by  a  new  difficulty,  Informational  Complexity** .  A3  soon  as  information 
and  control  is  decentralized,  the  information  available  and  the  effec¬ 
tiveness  of  the  decision  taken  by  each  controller  depend  upon  the  deci¬ 
sions  taken  earlier  by  other  controllers.  This  "simultaneity"  repro¬ 
duces  the  Combinatorial  Complexity  at  another  level  and  renders  inappli¬ 
cable  existing  techniques  of  controller  design.  The  difficulties 
created  by  Informational  Complexity  are  not  well  understood.  On  the  one 


*Systematic  attention  to  the  control  of  large  systems  began  less  than 
ten  years  ago.  An  indication  of  this  is  given  by  the  fact  that  the  IEEE 
Control  Systems  Society  established  its  Technical  Committee  on  Large 
Scale  Systems  in  1969*  its  first  editorial  on  large  systems  appeared  in 
197^,  and  its  first  Special  Issue  on  Large  Systems  In  April  1978. 

**We  are  indebted  to  Professor  Y.  C.  Ho  (personal  communication)  for 
his  conceptualization  of  these  two  kinds  of  complexity. 
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hand  we  have  several  counterexamples  of  decentralized  control  problems 
which  appear  quite  innocent,  but  whose  optimal  strategies  are  extremely 
complex? *8i 16,  on  the  other  hand,  we  now  know  a  (small)  class  of  decen¬ 
tralized  information  structures  for  which  optimal  control  strategies  do 
have  a  simple  structure11. 

The  second  direction  of  research  has  led  to  the  search  of  optimum 
decentralized  controllers  for  what  have  been  called  quasi-static  control 
problems^.  For  these  problems  optimum  controllers  can  have  a  simpler 
structure  since  only  "steady  state”  behavior  is  of  interest.  Some  very 
interesting  solutions  have  already  been  proposed  for  the  control  of  net¬ 
work  flows™  and  for  routing  messages  in  a  communication  network11. 

Whereas  the  second  direction  described  above  seeks  to  reduce  com¬ 
plexity  by  simplifying  the  control  task  (attention  to  steady  state  and 
ignoring  transient  behavior),  the  third  direction  reduces  complexity  by 
restricting  a  priori  the  structure  of  the  class  of  decentralized  con¬ 
trollers,  In  designing  fixed  structure  controllers  which  guarantee  ade¬ 
quate  response  against  random  disturbances,  the  attempts  to  date  have 
been  ad  hoc  it  seems  clear  that  considerable  advance  can  be  made 
if  fixed  structure  decentralized  controllers  are  investigated  more  sys¬ 
tematically.  This  prediction  is  bolstered  by  the  results  recently 
obtained  in  the  study  of  decentralized  linear  feedback  controllers  which 
guarantee  adequate  response  against  "deterministic  but  unknown"  distur¬ 
bances  . 

With  the  goal  of  designing  robust  controllers  for  large  scale  sys¬ 
tems,  the  following  approache3  are  relevant.  First  there  is  a  very  con¬ 
siderable  literature  on  the  stability  of  interconnected  systems  based  on 
the  Lyapunov  method.  In  terms  of  generality  of  the  results  and  of 
technical  sophistication,  the  paper  of  Rasmussen  and  Michel  deserves 
special  mention1?.  Second,  the  problem  can  be  approached  from  the 
input-output  point  of  view:  given  the  input-output  properties  of  the 
subsystem  and  the  topology  of  the  1 nterconnection ,  describe  the  I/O  pro¬ 
perties  of  the  overall  system.  Under  the  joint  sponsorship  of  NSF  and 
JSEP  we  conducted  a  thorough  investigation  of  this  problem'^.  We  refer 
the  reader  to  its  46  item  reference  list.  Recently  we  produced  an 
improved  formulation  of  the  problem  which  led  us  to  nece33ary  and  suffi¬ 
cient  conditions  for  stability,  even  for  the  nonlinear  case3?.  The  goal 
of  this  research  is  to  study  the  input-output  properties  of  the  overall 
system.  The  advantage  of  the  approach,  namely  graphical  decomposition, 
is  that  it  applies  equally  to  nonlinear  time-varying  just  as  well  as 
linear  time-invariant  systems,  either  lumped  or  distributed.  The  graph¬ 
ical  decomposition  technique  operates  as  follows:  the  digraph  of  the 
interconnected  subsystem  is  decomposed  Into  an  acyclic  interconnection 
of  strongly  connected  components  and  then  a  minimal  essential  set  is 
extracted  to  further  simplify  the  analysis  of  the  strongly  connected 
components. 

Oi  the  decentra] 1  zed  control  of  large  systems,  especially  as  it 
relates  to  tracking  and  disturbance  rejection,  the  work  of  Davison 19*?0 
deserves  special  mention.  We  used  it  as  well  as  graphical  decomposi t ion 
techniques  to  Investigate  the  eigenvalue  assignment  and  stabilization 
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problem  of  linear  time-invariant  interconnected  systems21;  one  main 
result  Is  that  we  can  achieve  eigenvalue  assignment  provided  we  operate 
«  at  the  column  subsystem  level,  whereas  if  wc  restrict  ourselves  to  sum¬ 

mation  node  level,  examples  show  that  eigenvalue  assignment  cannot  be 
achieved  In  general,  whereas  we  prove  that  we  can  achieve  stabiliza¬ 
tion21. 

On  the  more  specific  topic  of  tracking  and  disturbance  rejection, 
the  considerable  literature  and  several  approaches  used  are  reviewed  in 
reference  22.  Furthermore,  the  principal  results  are  given  a  stream¬ 
lined  self-contained  derivation  in  reference  23.  Unfortunately,  all 
this  literature  considers  only  linear,  usually  lumped  systems  and 
discusses  robustness  of  the  servocompensators  only  within  that  class22*. 
The  effect  of  nonlinear  plant  perturbation  on  optimal  control  design  has 
recently  been  examined  in  references  25,26;  in  the  broader  aspect  of  I/O 
stability  in  references  27-29* 

Onlv  one  recent  report  treats  the  tracking  problem  for  nonlinear 
systems 3*. 

Adaptive  control  strategies  are  increasingly  being  suggested  as  an 
appropriate  solution  to  the  situation  where  the  system  to  be  controlled 
is  not  known  and  is  subject  to  random  disturbances.  Adaptive  control  is 
a  compromise  between  two  extreme  alternatives.  The  first:  assume  the 
unknown  system  is  equal  to  a  ’'nominal”  system,  may  give  a  performance 
which  is  unacceptably  degraded;  the  second:  formulate  the  problem  as  a 
Bayesian  decision  problem  and  solve  it  by  Dynamic  Programming,  while 
formally  optimal,  it  cannot  be  "solved”  except  in  the  most  trivial 
instances. 

In  adaptive  control,  one  identifies  the  unknown  system  using  previ¬ 
ous  data;  and  then  uses  the  control  assuming  the  current  estimate  of  the 
unknown  system  to  be  the  true  system.  As  time  progresses,  and  more  data 
is  collected,  the  uncertainty  regarding  the  system  declines  and  the  per¬ 
formance  improves.  The  crucial  question  is  how  fast  this  improvement  in 
the  knowledge  and  performance  of  the  system  occurs.  Previous  research 
has  been  primarily  focused  on  linear  systems. 
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Progress  and  Relevant  Publications  Since  Last  Major  Proposal 


_.  Under  Combined  Sponsorship 

1 .  Control  of  Markov  Chain  Models 

The  research  reported  here  is  jointly  sponsred  by  NSF  under  Grant 
ENG79-03879.  We  have  extended  initial  work35  to  propose  a  multi-layer 
structure  for  controlling  large  Markov  chalns36.  The  name  ideas  can  he 
used  to  construct  a  finite  clain  model  for  n  one-dimensional  continuous 
diffusion^*!.  Since  the  finite  model  Is  obtained  from  the  continuous 
one  by  optional  sampling,  therefore  at  these  sampling  times  the  model  In 
exact .  This  feature  makes  the  model  more  attractive  than  the  more  usual 
grid  approximations.  However,  our  approach  suffers  from  a  major  defect. 
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not  present  in  other  methods.  In  that  it  is  limited  at  present  to  scalar 
processes.  It  is  an  important  difficult  task  of  future  research  to 
extend  the  results  in  to  the  multidimensional  case. 

Adaptive  control  has  become  a  major  direction  of  research  in  con¬ 
trol  engineering.  (One  measure  of  its  importance  is  the  fact  that  six 
sessions  of  the  1980  JACC  were  devoted  to  it.)  The  bulk  of  the  effort 
in  adaptive  control  continues  to  be  devoted  to  linear  system  models. 
However,  our  studies  are  directed  to  general  Markov  chain  models.  In 
our  first  set  of  results^ »c J  we  have  essentially  completely  analyzed 
the  case  where  it  is  assumed  that  the  unknown  system  parameter  belongs 
to  an  a  priori  specified  finite  set.  In  this  case  we  have  shown  that 
(i)  the  parameter  estimate  converges  in  a  finite  random  time  and  (ii) 
the  estimate  correctly  identifies  the  closed  loop  transition  probabili¬ 
ties  but  may  not  be  able  to  identify  the  open  loop  transition  probabili¬ 
ties.  This  result  inspired  Doshi  and  Shreve^  to  propose  a  randomiza¬ 
tion  of  the  control  policy,  since  such  "noisy”  input  achieves  open  loop 
identifiability . 

In  we  have  attempted  to  study  the  case  where  the  state  space 
of  the  chain  is  countable  and  where  the  unknown  system  parameter  is  in  a 
compact  set.  The  change  from  finite  to  the  compact  case  creates  diffi¬ 
culties.  It  is  no  longer  possible  to  prove  that  the  parameter  estimate 
converges.  Indeed,  recently  Kumar has  given  an  exaropl  *  where,  with 
probability  one,  the  estimate  does  not  converge.  However,  the  sequences 
of  estimates  always  have  certain  (so-called  "frequent")  limit  points 
where  the  closed-loop  transition  probabilities  are  correctly  identified. 
Finally,  two  versions  of  randomized  control  policies  are  shown  to 
achieve  open  loop  identifiability. 

2.  Design  of  Deterministic  Control  Systems 

Fundamental  studies  in  nonlinear  feedback  systems  have  been  com¬ 
pleted:  a  general  theory  of  multi -input  multi-output  nonlinear  feedback 
systems  has  been  completed  and  publ ished^ .  It  gives  a  complete  gen¬ 
eralization  of  the  theory  that  H.  S.  Black  presented  for  the  linear 
time-invariant  single-input  single-output  feedback  system.  In  addition, 
it  gives  a  completely  general  demonstration  that  the  linearized  Inverse 
return  difference  is  the  quantity  that  controls  the  degree  of  desensiti¬ 
zation  to  plant  and/or  output  disturbances.  It  is  also  shown  to  play  a 
key  role  In  linearizing  the  I/O  map  of  the  closed-loop  system. 

The  streamlining  of  the  theory  of  the  generalized  Nyqulst  criterion 
and  Its  generalization  to  the  distributed  case  has  been  pub l Inherit  11 . 
We  hear  that  the  method  has  been  programmed  and  is  being  used  at  TRW  In 
Los  Angeles.  The  features  of  our  generalized  Nyquist  criterion  is  that 
(a)  it  applies  to  the  distributed  case,  and  (b)  it  does  not  require  any 
reference  to  Riemann  surfaces,  branch  cuts,  etc.  that  previous  deriva¬ 
tions  required. 

An  Interesting  approach  to  the  study  of  linear  feedback  systems  due 
to  Zames  has  been  adapted  to  the  problem  of  designing  multivariable 
feedback  systems  with  stable  plants.  This  design  method  is 
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revolutionary  in  that  (a)  it  is  based  on  the  global  parameterization  of 
the  four  transfer  functions  pertaining  to  the  problem,  and  (b)  the  sta¬ 
bility  of  the  design  is  easily  checked  without  having  recourse  to  such 
things  as  Nyquist  diagrams  etc.  The  design  algorithm  allows  complete 
freedom  -  modulus,  of  course,  the  -  zeros  of  the  plant  -  in  the 
choice  of  I/O  map:  therefore  we  can  diagonalize  the  I/O  map,  assign 
zeros  independently  to  each  channel,  assign  poles  Independently  to  each 
channel  and  achieve  tracking  requirements.  The  simplest  version  of  the 
algorithm  was  presented  at  the  1980  JACcUl,  an  extension  of  it  will  be 
presented  at  the  1980  CDC^.  The  first  paper  describing  the  algorithm 
has  been  accepted  by  the  IEEE  Trans,  on  Automatic  Control.  We  are  work¬ 
ing  out  several  extensions:  one  of  our  goals  is  to  obtain  a  completely 
algebraic  theory  of  design^1!.. 

Many  studies  devoted  to  the  effect  of  the  plant  perturbation  on  the 
stability  of  the  system  require  that  the  perturbation  be  represented  by 
stable  operators.  This  is  toally  unrealistic  in  the  case  of  an  unstable 
plant:  indeed,  any  perturbation  in  an  unstable  pole  requires  an 

unstable  multiplicative  perturbationk*l . 

B.  Under  NSF  Exclusive  Sponsorship 

In  NSF  grant  ENG78-09032  Desoer  is  investigating  control  problems 
for  discrete-timetrol  and  continuous  timefnl  systems.  Also,  into], 
Desoer  et  al.  demonstrated  that  the  design  algorithm  developed  under 
JSEP  sponsorship  did  not  impose  unnecessary  constraints  on  the  choice  of 
the  I/O  map.  With  a  master's  student,  we  are  tying  the  design  algorithm 
technique  with  Professor  Polak's  inequality  constraint  approach  to 
design.  Simple  examples  show  that  the  method  works  and  gives  a  rational 
technique  for  selecting  parameters^ . 

C.  Significant  Accomplishments 

The  spread  of  computer -assisted  control  has  increased  the  possibil¬ 
ity  of  Implementing  adaptive  controllers  which  "tune"  the  control  inputs 
to  automatically  detected  changes  in  operating  conditions.  Considerably 
more  work  needs  to  be  done,  and  this  will  take  several  years,  before 
adaptive  controllers  can  be  implemented  on  a  routine  basi3. 

Work  on  statistical  identification  has  shown  us  that  it  is  possible 
to  design  experiments  to  identify  an  unknown  system,  and  engineers  have 
developed  many  methods  for  controlling  a  known  system.  The  aim  of  adap¬ 
tive  control  is  to  carry  out  both  tasks  -  identification  and  control  - 
concurrently.  Our  investigations,  and  the  work  of  others,  have  proved 
that,  in  principle  and  under  precisely  stated  conditions,  this  is  indeed 
possible.  Others  have  shown  where  adaptive  controllers  will  operate 
satisfactorily  for  linear  systems  of  the  form: 

z  aiyt-i  s  *  Vt-i  ♦  2  ciwt-i 

where  yt,  ut,  w^  are  the  output,  input  and  white  noise  and  the  ,  b* , 
are  unknown  parameters.  Our  work  has  resulted  In  conditions  for  the 
satisfactory  control  of  systems  of  the  form: 
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xfi  *  f(xt.wt,ut,theta) 

where  xt  is  a  discrete  state,  and  theta  is  the  unknown  parameter. 

A  very  exciting  development  in  control  system  design  occurred:  by 
adapting  a  global  parameterization  of  Zames,  we  have  developed  an  algo¬ 
rithmic  method  of  design  that  guarantees  system  stability,  that  quaran- 
tees  a  strictly  proper  controller  and  that  can  accomplish  the  following; 

(a)  decouple  the  system,  (b)  assign  the  poles  individually  in  each  chan¬ 
nel,  (c)  except  for  the  C^-zero  imposed  by  the  plant,  assign  the  zeros 
in  each  channel,  and  (d)  achieve  tracking  requirements.  It  is  apparent 
that  the  technique  is  very  general  and  applies  to  any  algebraic  system 
of  a  certain  kind  (e.g.,  continuous-time  or  discrete-time,  lumped  or 
distributed).  The  method  lends  itself  naturally  to  interactive 
computer-aided  design  methodologies. 
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Scientific  Objective 

This  project  is  a  continuation  of  Project  ISS-81-3.  Computer  Optim¬ 
ization  of  Electronic  Circuits.  Our  chief  objective  remains  the 
development  of  an  interactive,  optimization-based,  computer-aided  design 
methodology  for  the  design  of  electronic  circuits.  To  this  end,  we  pro¬ 
pose  to  continue  developing  specialized,  dovetailing,  single  and  multi¬ 
objective  optimization  algorithms,  circuit  simulation  techniques  and 
device  models.  In  addition,  since  the  development  of  progressively  more 
complex  systems  makes  sense  only  if  these  systems  can  be  maintained,  we 
propose  to  perform  research  on  system  diagnosabil i ty  and  fault  analysis. 

The  optimization  algorithms  to  be  developed  must  be  capable  of 
solving  complex,  nondi fferentiable ,  semi-infinite  optimization  and 
trade-off  problems  which  arise  in  electronic  circuit  design  when  toler¬ 
ances,  trimming,  and  distributed  constraints  are  taken  into  account. 
Our  simulation  procedures  must  be  capable  of  producing  not  only  circuit 
responses,  but  also  the  derivatives  of  these  responses  with  respect  to 
design  parameters.  In  addition,  they  must  have  special  variable  preci¬ 
sion  features  which  result  in  a  reduction  of  computing  time  at  low  pre¬ 
cision,  so  that  an  optimization  algorithm  can  call  such  a  procedure  with 
low  precision  when  far  from  a  solution  and  increase  precision  adap¬ 
tively,  as  a  solution  is  approached.  The  effect  of  such  variable  preci¬ 
sion  computation  is  to  considerably  reduce  computing  times.  Since  pro¬ 
duction  tolerances  can  be  imposed  on  only  a  very  small  number  of  device 
parameters,  and  since  the  optimization  algorithms  become  extremely  inef¬ 
ficient  when  tolerances  on  a  large  number  of  parameters  are  specified, 
device  models  must  be  developed  in  terms  of  production  parameters  by 
making  use  of  the  fact  that  the  variation  of  many  circuit  parameters  can 
be  expressed  in  terms  of  variations  of  production  parameters. 

Finally,  our  theoretical  work  will  be  implemented  and  experimen¬ 
tally  evaluated  in  an  interacti ve-computer  aided  design  facility  which 
we  are  developing. 

State  of  the  Art 

Optimization-based  computer-aided  design  is  certainly  not  a  new 
phenomenon,  see  for  example  [B9, B1 1 ■ B13» B15, B20, Cl ,K 1 ,K2, P7* H16] .  How¬ 
ever,  until  quite  recently,  its  use  was  severely  limited  by  the  fact 
that  the  available  optimization  algorithms  were  not  able  to  solve  design 
problems  with  tolerances,  trimming  and  distributed  constraints,  such  as 
those  imposed  on  frequency  or  transient  responses  [B20].  The  reason  for 
this  is  that  such  a  design  problem  frequently  assumes  the  form 

min{f(x) igi(x)  <  0,  i  =  1.2,. ..I;  max  hi(x,w)  <  0,  j  =1,2. ...J; 

~  wOJ 

max  max  min  max  zh  x  ,e ,  s ,  t , )  <  0}  (1) 
eCE  s£S  t^T  kCK 

where  r^(.)  represents  simple  constraints  on  the  nominal  design  vector 
x,  hJ(.,.)  represents  nominal  frequency  (or  time  )  dependent  constraints 
on  impedances  or  gains,  etc.,  and  z  represents  the  requirement  that  per¬ 
formance  be  maintained  in  the  presence  of  production  errors  e  which  can 
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Progress  and  Publications  Since  La3t  Major  Proposal 
A.  Under  JSEP 

1.  Professor  A.  Sangiovanni  Vincentelli  and  E.  Lelarasme,  a  graduate 
student,  have  developed  a  new  inner  approximations  method  for  the  design 
of  electronic  circuits.  Their  results  were  presented  in  the  invited 
paper  "An  Inner  Approximations  Algorithm  for  the  Design  of  Electronic 
Circuits, "  at  the  1979  International  Symposium  on  Mathematical  Theory  of 
Networks  and  Systems,  Delft,  Holland,  July  3-6.  1979. 

2.  Professor  A.  Sangiovanni  Vincentelli  has  developed,  with  N.  B.  Rab- 
bat  and  H.  Y.  Hsieh  of  IBM,  a  new  algorithm  for  the  analysis  of  larg^ 
scale  nonlinear  circuits  in  the  time  domain.  Their  results  were 
reported  in  the  paper  "A  Multi-level  Newton  Algorithm  for  the  Analysis 
of  Large  Scale  Nonlinear  Circuits  in  the  Time  Domain,"  IEEE  Trans .  on 
Circuits  and  Systems,  Special  Issue  on  Computational  Methods,  September 
1979. 


3.  Professor  A.  Sangiovanni  Vincentelli  has  started  a  research  effort 
in  computer  aided-design  of  power  electronic  circuits.  The  first 
results  of  this  effort,  co-authored  with  M.  Venturini  and  P.  Wood, 
"Models  for  Tv.v*' i stors  and  Diodes  in  Digital  simul ation , "  was  presented 
at  the  Power  El  eel.  Specialists  Conference,  San  Diego,  June  18-?3.  1979. 


4.  Professor  A.  Sangiovanni  Vincentelli  and  V.  Visvanathan,  a  graduate 
student,  have  developed  a  theory  for  the  fault  diagnosis  of  nonlinear 
memoryless  systems  and  a  technique  for  the  design  of  diagnosable  sys¬ 
tems.  Their  work  was  reported  in  "DC  Testing  of  Parametric  Faults  in 
Nonlinear  Systems,"  IEEE  International  Symposium  on  Circuits  and  Sys¬ 
tems,  Houston,  Texas,  April  <?8-30,  1980. 


B,  Professor  A-  Sangiovani  Vincentelli  and  E.  Lelarasme,  a  graduate 
student,  have  submitted  the  paper  "Time  Domain  Sensitivity  Analysis  by 
the  Perturbation  Method,"  to  IEEE  Trans .  on  Circuits  and  Systems . 

6.  Professors  E.  Polak,  A.  Sangiovanni  Vincentelli,  et  al ,  have 
described  our  interactive  optimization  software  system  in  "A  Software 
for  Optimization-Based  Computer-Aided  Design,"  TEEF  Internal ional  Sympo¬ 
sium  on  Circuits  and  Systems,  Houston,  Texas,  April  ?8-30,  1080. 


.  Professors  F.  Po 
have  made  a  consider 
design  centering,  to 
[■resented  by  Polak 
of  the  Optimal  Desig 
Inti  .  Symp.  on  Ma 
of  Technology,  Holla 
tional  Aspects  of 
Problem,"  IEEE  Trans 
A  considerable  imp 
Implementabie  Algori 
and  Tuni  ig  FroM  mm 


lak,  A.  Sangiovanni  Vincentelli  and  D.  Q.  Mayno 
able  amount  of  progress  in  developing  algorithms  for 
lernneing  and  tuning.  An  important  breakthrough  was 
and  Sangiovanni  Vincentelli  in  "Theoretical  Aspects 
n  Centering,  Tolerancing  and  Tuning  Problems,"  *4th 
th.  Theory  of  Networks  and  Systems,  Delft  University 
nd ,  July  3-6,  1<)79,  and  in  "Theoret  i  ml  and  Computa¬ 
ble  Optimal  Design  Centering,  Tolerancing  and  Tuning 
on  Circuits  amd  Systems ,  Vol .  0AS-P6 ,  No.  9,  1079. 
rovement  on  this  method  is  developed  by  Polak  in  "An 
thm  for  the  Optimal  Design  Centering,  Tolerancing 
presented  at  the  IRTA  International  Symposium  on 
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be  counteracted  by  adjusting  the  trimming  parameters  t.  The  variable  s 
may  represent  time  or  frequency  or  temperature,  etc.  The  above  problem 
is  easily  recognized  as  a  non-differentiable,  semi-infinite  optimization 
problem,  i.e.,  about  as  difficult  an  optimization  problem  as  one  cares 
to  think  of. 

Although  there  is  a  fairly  large  literature  on  algorithms  for  vari¬ 
ous  classes  of  semi-infinite  and  nondifferentiable  optimization  prob¬ 
lems,  see  e.g.,  [B25, B26.B27, D8, D9,G6,G7,G8, H12.M10, P14, Pt6, P17a ,H1 ,T1 ] , 
we  find  that  many  of  these  algorithms  are  either  conceptual  (i.e., 
involve  operations  that  one  is  not  able  to  perform)  or  address  them¬ 
selves  to  very  restrictive  classes  of  problems,  such  as  those  with  one- 
dimensionally  convex  constraint  sets  [ R9. B1 1 , B1 3.B15 ] .  As  far  as  gen¬ 
eral  purpose,  implementable  algorithms  are  concerned,  there  appear  to  be 
only  three  [P1M,G6,T1]  which  solve  the  simpler  case  of  (1)  when  there 
are  no  tolerance  and  trimming  constraints  present,  two  [G6,P4]  which 
allow  for  tolerances  and  only  one  [P5]  which  solves  (1)  in  its  complete 
generality.  All  of  these  three  algorithms  were  co-authored  by  E.  Polak. 

In  addit^ii  to  the  difficulties  mentioned  above,  the  use  of  optimi¬ 
zation  techniques  in  electronic  circuit  design  was,  and  still  is, 
severely  handicapped  by  the  fact  that  existing  simulation  codes,  e.g. 
the  ones  described  in  [K1,J6#Z2],  which  must  be  called  for  function  and 
derivative  evaluation,  have  inadequate  provisions  for  derivative  calcu¬ 
lations.  For  example,  if  we  consider  a  sample  of  the  best  known  com¬ 
puter  aided  analysis  programs,  such  as  ASTAP  [A6],  CIRCUS  [D7],  ECAP 
[E2] ,  NET  2  [M7 ] ,  SCEPTRE  [B61],  SUPERSCEPTRE  [S12],  SPICE  2  [M3]  and 
BELAC  [G12],  we  find,  quite  typically,  that  ASTAP  computes  no  sensitivi¬ 
ties,  SPICE  2  computes  only  small-signal  sensitivities,  and  none  of  them 
compute  temperature  sensitivity.  We  are  now  beginning  to  develop  an 
effort  to  extend  these  programs  so  that  required  derivatives  will  also 
be  computed.  In  addition  to  the  simulation  programs  mentioned  above, 
programs  based  on  tearing  algorithms  [B28,C9,H5,H6,K8,M1M,R1,S3- 
S8,S10,S1 1 ,W5]f  and  on  model  approximation  and  decoupling  [C6.F1]  are 
also  being  considered  for  extension  since  they  are  more  efficient  in 
large-system  applications. 

As  specialized  optimization  algorithms  began  to  be  developed  for 
electronic  circuit  design,  such  as  [Pi4,P5],  it  became  obvious  that 
existing  device  models  were  completely  incompatible  with  the  need  of 
keeping  the  total  number  of  toleranced  variables  small.  Quite  recently 
R.  G.  Meyer  was  able  to  show  that  more  suitable  models  can  be  developed, 
because  many  of  the  important  variables  are  correlated,  and  that  their 
behavior  can  be  expressed  in  terms  of  a  reasonably  small  number  of  vari¬ 
ables  the  tolerances  on  which  can  be  controlled  in  the  production  pro¬ 
cess,  such  as  sheet  thickness,  resistivity,  device  scaling  factors,  epi¬ 
taxial  doping,  etc. 

As  far  as  fault  diagnosis  is  concerned,  we  find  theoretical  results 
for  the  diagnosability  of  linear  circuits  and  .systems 
[ R1 , B3*B4, E3* LI , N 1 ,  R2, R3*S1  ,S2,S1 1  ,T1  ,T*I-T6] .  Some  preliminary  results 
for  the  diagnosability  of  memoryless  nonlinear  circuits  and  systems  were 
co-authored  by  A.  Sangiovanni  Vincentelli  [SI?].  However,  there  appear 
to  be  no  diagnosability  results  for  dynamic  nonlinear  circuits  and  sys¬ 
tems  . 
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Numerical  Methods  in  Engineering,  Versailles,  France,  December  10-15, 
1979,  and  also  to  appear  in  JOTA.  Mayne,  Polak  and  Voreadis  will 
present  the  paper  "A  Cut  Map  Algorithm  for  Design  Problems  with  Parame¬ 
ter  Tolerances"  at  the  19th  IEEE  Conference  on  Decision  and  Control, 
Albuquerque,  N.M.,  December  10-12,  1980.  This  is  a  highly  efficient 
algorithm  for  use  when  only  feasibility  and  not  optimality  need  to  be 
obtained.  Finally,  Mayne,  Polak  and  Sangiovanni-Vincentell i  have  sum¬ 
marized  their  work  in  the  survey  paper  "Computer  Aided  Design  via  Optim¬ 
ization"  which  was  presented  at  the  IFAC  Workshop  on  Control  Applica¬ 
tions  of  Nonlinear  Programming,  Denver,  Co.,  August  1979.  This  This 
paper  has  also  been  provisionally  accepted  by  Automatics. 

8.  Professor  E.  Polak  and  A.  N.  Payne,  a  graduate  student,  have 
developed  an  optimal  process  for  multi-objective  decision  making.  Their 
work  was  reported  in  the  paper  "An  Interactive  Rectangle  Elimination 
Method  for  Bi-objective  Decision  Making",  IEEE  Trans,  on  Control,  vol . 
AC -2 5,  No.  3,  1980. 

9.  Professor  E.  Polak  presented  an  invited  paper  "On  the  Nature  of 
Optimization  Problems  in  Engineering  Design"  at  the  10th  International 
Symposium  on  Mathematical  Programming,  Montreal,  Canada,  August  27-31, 
1079. 


10.  Professors  E.  Polak  and  D.  0.  Mayne  have  developed  two  super- 
linearly  convergent  algorithms  for  constrained  optimization  problems. 
The  first  is  of  the  secant  type  and  is  described  in  "A  Robust  Secant 
Method  for  Optimization  Problems  with  Inequality  Constraints,"  to  appear 
in  JOTA.  The  second  one  is  of  the  recursive  quadratic  programming  type, 
with  an  automatic  feature  for  setting  the  exact  penalty  function,  is 
described  in  "A  Superl inear ly  Convergent  Algorithm  for  Constrained 
Optimization  Problems"  at  the  10th  International  Symposium  on  Mathemati¬ 
cal  Programming,  Montreal,  Canada,  August  27-31*  1979.  This  paper  is 
also  under  consideration  for  inclusion  in  a  Mathematical  Programming 
Studies  volume. 

11.  Professors  E.  Polak  and  D.  Q.  Mayne  have  developed  several  methods 
for  the  solution  of  nonlinear  inequalities  in  a  finite  number  of  itera¬ 
tions.  Their  work  is  reported  in  "On  the  Finite  Solution  of  Nonlinear 
Inequalities",  IEEE  Trans .  on  Automatic  Control ,  vol.  AC-24,  No.  3, 
1979,  with  A.  J.  Heunis,  in  "Solving  Nonlinear  Inequalities  in  a  Finite 
Number  of  Iterations",  to  appear  in  JOTA ;  and  in  "Algorithms  for  the 
Design  of  Control  Systems  by  the  Method  of  Inequalities"  and  in  "A  Qua- 
dratically  Convergent  Algorithm  for  Solving  Infinite  Dimensional  Ine¬ 
qualities",  both  presented  at  the  18th  IEEE  Conference  on  Decision  and 
Control,  Fort  Lauderdale,  Florida,  December  12-14,  1979. 

12.  Professor  E.  Polak  and  Dr.  C.  Gonzaga  have  developed  a  special 
outer  approximations  theory  for  computer-aided  engineering  design. 
Their  work  is  reported  in  "On  Constraint  Dropping  Schemes  and  Optimality 
Functions  for  a  Class  of  Outer  Approx imat ions  Algorithms,"  SIAM  ,J.  Con- 
trol  and  Optimi zat i on ,  vol.  17,  No.  4,  1079. 


13.  Professor  F.  Polak  and  A.  Tits,  a  graduate  student,  have  developed 
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a  new  multiplier  method  which  is  described  in  "A  Globally  Convergent 
Multiplier  Method  with  Automatic  Penalty  Limitation,”  to  appear  in  Jour¬ 
nal  of  Applied  Mathematics  and  Optimization. 

Professor  A.  Sangiovanni  Vincentelli  and  his  students  have  been 
testing  some  of  our  new  optimization  algorithms  to  determine  if  they  can 
be  used  to  improve  significantly  on  a  competent  design.  In  the  case  of 
a  wide-band  amplifier,  designed  by  Prof.  R.  G.  Meyer,  they  have  found 
that  they  could  increase  the  bandwidth  from  900  MHz  to  1300  MHz  while 
meeting  input  impedance,  output  impedance,  insertion  gain,  and  reflec¬ 
tion  coefficient  constraints.  This  is  a  very  significant  improvement 
which  shows  that  our  optimization  algorithms  are  well  suited  to  elec¬ 
tronic  circuit  design. 

15.  In  the  past,  optimization  with  tolerancing  of  practical  electronic 
circuits,  has  been  seriously  limited  by  the  number  of  variables 
involved.  For  example,  a  reasonably  accurate  model  of  a  bipolar  transis¬ 
tor  or  MOS  device  requires  at  least  8  parameter  values  for  its  specifi¬ 
cation.  These  parameters  parameters  are  all  functions  of  bias  condi¬ 
tions,  temperature  and  process  variation.  Since  the  smallest  sub¬ 
circuit  for  which  optimization  can  be  usefully  employed  will  typically 
contain  ten  or  twenty  components,  the  number  of  variables  becomes  prohi¬ 
bitive. 

Prof.  R.  G.  Meyer’s  initial  effort  in  the  area  of  modelling  for 
optimization  has  been  directed  towards  reducing  the  number  of  variables 
required,  while  still  maintaining  model  accuracy.  Initial  work  has  been 
directed  towards  bipolar  integrated-circuit  technology,  where  a  new 
modelling  scheme  has  been  developed.  This  allows  the  modelling  of  arbi- 
tarily  large  circuits  using  only  four  basic  process  parameters  to 
represent  the  effect  of  process  tolerances.  This  advance  has  been  made 
possible  by  recognition  of  the  excellent  tracking  characteristics  which 
are  inherent  in  modern  integrated-circui t  technology  together  with  a 
close  attention  to  the  capabilities  and  requirements  of  the  optimization 
algorithms  being  developed  by  our  research  team.  Using  these  models, 
realistic  optimization  of  complex  circuits,  including  tolerancing,  is 
now  possible. 
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Significant  Accomplishments 

The  significant  accomplishments  of  our  research  can  be  summarized 
as  follows: 

1.  We  have  developed  a  prototype  interactive,  optimization  based  com¬ 
puter  aided  design  package,  capable  of  accepting  semi-infinite  con¬ 
straints  on  frequency  and  time  responses,  but  not  tolerancing  or  post 
manufacturing  tuning.  This  package  serves  as  a  test  bed  for  the 
development  of  ideas  and  concepts  for  a  much  more  sophisticated  package 
which  is  now  being  constructed. 

2.  We  have  shown  that  device  models  can  be  obtained  which  make  it  prac¬ 
tical  to  solve  the  optimal  design  centering,  tolerancing  and  tuning 
problem . 

3.  We  have  developed  a  number  of  optimization  and  inequality  solving 
algorithms  for  computer  aided  design  which  either  solve  problems  that 
were  hetherto  unsolvable,  or  which  converge  much  more  rapidly  than 
currently  available  algorithms. 

4.  We  have  obtained  a  number  of  important  results  in  circuit  and  system 
fault  diagnosis. 
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Scientific  Objective 

The  objective  of  the  research  completed  last  year  was  to  study  arch¬ 
itecture,  design  procedure  and  possible  circuit  implementations  for  very 
large  scale  integrated  (VLSI)  circuits  to  be  used  in  computing  systems  of 
the  future.  Such  chips  will  help  to  construct  digital  computers  that  are 
more  reliable,  less  expensive,  have  higher  performance  and  handle  a  level 
of  complexity  well  beyond  the  capabilities  of  today’s  computers. 

We  used  an  interdisciplinary  team  to  pool  expertise  from  different 
fields  such  as  innovative  integrated  circuit  design,  dynamic,  non-Von-Newmann 
computer  architecture,  software  directed  hardware,  computer  aided  design  and 
optimization.  With  this  team  as  a  whole  or  with  subgroups  of  it,  we  studied 
the  basic  issues  in  the  relationship  between  VLSI  chips  and  future  computing 
systems.  Out  of  this  research,  we  got  an  understanding  for  the  kind  of  chips 
that  will  be  needed  to  build  the  computing  systems  of  the  next  decade.  A 
first  phase  has  been  concluded.  More  detailed  research  into  specific  related 
topics  is  now  being  pursued  under  separate  grants  from  NSF  and  ARPA. 

Abstract  of  Major  Findings  and  Publications 

The  research  during  the  two  years  of  this  grant  focused  on  the  organ¬ 
ization  of  a  large  number  of  future  single-chip  computers  into  a  well-org¬ 
anized  network.  In  particular,  optimal  interconnection  topologies  and  the 
relationship  between  topology  and  applications  have  been  investigated.  A 
second  major  issue  concerned  the  switching  hardware  necessary  to  make  message 
packet  communication  between  these  processors  possible.  The  third  study  area 
concerned  potential  Implementations  of  such  single-chip  computers  in  the  time 
frame  of  the  mid-1980’s.  [Sequin  79] 

Interconnection  topologies  for  single-chip  computers 


Various  interconnection  topologies  based  on  binary  tree  structures  have 
been  studied.  The  binary  tree  is  a  particularly  attractive  interconnection 
topology  in  the  context  of  large  numbers  of  very  large  scale  integration 
(VLSI)  single-chip  computers.  The  average  path  length  between  processors 
grows  only  logarithmically  with  the  number  of  computers  in  the  network, 
and  the  number  of  communication  ports  per  node  is  kept  reasonably  low  as  re¬ 
quired  by  the  limited  number  of  pins  available  on  future  VLSI  chip  packages. 
For  improved  total  bandwidth,  a  more  uniform  message  traffic  distribution 
cud  the  provision  of  redundant  paths  which  is  the  basis  for  fault  tolerance, 
the  basic  binary  tree  needs  to  be  enhanced  with  additional  links  between 
nodes.  Interconnecting  nodes  on  the  same  level  with  horizontal  links  be¬ 
tween  nearest  neighbors  leads  to  half-ring  and  full-ring  structures  which 
have  simple  message  routing  algorithms  and  a  reasonably  good  traffic  distri¬ 
bution  [Sequin  et  al  78]. 

Half-ring  and  full-ring  trees,  however,  have  certain  shortcomings  for 
computational  problems  which  require  a  lot  of  communication  between  distant 
nodes  in  the  tree,  such  as  sorting  of  a  very  large  data  base  distributed 
over  all  leaves  of  the  tree.  A  study  to  find  ways  to  alleviate  these  prob¬ 
lems  has  led  to  a  new  structure,  named  "Hypertree, ”  which  combines  some  of 
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the  best  features  of  the  binary  tree  with  the  compactness  of  an  n-dimension- 
al  hyper-cube.  The  addition  of  suitably  selected  groups  of  n-cube  links  at 
each  level  of  the  binary  tree  provides  short  paths  connecting  those  nodes 
the  address  of  which  differ  in  only  one  bit,  i.e.,  pairs  of  nodes  which 
need  to  exchange  data  rather  frequently  in  algorithms  such  as  sorting  or 
Fast  Fourier  Transforms.  A  paper  to  report  these  findings  has  been  sub¬ 
mitted  to  IEEE  Trans,  on  Computers. [Goodman,  Sequin] 

Dynamic  simulation  of  message  traffic 

The  various  interconnection  topologies  mentioned  above  have  been  analyzed 
analytically  with  respect  to  average  message  path  length  for  certain  simple 
traffic  patterns.  For  the  more  realistic  and  more  complicated  cases  of  messages 
originating  with  some  statistical  properties  and  communicating  with  different 
targets  with  different  probabilities,  e.g.,  because  of  a  dependence  on  lo¬ 
cality,  the  closed  form  analysis  becomes  too  difficult.  For  these  cases  message 
traffic  simulation  for  networks  with  up  to  255  nodes  have  been  performed  on  a 
VAX  11/780.  Even  though  the  simulation  could  run  up  to  a  million  times  slower 
than  expected  real  time  performance  for  the  largest  networks,  significant  in¬ 
sights  could  be  gained  with  respect  to  potential  bottlenecks  in  the  network, 
suitable  routing  algorithms,  the  degree  of  time-multiplexing  desirable  on  each 
link,  and  the  buffering  requirements  in  each  node  [Suhler  and  Despain  79], 
[Goldberg  79]. 

Relationship  between  topology  and  application 

The  performance  of  these  networks  of  single-chip  computers  will  depend  to 
a  significant  degree  on  how  well  a  particular  problem  can  be  mapped  onto  the 
basic  tree  structure.  A  first  problem  being  investigated  is  the  design  of 
very  large  data  base  machines.  Because  massive  amounts  of  data  must  be  moved, 
the  tppology  of  the  processor  interconnections  is  important,  To  determine  an 
appropriate  interconnection  scheme,  a  simple  model  of  a  very  common  but  dif¬ 
ferent  data  base  operation  is  used  to  determine  an  optimal  interconnection 
scheme.  This  is  the  elimination  of  duplicate  information  in  a  collection  of 
data  elements.  In  particular,  four  methods  are  considered  which  can  perform 
the  elimination  of  duplicates.  Two  involve  the  exhaustive  comparison  of  all 
elements  against  each  other,  and  two  involve  global  sorting.  These  methods 
and  their  corresponding  optimal  interconnection  topologies  are  analyzed  and 
compared  to  help  determine  a  suitable  multiprocessor  topology  and  computer 
architecture  for  a  data  base  machine,  A  hybrid  architecture  is  shown  to  be 
near-optimal  under  our  assumptions  [Goodman  and  Despain  80]. 

A  second  invest igat ion  has  been  started  to  determine  how  a  database  such 
as  INGRES  would  fit  onto  the  binary  tree  structure.  Early  results  show,  for 
example,  if  raul ti-relational  queries  are  sent  at  the  rate  of  one  every  five 
seconds,  under  a  reasonable  assumption  for  disks  and  communication  rate,  in¬ 
creasing  the  size  of  the  tree  from  15  to  63  processors  increases  the  through¬ 
put  twenty-fold.  We  have  also  found  that  for  a  tree  of  255  processors,  de¬ 
creasing  the  number  of  disks  from  64  to  16  disks  only  reduces  the  throughput 
twenty  percent.  Static  placement  of  processes  on  processors  gives  better 
performance  than  random  placement  of  processes  on  processors  Jlespain  80], 
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Processor  architecture  issues 

A  time  frame  for  the  potential  implementation  of  the  mentioned  single¬ 
chip  computers  is  the  mid  1980fs,  By  that  time  it  will  be  possible  to  place 
about  one  million  transistors  on  a  single  chip  of  silicon.  While  this  number 
may  seem  large  by  present-day  standards,  it  is  certainly  not  infinite,  and  a 
judicious  allocation  of  these  transistors  has  to  be  made  to  the  various  dif¬ 
ferent  functions  required  in  a  processor  in  order  to  obtain  best  overall  re¬ 
sults.  An  investigation  has  been  carried  out  to  perform  a  tentative  alloca¬ 
tion  of  about  a  million  transistors  in  such  a  manner  that  a  powerful  single- 
chip  computer  results. 

One  of  the  primary  goals  in  this  study  was  to  realize  as  much  primary 
memory  of  short  access  time  on  the  chip  as  possible.  The  solution  was  to 
use  an  on-chip  memory  hierarchy  composed  of  a  substantial  amount  of  high 
density  dynamic  random  access  memory,  buffered  through  high-speed  caches 
built  from  faster  (and  also  large  and  more  power  consuming)  static  memory 
cells.  Separate  caches  are  being  used  for  data,  programs,  microcode,  ad¬ 
dresses  and  descriptors.  This  results  in  a  multiport  memory  which  permits 
effective  sharing  of  the  limited  amount  of  memory  by  these  various  types  of 
information,  resulting  in  increased  flexibility  and  efficiency.  On  the  whole, 
it  turned  out  that  most  of  the  chip  area  is  used  by  memory  functions,  since 
even  t\\e  finite  state  machines  used  for  control  and  other  "random  logic11 
functions  are  best  implemented  in  a  regularly  structured,  compact  manner  by 


Communication  between  single-chip  computers 


The  hardware  required  to  provide  efficient  message  packet  communication 
between  different  nodes  in  the  described  network  has  been  explored.  A  top 
level  description  has  been  outlined  In  [Sequin  et  al  78],  A  prototype  of  a 
communication  switch  has  been  constructed  with  standard  MSI  and  LSI  TTL  parts 
to  verify  these  ideas.  From  the  study  of  this  prototype,  the  crucial  building 
blocks  within  the  coramunicat ion  hardware  were  identified  and  a  suitable  allo¬ 
cation  of  the  available  number  of  transistors  on  future  VLSI  chips  to  the 
various  functions  within  the  communication  switch  was  proposed.  This  has  led 
to  the  definition  of  general  purpose  VLSI  components  useful  to  efficiently 
interconnect  single-chip  computers  of  the  raid-1980?s  [Sequin  80]. 


The  main  emphasis  in  the  design  of  the  interconnection  networks  and  of 
the  switching  hardware  has  been  to  provide  as  much  bandwidth  as  possible  for 
communication  hetween  processors.  Recently,  it  has  become  apparent  that  short 
latency  time  in  establishing  a  message  path  to  another  node  is  of  equal  im¬ 
portance.  In  the  first  design,  a  substantial  Overhead  was  paid  for  the  rout¬ 
ing  of  the  header  of  a  new  message  and  the  setting  up  of  a  new  virtual  com¬ 
munication  channel.  To  reduce  latency  time,  the  most  important  and  time- 
consuming  operations  are  now  performed  with  special  purpose  hardware.  A 
preliminary  analysis  indicates  that  routing  of  a  message  header  through  a 
node  should  ultimately  be  possible  in  less  than  ten  clock  cycles. 
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Support  for  languages  and  systems 

Architectural  support  of  operating  systems  and  languages  is  clearly  an 
important  issue  in  multiprocessors.  One  of  the  drawbacks  of  the  field  has 
been  the  lack  of  a  useful  definition  of  high-level  computer  architectures; 
i.e.,  is  the  Intel  8080  a  high-level  language  computer  architecture  because 
you  can  run  a  high-level  language  program  on  it?  A  review  of  the  salient 
features  of  high-level  computer  language  has  clarified  the  fact  that  hard¬ 
ware  and  software  have  to  be  viewed  together.  This  has  led  to  some  new  in¬ 
sights  of  what  constitutes  a  truly  high-level  computing  system  [Ditzel  and 
Patterson  80]. 

Conclusions 


During  the  two  years  of  this  research,  insight  has  been  gained  into 
some  of  the  relationships  of  the  emerging  VLSI  technology  and  of  mitable 
architectures  for  large  distributed  computing  systems.  While  studies  of 
various  network  topologies  have  found  a  conclusion,  other  issues  are  fur¬ 
ther  pursued  with  separate  funding.  In  practicular,  research  continues 
in  the  area  of  switching  hardware,  implementation  studies  of  a  VLSI  pro¬ 
cessor  with  on-chip  cache  memory,  and  database  application  studies  on 
tree  structured  networks. 

The  two  year  study  concluded  with  this  report  has  been  very  valuable 
to  identify  and  formulate  more  specialized  research  issues  in  the  vast 
field  of  VLSI  computer  architecture. 
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Scientific  Objective 

Image  Processing,  in  particular  Feature  Extraction  and  Image  Coding, 
have  been  studied  extensively  for  at  lea9t  15  -  20  years.  Both  prob¬ 
lem  areas  are  extremely  important.  Feature  extraction  is  the  crux  of 
successful  processing  in  diagnosis  of  medical  images,  analysis  of  re¬ 
connaissance  photos,  or  monitoring  part9  assembly.  Image  encoding  is 
similarly  crucial  in  transmitting  data  back  from  a  satellite  or  recon¬ 
naissance  vehicle  or  in  storing  service  manual  diagrams ,  or  a  wealth  of 
material  In  other  applications.  In  spite  of  the  importance  of  these 
areas,  the  problems  have  not  really  been  solved.  While  clever  ad  hoc 
extraction  and  coding  algorithms  have  been  developed,  it  is  safe  to 
make  two  assertions: 

1)  The  best  of  the  extraction  algorithms  are  not  good  enough  to 
segment  ordinary  grey-scale  images  which  contain  arbitrary  (non-test) 
obj ects. 

2)  In  both  the  case  of  extraction  or  coding,  not  only  is  an 
"optimum"  method  not  known,  but  the  performance  of  the  optimum  ex¬ 
tractor  or  coder  is  not  known  for  comparison  with  computationally 
feasible  methods. 

The  reason  for  this  incomplete  state  of  these  two  problem  areas  is  the 
lack  of  any  sort  of  model  for  the  shape  of  image  objects.  Such  an  ef¬ 
fective  model  would  allow  one  to  derive  optimum  algorithms  and  evalu¬ 
ate  the  effectiveness  of  practical  algorithms  relative  to  the  optimum 
obtainable.  Our  objective  is  to  develop  such  a  model  or  model  with 
the  following  attributes: 

1)  by  appropriate  selection  of  a  modest  number  of  parameters, 

one  should  be  able  to  generate  sample  objects  representa¬ 
tive  of  the  shapes  of  real  objects  in  a  variety  of  con¬ 
texts:  contour  lines,  structures  in  a  reconnaissance 

scene,  or  symbols  in  a  parts  manual. 

2)  the  model  should  allow  one  to  determine  an  optimum  pro¬ 
cessing  (encoding  or  extraction)  algorithm  (or  at  least 
its  structure)  and  obtain  bounds  on  performance.  Further, 
the  model  should  allow  one  to  determine  to  what  extent  ad¬ 
vantage  can  be  taken  of  global  structure;  i.e.,  how  much 
better  does  an  algoritlun  do  that  works  with  the  whole  raster 
than  an  algorithm  that  processes  only  one  or  several  lines 
at  a  time?  This  will  address  the  important  issue  of  the 
relation  between  effectiveness  of  the  processing  and  its 
computational  complexity. 

State  of  the  Art 


What  is  the  current  state  of  trying  to  model  images  and  the  objects  con¬ 
tained  therein?  The  problem  of  representing  two-dimensional  objects  is 
basically  different  than  representing  a  function  of  an  independent  variable. 
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and  the  standard  methods  of  one-dimensional  analysis  and  representation 
do  not  carry  over.  Grenander1  has  produced  a  scholarly  treatise  which 
provides  an  elegant  structure  for  representing  images,  based  on  genera¬ 
tors  or  primative  elements*  However,  this  structure  does  not  really 
solve  the  problem,  since  it  provides  no  insight  or  knowledge  concerning 
what  constitutes  a  reasonable  or  useful  set  of  generators  for  image  ob¬ 
jects.  Without  this,  the  structure  is  useless.  This  choice  of  gener¬ 
ators  must  be  arrived  at  by  a  sense  of  empirical  observations  and  mathe¬ 
matical  gestalt  that  the  community  is  currently  lacking. 

A  number  of  mathematical  models  have  been  proposed  on  a  line-by-line 
basis,  usually  assuming  a  Markov  structure  on  a  TV  raster  scan.  How¬ 
ever,  these  really  beg  the  issue  of  representing  the  two-dimensional 
structure,  and  some  produce  sample  images  which  reflect  this  one-di¬ 
mensional  ordering  in  a  way  that  is  quite  disconsonant  with  the  appear¬ 
ance  of  real  images. 

The  only  really  satisfactory  two-dimensional  model  to  have  been  precise¬ 
ly  formulated  and  extensively  studied  is  the  Tsing  model^*^.  This  study 
was  initiated  hv  physicists  and  has  recently  received  some  attention  by 
mathematicians.  Tills  model  generates  objects  of  maximum  entropy  for  a 
given  average  linear  combination  of  perimieter  and  area.  This  model  is 
Interesting  and  could  in  fact  realistically  represent  natural  objects 
such  as  lakes  or  contour  lines.  However,  it  has  some  drawbacks.  First, 
It  does  not  represent  more  structured  objects  such  as  faces  or  man-made 
objects  that  appear  in  scenes.  Secondly,  it  does  not  lend  itself  to 
calculations;  e.g.,  calculation  of  the  entropy  rate  of  the  Ising  model 
is  possible  only  numerically  and  only  in  certain  sub-cases.  Third,  it 
is  not  a  canonical  representation  of  object  shape:  to  obtain  a  natural 
or  satisfactory  appearing  representation  of  a  smooth  object  requires  a 
very  high  sampling  density  (relative  to  the  curvature  of  the  object). 
Fourth,  generation  of  sample  objects  with  this  distribution  requires 
an  iterative  algorithm  that  uses  extensive  computer  time. 

A  model  represent  ing  objects  as  ruled  surfaces  has  resulted  in  a  useful 
code,  but  this  model  did  not  really  represent  or  exploit  the  two-di¬ 
mensional  structure  of  the  objects  within  the  image'*.  Pavlides^  consi¬ 
dered  a  similar  model  constructed  of  planes  and  deduced  some  interesting 
properties  of  computing  this  representation  from  a  sample  image,  hut  tl^d 
not  address  the  central  problem  of  how  to  model  t  he  shapes  (boundaries) 
of  the  planes. 

Again  we  note  that  many  useful  encoding  and  extraction  algorithms  now 
exist,  such  as  the  Hucckel  edge  detection  algorithm***^.  While  such 
algorithms  are  clever  and  reasonably  effective,  they  do  not  satisfac¬ 
torily  separate  objects  in  real  images,  nor  do  they  answer  the  question 
of  whether  such  separation  is  possible  and  the  computational  cost  re¬ 
quired  to  achieve  it. 

Progress  to  Date 

Our  past  and  proposed  work  focuses  around  three  approaches  which  we 
describe  in  sequence. 


ISS-81-6  -  p.  4 


Two  Component  Coding  of  Grey-Scale  Images 

We  have  developed  a  coding  algorithm^  based  on  an  image  model  which  re¬ 
presents  log | LOG- intensity  across  a  scan  line  as  a  discontinuous  component 
composed  of  straight-line  segments  plus  a  continuous  Gaussian  component 
representing  texture.  This  method  is  quite  effective,  yielding  better 
results  than  other  codes  we  have  seen  and  notably  better  results  than 
a  Transform  Code.  The  method  is  limited  in  that  the  discontinuous  com¬ 
ponent  is  modeled  and  processed  line-by-line.  This  leave  objection¬ 
able  streaks  and  edge  business  and  requires  extra  bits  to  encode  the 
location  of  the  break  points  of  the  line  segments. 

The  algorithm  for  extracting  the  discontinuous  or  object  component  from 
a  grey  scale  image  has  been  modified  to  identify  and  align  breakpoints 
on  successive  scan  lines,  identify  planar  regions  and  approximate  them 
as  such,  and  use  a  nonlinear  estimator  to  smooth  the  breakpoint  location. 
These  improved  descriptions  have  been  used  to  improve  the  two-component 
image  code. 

Tree  Generated  Image  Objects 

We  have  developed  and  studied  a  model  which  generates  objects  by  gener¬ 
ating  a  tree  or  skeleton  and  then  filling  in  smooth  lines  about  the 
’’bones"  of  the  skeleton.  To  date  we  have  generated  sample  objects  from 
this  model;  these  are  similar  in  shape  to  objects  such  as  tools  or  air¬ 
planes.  We  have  also  related  errors  along  the  edge  to  errors  in  the 
generating  points  and  calculated  the  entropy  necessary  to  transmit  the 
objects’  shape  to  within  a  specified  accuracy.  We  are  now  in  the  pro¬ 
cess  of  calculating  the  rate  required  to  describe  these  objects  by  the 
chain-link  code  discussed  by  Freeman®  and  Morin^.  This  work  was  pri¬ 
marily  supported  by  JSEP. 

The  chain-link  code  has  been  investigated  and  better  bounds  were  de¬ 
termined.  Ry  scanning  sample  objects,  algorithms  were  determined  which 
could  be  used  for  improved  encoding  techniques. 

Curvature  Description  of  Object  Shape 

We  have  completed  an  algorithm  that  can  scan  black-and-white  images 
and  determine  a  sequence  of  curvature  pulses  representing  an  object 
boundary.  These  pulses  are  used  to  drive  a  dynamic  model  which  gen¬ 
erates  a  smooth,  pleasing  approximation  to  the  object  boundary.  The 
model  incorporates  feedback  so  that  errors  in  position  and  angle  do 
not  accumulate. 

This  curvature  based  model  has  been  used  in  two  ways.  First,  curv¬ 
ature  sequences  for  a  variety  of  objects  have  been  generated  for  the 
purpose  of  extracting  a  useful  statistical  model  for  these  sequences. 
Second,  the  problem  of  coding  these  sequences  taking  advantage  of 
the  statistical  properties  has  been  examined. 
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Scientific  Objectives 

The  long  term  objective  of  this  research  is  to  achieve  a  dramatic  reduc¬ 
tion  in  the  complexity  of  certain  signal  processing  functions  Implemented  in 
MOS  Integrated  circuit  technology.  In  particular,  analog  sampled  data  and 
switched  capacitor  filters  will  be  used  for  realizing  certain  non-adaptive 
and  adaptive  filters  which  have  previously  been  realized  by  digital  recur¬ 
sive  and  transversal  filter  structures.  To  achieve  this  goal,  the  adaptive 
filtering  algorithms  used  heretofore  will  have  to  he  fundamentally  restruc¬ 
tured.  Specific  attention  will  be  directed  at  applications  in  speech  pro¬ 
cessing,  including  analysis,  synthesis,  and  recognition  of  speech. 

State  of  the  Art 


This  proposal  touches  on  several  diverse  areas  of  past  research,  includ¬ 
ing  algorithms  for  the  analysis,  synthesis,  and  recognition  of  speech,  as  well 
as  the  implementation  of  signal  processing  algorithms. 

Methods  of  synthesis  of  speech  are  closely  tied  to  the  corresponding  analysis 
method,  so  we  concentrate  on  the  latter  here.  Speech  analysis  has  been  at¬ 
tempted  by  representation  with  suitable  basis  functions^ and  by  short-time 
spectral  analysis  techniques'”-^.  The  latter,  which  have  been  more  success¬ 
ful  and  much  more  extensively  pursued,  include  the  channel  vocoder^”^ ► ^ ,  es¬ 
timation  of  format  frequencies  by  Fourier  transform  techniques^ > ^ , 6 ^  and  the 
homomorphic  vocoder^ The  channel  vocoder  is  usually,  but  not  always, 
implemented  by  a  bank  of  bandpass  filters,  while  the  remaining  techniques  are 
almost  always  implemented  digitally.  Gold^  gives  an  excellent  summary  of  the 
current  state  of  the  art  in  digitally  implemented  vocoders.  A  typical  complex¬ 
ity  for  an  implementation  with  general  purpose  devices  is  about  150  integrated 
circuits,  with  a  speed  of  operation  consistent  with  only  high  speed  bipolar 
logic  families. 

In  LPC,  which  has  been  extensively  Investigated  because  of  Its  very  ef¬ 
ficient  representation  of  the  vocal  tract  parameters,  It  Is  most  common  to  de¬ 
rive  the  predictor  coefficients  or  reflection  coefficients  on  a  block  basis 
using  recursive  algorithms® However,  it  is  also  possible  to  update  pre¬ 
dictor  coefficients  continuously  in  time  by  an  adaptive  gradient  algorithm 
similar  to  that  used  in  adaptive  noise  cancelling^. 

A  typical  speech  recognition  system  is  con^osed  of  three  separate  ^no¬ 
tional  blocks^®*^ » The  first  block  is  the  signal  processing  front  end 
which  attempts  to  represent  the  speech  in  terms  of  a  minimal  set  of  para¬ 
meters,  such  as  the  LPC  coefficients  discussed  earlier.  Second,  a  measure 
of  the  similarity  between  the  speech  sound  being  analyzed  and  a  set  of  tem¬ 
plates  is  calculated.  Finally,  a  matching  technique  is  used  to  determine 
the  template  which  best  fits  the  speech  sound. 


The  two  parts  of  a  speech  recognition  system  which  are  most  conplex  are 
the  front  end  speech  processing  and  the  template  matching  . 


Previous  Work  at  Berkeley 

At  Berkeley  a  variety  of  analog  integrated  circuits  have  been  developed 
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to  perform  specific  signal  processing  functions.  These  circuits  include  A/D 
converters,  multiplying  D/A  converters,  and  a  variety  of  approaches  to  mono¬ 
lithic  MOS  filtering.  The  most  successful  development  in  the  latter  category  * 
has  been  the  switched  capacitor  filter  technique,  with  which  low  pass,  band 
pass,  and  high  pass  filters  can  be  implemented.  These  filters  require  a 
small  silicon  area,  and  thus  can  be  fabricated  on  the  same  integrated  circuit  ♦ 
with  other  analog  and  digital  circuitry. 

In  addition,  at  this  time  there  is  an  on-going  effort  in  applications  of 
MOS-LSI  to  low  cost  implementations  of  narrow  band  vocoding.  This  research 
is  concerned  with  applying  the  new  MOS  analog  sampled  data  techniques  to  the 
implementation  of  existing  speech  processing  algorithms21.  This  effort  will 
complement  the  research  proposed  here,  which  is  concerned  with  developing  en¬ 
tirely  new  algorithms  to  more  fully  exploit  the  capabilities  of  MOS-LSI,  par¬ 
ticularly  as  it  is  expected  to  evolve  into  VLSI  in  the  future. 

Progress  and  Publlcat icms 

A.  Under  JSEP 

A  new  generalized  form  of  the  lattice  filter  in  which  the  delay  element 
’  has  been  replaced  by  an  arbitrary  all-pass  filter  has  been  discovered.  This 

has  been  written  up  in  the  following  paper,  which  has  been  accepted  for  pub¬ 
lication:  D.  G.  Messerschmitt ,  “A  Class  of  Generalized  Lattice  Filters,” 

1EE  Trans.  AccousticsT  Speech  and  Signal  Processing,  April  1980. 

Abs t  rac  t 

It  is  shown  that  the  digital  lattice  filter  utilized  as 
a  predictor  in  speech  processing  is  a  special  case  of  a 
broad  class  of  filters  with  identical  properties.  In 
which  the  delay  element  is  replaced  by  an  arbitrary  an¬ 
alog  or  digital  all-pass  filter.  The  performance  of 
three  filters  in  this  class,  the  delay  element  and 
single-pole  analog  and  digital  all-pass  sections,  are 
compared  for  the  spectral  estimation  application.  An¬ 
other  application  is  a  realization  of  analog  transfer 
functions  in  a  lattice  realization  using  a  switched  cap¬ 
acitor  all-pass  section. 

•  ^  significant  issue  in  the  choice  of  an  all-pass  filter  for  the 

generalized  lattice  is  the  effect  on  filter  coefficient  sensitivity. 

' t  "is  senslt ivity  relates  directly  to  the  hit  rate  required  when  filt- 

:  er  coefficients  are  used  to  characterize  the  vocal  tract  in  speech 

j  applications.  I he  first  step  was  to  obtain  for  the  first  time  ana- 

>|  vticai  expressions  for  t  lie  sensitivity,  which  have  been  published 

In  the  following  paper:  Peter  I,.  Chu  and  David  C.  Messerschmitt, 

Zero  Sensitivity  of  the  Digital  Lattice  Filter,”  1  KICK  Int.  Con  f .  on 
Audio,  Acoustics  and  Signal  Processing,  April  1980. 
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Abstract 

Recursive  equations  for  the  first-order  zero  move¬ 
ment  sensitivities  of  the  all-zero  digital  lattice 
filter  are  given.  In  addition  to  being  useful  for 
numerical  evaluation,  these  equations  allow  general 
conclusions  to  be  drawn  concerning  the  quite  differ¬ 
ent  sensitivity  properties  of  early  and  final  stages. 

The  approach  differs  from  previous  sensitivity  stud¬ 
ies,  which  considered  measures  consisting  of  integ¬ 
rals  of  the  ratio  of  the  deviated  spectrum  to  the 
original  spectrum,  and  is  therefore  more  useful  in 
directly  studying  the  sensitivity  of  format  band¬ 
width  and  center  frequencies  when  the  filter  is  used 
in  speech  synthesis. 

It  is  planned  to  publish  an  expanded  version  of  this  paper  in  the 
ASSP  Transactions  in  the  near  future. 

More  recently,  an  important  new  application  of  the  generalized  lat¬ 
tice  has  been  found.  It  has  been  found  that  by  choosing  the  all-pass 
filter  judiciously,  the  sensitivity  characteristics  can  be  tailored 
closely  to  the  human  ear.  Specifically,  the  low-frequency  poles  can 
be  made  less  sensitive  to  filter  coefficient  changes,  and  the  high 
frequencey  poles  correspondingly  more  sensitivie,  with  an  attendant 
subjective  improvement.  The  result  should  be  a  reduction  in  the  bit 
rate  required  for  LPC  speech  comparable  to  that  previously  obtained 
by  use  of  pre-emphasis. 

Current  work  in  this  area  is  concentrating  on  proper  choice  of  an 
all-pass  filter,  subjective  tests  to  confirm  the  benefit  derived  from 
the  generalized  lattice,  and  further  understanding  of  the  sensitivity 
issues.  In  the  latter  area,  effort  is  being  devoted  to  understanding 
of  the  sensitivity  properties  of  other  sets  of  parameters  useful  for 
representing  the  vocal  tract  conf igurabion,  such  as  autocorrelation 
coefficients. 

The  second  major  area  of  study  has  been  application  of  these  tech¬ 
niques  to  a  speech  recognition  system,  as  well  as  other  components  of 
such  a  system  including  template  matching. 

Although  template  matching  has  been  done  on  high  speed  digital  com¬ 
puters  in  the  past,  it  appears  to  be  particularly  advantageous  to  imple¬ 
ment  these  time-consuming  and  regular  algorithms  in  hardware.  One  part¬ 
icularly  difficult  aspect  of  template  matching  is  the  time  alignment  of 
the  analysis  information.  A  dynamic  programing  algorithm  has  been  in¬ 
vestigated  for  use  in  providing  an  optimal  method  of  alignment  and  time 
warping.  Extensive  computer  simulations  were  made  to  determine  the 
performance  of  the  algorithm.  Accuracies  which  were  achieved  in  these 
simulations  were  from  98-99%  for  a  vocabulary  of  50  words. 
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Progress  and  PuM  i  oat  i  ons  ,  cont'd 

A  breadboard  ti-  simulate  an  integrated  circuit  which  would  implement 
che  dynamic  programing  algorithm  in  real  time  was  then  constructed  and 
its  poiiormanee  amaars  consistent  with  the  computer  simulations.  Care¬ 
ful  .  ons ider at  ion  oi  the  archit ccture  has  made  it  possible  to  consider 
implementing  tin*  entire  dynamic  programing  function  on  a  single  chip. 
Preliminarv  results  have  been  reported  in  the  following  paper: 

H.  Murvelf  .  M.  howrv  and  R.  W.  Brodersen,  "how  Cost  implement  at  ion  of 
speech  Recognition  System, M  Proc .  of  COMECON,  San  Francisco*  March  l‘)80. 

H.  Coder  Other  _Spo:  orshiy 

I'ndor  NSF  spon.se:  'p,  wv  are  studying  the  applcat  ion  of  the  lattice 

t  liter  to  I. PC  and  to  t  ne  prediction  filter  in  AiMVM.  01  particular  ivn- 
corn  is  the  speed  of  adapt  at  Lon  of  adaptive  forms  of  the  lattice  filter* 
and  unde  rs  t  and  i  ng  of  this  I -sue  is  being  developed. 

Also  under  NS!  anc  ! nc.i.st  r  i  al  support ,  we  are  exploring  various 
hardwurd  implement  at  ion  issues  relating  to  adaptive  filters.  Recently, 
breadboard  ing  has  heei  a  new  approach  to  the  implementation  of  an  ana¬ 
log,  adaptive  transversal  filter  and  an  AD  PCM  speech  coder  with  adaptive 
tpiani  i:*  i at  ion  which  uses  svi t  ched-eapacitor  techniques. 

Additional  work  oa  imp :  ement  at  ion  of  adaptive  filters  is  being  per¬ 
formed  under  a,  cord  :  .1,  *  tf  DAKPA  for  implementation  of  a  low  cost  vocoder 

for  socute  speech  ;  ssi  on .  Idlest*  Implementations  use  conventional 

algorithms  Inst  e  :u  ei  '  new  generalized  lalt  Ice  1  liters  described  above 

S ign it  i i_.\n t  Ac c ore 1 1  i  shmen t  s 

1.  Discovery  u'  the  generalized  form  of  the  lattice  filter,  of  in¬ 
terest  particularly  in  swi t ched-capac itor  implementation. 

2 .  The  real!:  .:;  iv'n  that  the  lattice  filter  can  he  used  to  realize 
arbitrary  eont  iniions -*  i,v,e  as  well  as  discrete-time  transfer  functions. 

1.  \  simplified  proof  of  lattice  filter  stability  properties. 

4.  I  r.ders  t  and  i  ng  of  the  sensitivity  properties  of  the  lattice  filter 

5.  Extension  of  "his  sensitivity  understanding  to  the  generalized 

I  at  t ice  filter. 

b.  Discoverv  that  the*  generalized  lattice  can  he  used  to  reduce 
hit  rate  of  RFC  encoded  speech  for  a  fixed  subjective  quality. 

Development  of  a  hardware  approach  to  dynamic  programing  t  l me 
alignment  in  speech  reeognlt  ion  template  matching. 


